seria BAZELE MATEMATICE ALE CERCETĂRII OPERAŢIONALE APPLIED QUEUEING THEORY ALEC M LEE, M A Director of Operaţional Research, Air Canada The Macmillan Press Limited London © Alee M I ee 1966 Alee M Lee TEORIA AŞTEPTĂRII CU APLICAŢII Traducere din limba engleză de dr VINCENŢIU DUMITRU, ANDREI ALEXIU ţi TIBERIU SPIRCU Editura Bucureţti - 1976 PREFAŢA LA EDIŢIA ROMÂNĂ Cu toate că n-am vizitat România pînă în prezent, îmi imaginez că sistemele de aşteptare sînt la fel ca oriunde; ca şi canadienii, germanii sau japonezii, probabil că românii sînt la fel de iritaţi de aspectele neplăcute ale sistemelor de aşteptare De aceea, voi fi foarte mulţumit dacă cele scrise în cartea mea vor ajuta la reducerea inconvenientelor inerente sistemelor de aşteptare ALEG LEE ALEG LEEîn ultimii ani articolele şi cărţile avînd ca temă teoria matematică a aşteptării s-au înmulţit în mod alarmant Totuşi puţine sînt acelea care au dovedit că au o utilitate practică; de aceea n-am considerat necesar să mai adaug referinţe bibliografice Trebuie însă să scot în evidenţă faptul că în anii care au trecut din momentul scrierii cărţii, tehnicile de simulare au devenit mult mai importante şi mai folositoare decît am presupus atunci Aceasta se datoreşte doar în parte progresului acestor tehnici în sine, fiind mai mult o consecinţă a tendinţei foarte favorabile a costului calculelor, ea urmare a progresului rapid al tehnologiei echipamentelor de calcul România, care a făcut din proiectarea şi fabricarea echipamentelor de calcul o componentă semnificativă a dezvoltării ei tehnologice, este bine plasată pentru a fi beneficiara aeestor tendinţe Ednaston, Derbyahlre Februarie 1976 PREFAŢĂ Revistele de cercetare operaţională din toată lumea continuă să fie pline cu lucrări legate de teoria firelor de aşteptare, iar un cititor ocazional ar putea să-şi formeze părerea că problemele acestui subiect au fost deja epuizate Nu este însă aşa Majoritatea acestor lucrări îşi au originea, se pare, în zona crepusculară a cercetării academice în ele sînt explorate curajos misterele modelelor simple ale celor mai familiare procese de aşteptare şi sînt obţinute binecunoscutele proprietăţi ale modelelor clasice î n timp ce aparatul matematic devine din ce în ce mai bine pus la punct, rămîne într-adevăr o dificultate a găsi relatări ale muncii experimentale cu procesele de aşteptare, sau a face observaţii empirice ale acestor procese, sau chiar aplicaţii ale unor teorii existente Se pare că o mare parte a muncii de cercetare se depune pentru teorie şi o mult mai mică parte pentru aplicaţii Există o aparentă inactivitate în aplicaţiile teoriei aşteptării, dar realitatea este alta Problemele legate de aglomeraţie constituie un sector majorîn care este realizată foarte multă practică de cercetare operaţională, par pentru a o descoperi este necesar — exceptînd sectoarele de transport şi telecomunicaţii — a asculta conversaţiile de informare necontrolate, care au loc în relaxantele ore ce urmează conferinţelor Sînt mai multe motive pentru această reticenţă, dar unul este dominant în rîndurile multor analişti de cercetare operaţională şi anume că profesiunea lor este fără îndoială de natură matematică, cu un foarte ridicat nivel de rigoare Ei ştiu că în propriile lor aplicaţii cotidiene ale teoriei aşteptării asemenea rigoare nu este menţinută, pentru că de obicei nu este posibil a formula modele ale proceselor de aşteptare în termeni matematici foarte riguroşi şi a le rezolva De aceea, aceşti analişti sînt refractari în a discuta experienţele lor pentru teama că ele vor fi expuse criticilor unor persoane competente Dar această teamă nu este justificată Mulţi dintre cei mai competenţi specialişti în teoria aşteptării sînt intens preocupaţi de aplicaţiile practice, chiar dnd aceste aplicaţii implică o mare ramificaţie şi diversificare a modelelor teoretice Mai mult, această atitudine de reticenţă îi separă pe studenţii în cercetare operaţională de realităţile fenomenelor pe care le analizează şi duce la consolidarea ideii că modelele teoretice sînt altceva decît aproximări conceptuale ale realităţii Lucrarea de faţă este o încercare de a abandona această reticenţă Ea nu tratează despre teoria firelor de aşteptare, ci despre problemele aplicării rezultatelor ei în practică Sînt multe cărţi adecvate, şi chiar cîteva excelente, legate de această teorie şi nu am avut niciodată nici interesul, nici competenţa să le completez numărul Altfel spus, scopul acestei cărţi este aplicarea teoriei aşteptării în găsirea de soluţii operaţionale problemelor de aşteptare din lumea reală Deoarece este dificil de găsit cazuri concrete în reviste (fac excepţie unul sau două din articolele apărute mai recent cum ar fi „clasicul" „Traffic Delays and Toll- -booths" al lui Edie), am fost nevoit să utilizez materiale din surse familiare sau de la prieteni, dintre care unele nu au fost pînă acum prezentate publicului larg Aceasta a însemnat să folosesc material obţinui din studii la care am fost eu însumi participant Am urmat această cale cu regret însă pot fi sigur de faptul că toate cazurile prezentate sînt reale, ca şi datele utilizate şi — cu anumite excepţii specificate în text — nefalsificate Pericolele falsificărilor sînt limpezi Nu este o carte cuprinzătoare Unele modele de aşteptare ce au fost foarte mult studiate nu sînt menţionate sau sînt prezentate doar în cîteva rînduri Am punctat numai ceea ce cred necesar pentru atingerea scopului propus Mai mult, nu există digresiuni legate de teoria probabilităţilor, statistica matematică sau de transformata Laplace Se presupune că cititorul are anumite cunoştinţe despre aceste subiecte Orice cititor care este familiarizat cu vreuna dintre cărţile teoretice mai bune, legate de aşteptare (cum ar fi cartea Elements of Queueing Theory a lui Saaty) poate omite capitolele 2 — 5 ale prezentei cărţi Aceste capitole au fost scrise pentru cei ce se iniţiază în acest domeniu sau pentru cei a căror memorie necesită o reîmprospătare printr-o recapitulare rapidă Oricum, sper ca această carte luată global să fie interesantă şi utilă Sînt profund îndatorat multora pentru ajutorul acordat 0 recunoaştere în acest sens este făcută în cadrul textului şi notelor capitolelor Trebuie oricum să-mi arăt recunoştinţa şi să aduc mulţumirile mele Companiei British European Airways, precum şi Companiei Air Canada pentru permisiunea de a utiliza extensiv materiale legate de studiile de cercetare operaţională la care am participat în perioada 1956 — 1964 Fără acest ajutor n-aş fi putut nici măcar să încep să scriu această carte Aş dori de asemenea să mulţumesc secretarei mele, d-nei Jeannine Lacelle pentru magistrala transpunere a indescifrabilelor mele notiţe; de asemenea soţiei mele, Barbara, pentru realizarea celor mai multe dintre grafice, a indexului şi a unei mari cantităţi de cafea ALEC M LEE Montreal Place Vllle-Marie INTRODUCERE Pe măsură ce societăţile umane devin din ce în ce mai complexe, se pare că devine tot mai dificil pentru membrii lor de a evita participarea involuntară la procesele de aşteptare Fiecare cititor al acestei cărţi trebuie să aibă o oarecare experienţă personală în acest domeniu Pentru a cumpăra un timbru de la un oficiu poştal este necesar să stai la rînd pentru a fi servit La aeroport, pentru a primi confirmarea unui zbor, de cele mai multe ori eşti nevoit să stai tot la rînd Dimineaţa, înainte de a pleca de acasă, pentru a primi prin telefon date privind starea vremii, eşti obligat să formezi de mai multe ori numărul, pînă a obţine ceea ce doreşti Toate acestea sînt procese de aşteptare Experienţa legată de aşteptare este foarte bogată Dacă experienţa pasivă a unei activităţi pe parcursul unei perioade mai lungi ar duce la calificativul de expert (aşa cum mulţipracticieni insistă că este), atunci noi toţi am fi experţi în fenomene de aşteptare De fapt, proporţia oamenilor ce pot fi găsiţi stînd la rînd la orice oficiu poştal şi care sînt experţi în materie de sisteme de aşteptare este redusă : şi s-ar putea admite, dar n-ar fi prea plăcut, că proporţia lor printre angajaţii oficiului poştal este încă şi mai mică Nu poate fi presupusă prea multă familiaritate în legătură cu caracteristicile sistemelor de aşteptare De aceea, ar fi nimerit a începe prin exemplificare, prin a stabili ceea ce înseamnă un sistem de aşteptare şi ce probleme se ridică cel mai semnificativ şi mai des în legătură cu asemenea sisteme Deşi cele trei exemple ce urmează sînt descrise umoristic, ele sînt reprezentative pentru problemele şi procesele foarte reale „Je me presse de rire de tout, de peur d’etre obligi d’en pleurer” Ele prezintă întîmplări frecvente 1 1 pr oblema oţel Ariilor Corporaţia „The Allied Coke and Iron” este un vast concern ce acoperă diverse activităţi economice Printre componentele sale sînt mai multe întreprinderi de colectare a fierului vechi, precum şi două oţelării Cea mai veche oţelărie foloseşte procedeul Siemens de obţinere a oţelului, avînd o nevoie constantă de fier vechi de bună calitate Fierul vechi se transportă în camioanele întreprinderilor de colectare ale companiei Camioanele sosesc la oţelărie şi, înainte de a merge la furnale ca să-şi descarce încărcătura, trec peste un cîntar-basculă, pentru a o cîntări şi înregistra Acesta este punctul critic La aproape orice oră, între 6 dimineaţa şi 10 seara, la cîntarul- basculă este un şir (rînd) format din camioane ce aşteaptă să fie cîntărite Lungimea şirului variază; dar aproape fiecare camion pierde o parte apreciabilă a timpului aşteptînd Unul dintre operatorii de la cîntarul-basculă are o idee : oare n-ar fi mai bine pentru companie să instaleze un al doilea cîntar-basculă şi astfel să reducă timpul de aşteptare al camioanelor? Oare costul instalării şi folosirii cintarului- basculă suplimentar nu va compensa posibilitatea de a spori utilizarea camioanelor? Operatorul (pe nume Szymanowski) semnalează în scris sugestia sa şi o înaintează conducerii companiei După cîteva luni sugestia sa ajunge să fie examinată de către Comitetul director Inginerul-şef adjunct este însărcinat să facă o evaluare sumară El replică imediat că : „este clar că prin dublarea capacităţii de cîntărire timpul mediu de aşteptare a fiecărui camion se va înjumătăţi Aceasta ne va aduce un cîştig relativ mic Eu cred că propunerea serviciului tehnic, de a ne vinde camioanele şi de a cumpăra unele mai mari, este mai valoroasă Ea ne va costa doar 250 000 $ Ceilalţi membri ai Comitetului director (puţin obosiţi) sînt impresionaţi de acest raţionament critic Domnului Szymanowski i se mulţumeşte pentru sugestia sa care, regretăm, este imposibil de materializat astăzi” şi i se acordă o primă de consolare de 20 $ PROBLEMA CASEI DE BILETE Acţiunea are loc într-o gară din micuţul oraş scoţian Paltries, într-o frumoasă dimineaţă de iulie Tocmai s-a oprit ploaia Murdo Macdonald intră să-şi cumpere un bilet pînă la Edinburgh, pentru trenul de la ora nouă El găseşte deschis, ca de obicei, un singur ghişeu destinat vînzării biletelor de clasa întîi (pentru care nu este nici un amator) şi un altul pentru clasa a doua (unde se află cincisprezece persoane aliniate aşteptînd să le vină rîndul) Macdonald se aşază la coadă pentru clasa a doua El este un om rezonabil şi, în ocaziile similare din trecut, a rezistat tentaţiei de a reclama nedreptăţile din organizarea căilor ferate Chiar şi astăzi el nu sugerează ca ambele case de bilete — cea pentru clasa întîi şi cea pentru clasa a doua — să vîndă acelaşi gen de bilete şi ca pasagerii pentru clasa întîi să fie constrînşi astfel să stea la rînd cu ceilalţi El îşi dă seama că un pasager ce va călători cu clasa întîi plăteşte mai mult şi este normal să i se ofere un serviciu superior calitativ, dar astăzi el este nespus de iritat de mişcarea înceată a funcţionarului de la ghişeu şi, aşa cum se în- tîmplă de cele mai multe ori, din iritarea sa îi vine o idee De ce să nu se vîndă biletele pentru ambele clase la amîndouă case de bilete, dar astfel încît să i se îngăduie fiecărui pasager ce doreşte bilet pentru clasa întîi să intre în faţa rîndului, de îndată ce a sosit? în acest fel privilegiile nu vor fi diminuate, în timp ce capacitatea de vînzare a biletelor va fi practic dublată (pentru că de fapt nu există pasageri pentru clasa întîi la Paltries) „Adoptînd această procedură” scria mai tîrziu Macdonald către editorul revistei „Paltries gazette and Stock- breeder’s Tribune”, „timpul mediu pierdut la rînd va fi redus în ceea ce mă priveşte, ca şi în ceea ce priveşte şi pe ceilalţi pasageri de clasa a doua, şi aceasta fără o cheltuială suplimentară a societăţii de căi ferate şi fără vreun inconvenient pentru cei doritori de bilete de clasa întîi” PROBLEMA SOŢIEI FURIOASE William K Overlander se găseşte în barul clubului de golf fellowridge, sorbind un Martini şi adresîndu-se unui coleg de club Sste indignat El tocmai a emis o instrucţiune în cadrul companiei ale, al cărei efect va fi că în viitor nici unul dintre angajaţi (dar abso- ît nici unul) nu va mai utiliza serviciile Companiei aeriene Transce- ental Airlines pentru călătoriile oficiale de afaceri; aceasta din auză că în ziua anterioară soţia lui a avut o experienţă teribilă „După cum îţi spuneam”, povesteşte el printre înghiţituri, „Marilyn a format de patru ori, în interval de 30 de minute, numărul de telefon al biroului lor de rezervări de bilete şi nu a reuşit să-l prindă Apoi, în sfîrşit, a cincea oară a reuşit Mi-a povestit că a sunat de trei ori şi apoi a răspuns o voce ce părea mai curînd o înregistrare pe bandă, care i-a spus că agenţii sînt foarte ocupaţi şi să aştepte puţin Ei bine, a aşteptat circa jumătate de minut, timp în care nu s-a petrecut nimic, dar absolut nimic, cu excepţia intervenţiei vocii înregistrate ce a repetat mesajul anterior Ca urmare, a închis” „Ce a urmat ? ” — îl întreabă prietenul său „A încercat din nou ? ’ ’ „Nici vorbă”, zise Overlander „Aceasta nu a fost prima dată cînd se întâmplă aşa ceva la Transcedental Nu, domnule, soţia mea a sunat la Compania aeriană Commonwealth şi Caribbean Airways şi ştii ce s-a întâmplat? A dat de ei de prima oară, fără nici un fel de aşteptare Iată ceea ce se numeşte într-adevăr un bun serviciu Băieţii de la Transcedental trebuie să capete noi ajutoare, după părerea mea Probabil ei au numai o singură bandă înregistrată şi trei femei bătrîne la panoul de comandă ” Overlander, ale cărui sentimente le înţelegem şi le respectăm, ar fi surprins să afle că Transcedental are angajaţi exact acelaşi număr de agenţi pentru biroul de rezervări de bilete, ca şi Commonwealth şi Caribbean, pentru a servi (foarte convingător şi coincidenţial din punctul nostru de vedere) exact acelaşi număr de pasageri Nu este vorba despre un număr prea mic de funcţionari sau de un refuz de a investi bani; nu acestea sînt cauzele pentru care cei de la Transcedental sînt mai puţin prompţi în efectuarea serviciilor Este vorba, după cum se va vedea, de ignoranţă NATURA SISTEMELOR DE AŞTEPTARE Aceste trei povestioare moraliste sînt toate despre procese d« aşteptare şi ilustrează modul în care se pun problemele legate dt astfel de procese în toate trei există clienţi care cer să fie serviţi Camioanele cu fier vechi trebuie să fie cîntărite; cetăţeanul Macdo nald vrea să i se vîndă un bilet de tren; stimabila cetăţeancă Over lander vrea să vorbească cu un agent al unei companii aviatice pentn a i se rezerva im bilet de avion Clienţii doresc să fie serviţi imediat, dar dorinţele lor nu sînt satisfăcute imediat, ei trebuind să aştepte la rînd Momentele în care ei se aşază la rînd se numesc momente de sosire Intervalul dintre momentul sosirii unui client şi momentul în care el este servit se numeşte timp de aşteptare Clienţii sînt serviţi de către staţii de servire, numite uneori şi canale sau ghişee Termenii : staţie de servire, canal sau ghişeu semnifică acelaşi lucru şi în această carte vor fi folosiţi alternativ Cîntarul-basculă din primul exemplu este un ghişeu; funcţionarii de la ghişeu din al doilea exemplu reprezintă staţii de servire; fiecare telefon de la biroul pentru rezervări de bilete pentru avioane este un canal Fiecare termen îşi are propria valoare ilustrativă Timpul scurs pentru satisfacerea serviciului cerut de client se numeşte timp de servire Suma dintre timpul de aşteptare al clientului şi timpul de servire este îndeobşte cunoscută ca fiind timpul total de aşteptare De cele mai multe ori este mai puţin important să ne gîn- dim la timpii totali de aşteptare, ci mai degrabă la timpii de aşteptare la rînd, dar apar frecvent şi excepţii Tipurile de întrebări ce se pun de obicei, în legătură cu sistemele de aşteptare, sînt următoarele : Cîte staţii de servire trebuie instalate suplimentar, dacă numărul clienţilor şi timpul pentru servirea fiecărui client nu se schimbă, pentru a reduce procentul clienţilor care au de aşteptat mai mult de 2 min la mai puţin de 5 % ? Ce efect are introducerea priorităţilor pentru anumite clase de clienţi asupra timpului de aşteptare? Cit de mult trebuie mărit un fir de aşteptare pentru a reduce pierderea clienţilor la un procent mai mic decît 10 %? Ce ar fi mai convenabil: să se instaleze o altă staţie de servire sau să se mărească camera de aşteptare ? Dacă există două clase de clienţi, o anumită calitate a servirii clienţilor (măsurată prin timpul mediu de aşteptare) poate fi realizată mai ieftin prin înzestrarea fiecărei clase cu cîte un grup separat de servire sau prin introducerea unui singur grup care să servească ambele clase ? Acestea sînt tipurile de probleme abordate şi tratate în cadrul acestei cărţi MODELE ALE SISTEMELOR DE AŞTEPTARE GENERALITĂŢI Nu există o ştiinţă a sistemelor de aşteptare Există o teorie a aşteptării legată de studiul unor modele matematice riguroase, dar care dă impresia că se bazează pe un surprinzător de şubred postament de cunoştinţe empirice Nu sînt multe rezultatele experimentale obţinute în legătură cu procesele de aşteptare Aşa după cum se va vedea, acest fapt bizar — din punct de vedere ştiinţific — are influenţe negative privind aplicabilitatea şi utilitatea rezultatelor obţinute pe baza teoriei Dar înainte de a studia acest aspect deosebit de important, precum şi consecinţele sale, va trebui să dăm un răspuns la două întrebări preliminare : ce este un model simbolic — sau matematic*’ şi care este scopul unui astfel de model? Pentru un model al sistemului de aşteptare din staţia de cale ferată Paltries putem considera următorii factori: Particularităţile funcţionarului de serviciu (James Cameron este conştiincios, dar lent; Angus Sutherland este un tip iute, dar înclinat să dea bilete greşite sau să nu dea corect restul; Alan Robertson este rapid şi eficient; Jan Macrae nu este deloc practic) Precizarea zilei din săptămînă (sînt mai mulţi pasageri de luni pînă vineri decît sîmbetele sau duminicile) ♦> Churchinan, C W , Akcoff ‘ R , Arnoff, L Introduction to Operations Research- New York, Wiley, 1957 Starea vremii — dacă a fost foarte frig în timpul nopţii sau nu (variaţia de temperatură de la Robbinston la Paltries este foarte mare şi dacă rouă a îngheţat pe linia ferată, trenul va întîrzia ; nu va mai exista graba în acest caz în procesul de aşteptare) Prezenţa sau absenţa călătoarei Fiona McTaggart (dacă ea se va afla la coadă pentru un bilet, atunci operaţia va întîrzia îngrozitor de mult — pînă la disperare — pentru că : a) nu-şi găseşte niciodată portofelul pentru a-şi putea plăti biletul; b) parlamentează cu fiecare funcţionar de la ghişeu despre o reducere de preţ al biletului, care de fapt nu poate fi obţinută decît joia în iulie, după ora 17) Natura călătoriei pasagerilor (Majoritatea pasagerilor merg la Edinburgh în partea de sus şi la Robbiston în partea de jos, aceştia necreînd nici o problemă Dar, întrucît nu există bilete tipărite pentru staţiile de după Berwick-on-Tweed, în cazul cînd un pasager doreşte totuşi un astfel de bilet, acesta trebuie scris de mînă pe un formular special, apoi înregistrat într-un registru negru Aceasta, evident, necesită timp ) Temperatura şi gradul de umiditate existent (Timpul de servire a unui pasager creşte proporţional cu ambele elemente, exceptînd cazul în care funcţionarul de serviciu este Alan Robertson care s-a născut la Madras şi se simte în elementul său atunci cînd este cald ) Pentru procesele de aşteptare în care sînt implicaţi clienţii-oameni şi oameni ce realizează serviciul cerut, oricine care a avut un contact cu astfel de probleme un timp suficient de lung poate alcătui o listă asemănătoare, cuprinzînd un şir de factori de acest gen Evident, nu toţi factorii au aceeaşi importanţă Cititorul va intui că numărul 4 — McTaggart — este cu mult mai puţin important decît numerele 1, 2 şi 6 Dacă elementul esenţial al procesului îl constituie aglomeraţia din orele de vîrf, atunci factorul 3 va fi semnificativ doar dacă ora de vîrf este în primele ore ale dimineţii sau momentul de vîrf este iarna; în celelalte cazuri nu Raţionînd în acest mod, cititorul poate scurta lista factorilor ce trebuie luaţi în considerare în orice moment al procesului Să ne oprim puţin asupra dezamăgirilor d-lui Karel Szymanowski Factorii ce influenţează operaţia de cîntărire pot fi următorii: L° frecvenţa (intensitatea) sosirilor camioanelor; 2°particularităţile jperatorului de la cîntar; 3° încărcarea medie a unui camion (camioa- îele excesiv încărcate se vor deplasa mai lent şi vor pierde astfel nai mult timp în deplasarea în rînd spre cîntar); 4° camioanelor de zece tone le este permis să se aşeze la rînd în faţa vagoanelor de cinci tone dacă şi numai dacă ele transportă deşeuri Etc După cum se poate remarca, aceşti factori nu sînt aceiaşi cu cei pe care i-am găsit în ilustrarea anterioară, nici nu putem presupune că ei interacţionează în acelaşi mod Totuşi unii sînt comuni pentru ambele cazuri Selectînd, reducînd şi comparînd listele de factori care influenţează funcţionarea altor sisteme de aşteptare, putem identifica un număr mic de factori comuni şi de semnificaţie majoră în toate cazurile; este vorba de factori primari Factorii care nu sînt comuni tuturor proceselor de aşteptare, apărînd doar în anumite procese de aşteptare sau chiar numai intr-unui, sînt factori secundari Vor trebui luate în considerare şi legăturile şi interdependenţele dintre factori Se poate conchide că factorii primari sînt : 1° frecvenţa şi modul de sosire al clienţilor; 2° repartiţia timpilor de servire; 3° ordinea în care clienţii sînt serviţi; 4° numărul de staţii de servire; 5° repartiţia timpilor de aşteptare ale clienţilor în general sîntem interesaţi în a afla efectele modificărilor factorilor 1 — 4 asupra factorului 5 Repartiţia timpului de aşteptare este o măsură a eficacităţii sistemelor de aşteptare (Vom vedea că mai există şi alte măsuri semnificative şi utile) în plus, putem trage concluzia că legăturile semnificative dintre factori sînt: Legătura dintre parametrii repartiţiei timpilor de servire şi frecvenţa medie de sosire a clienţilor Legătura dintre numărul de fire de aşteptare şi intensitatea medie de sosire Legătura dintre disciplina de aşteptare şi numărul de canale) în funcţiune Sîntem acum, într-o măsură mai mare sau mai mică, în situaţia de a începe formularea unui model general al sistemelor de aşteptare CARACTERISTICILE MODELELOR DE AŞTEPTARE Modelul matematic constă în cîteva legături simbolice între factorii primari Prin atribuirea unor forme particulare repartiţiiloi statistice asociate acestor factori primari şi prin atribuirea unor valori tuturor parametrilor incluşi în model, vom putea realiza previziuni Putem spune, de exemplu, care va fi media timpilor de aşteptare a clienţilor în caz că ar fi zece staţii de servire şi care ar fi această medie în caz că ar fi 12, toate celelalte elemente rămînînd aceleaşi Obiectivul unui model matematic simbolic este de a realiza previziuni posibile despre comportarea sistemelor sub influenţa modificării condiţiilor sau, uneori, despre sisteme care deocamdată nu există De aceea, teoria aşteptării este o teorie generală; este vorba nu de comportarea sistemelor de aşteptare propriu-zise, ci a unor modele simbolice abstracte ale acestora Este o idee importantă, ce trebuie reţinută Este uşor a cădea în greşeală crezînd că modelul ar fi însuşi sistemul Dezvoltarea teoriei aşteptării, aşa cum este indicată de natura principalelor modele studiate, pare a fi puternic influenţată de două lucruri, anume de originea ei, origine ce se află în studiul aglomeraţiei sistemelor telefonice, şi de aceea ce este uşor şi posibil în analiza matematică în cadrul acestui capitol vom observa că multe dintre presupunerile cel mai frecvent făcute, referitoare la studiile de aşteptare, sînt tocmai acelea care par realiste şi valide în realizarea şi studierea modelelor proceselor de aşteptare asociate sistemelor telefonice Acestea însă nu pot fi aplicate în toate ocaziile; ele sînt neadecvate mai cu seamă atunci cînd clienţii şi staţiile de servire sînt oameni Vom observa de asemenea că aceste presupuneri sînt chiar acelea care conduc la modele susceptibile de a le fi aplicată analiza matematică (Matematicienii nu sînt mai dornici decît alţii să-şi dedice preocupările unor lucruri inabordabile ) Vom analiza principalele caracteristici ale sistemelor de aşteptare şi în principal patru care sînt de cea mai mare utilitate în vederea unei clasificări Acestea sînt: fluxul de intrare, mecanismul de servire, disciplina de aşteptare şi numărul staţiilor de servire (canalelor) Terminologia a fost creată de D G Kendall într-o carte celebră prezentată la Societatea Regală de Statistică din Londra, în 1951 Această lucrare, poate mai mult decît altele, a avut ca urmare marele interes pentru teoria aşteptării, care a preocupat Anglia după anul 1950 Kendall şi-a extins nomenclatorul intr-un sistem de clasificare care este astăzi universal Clasificarea lui Kendall este de forma : A/B/s, unde A specifică modul de intrare, B este repartiţia timpilor de servire, iar s este numărul de staţii de servire In general se presupune că clienţii sînt serviţi în ordinea sosirii, deoarece codul clasificării nu permite definirea disciplinei de aşteptare Această schemă este însă prea vagă faţă de obiectivele acestei cărţi De aceea vom folosi o formă extinsă a notaţiei lui Kendall, anume A/B/s : (d/e), unde A, B, s au aceeaşi semnificaţie pe care le-a dat-o Kendall, d este numărul maxim de clienţi ce pot fi cuprinşi în sistem în orice moment (incluzînd pe cei ce sînt serviţi), iar e identifică disciplina de aşteptare Cîteva exemple vor lămuri aceste probleme M/M/l : (20/FIFO) reprezintă un model de aşteptare cu intrare aleatoare, timpi de servire aleatori, o staţie de servire; capacitatea maximă a sistemului este de 20 (adică în timp ce 19 persoane aşteaptă şi una este servită, toţi noii veniţi sînt refuzaţi şi pleacă), iar disciplina de aşteptare este primul venit, primul servit M/Ek/5- (oo/SIRO) înseamnă intrare aleatoare, timp de servire erlangian, cinci staţii de servire, capacitate a firului de aşteptare infinită (adică nici o intrare nu este refuzată), iar disciplina de aşteptare aleatoare Chiar şi această notaţie nu reuşeşte să clasifice complet modelele de aşteptare, dar este ceva mai utilă pentru scopul nostru, decît cea folosită de către Kendall FLUXUL DE INTRARE Apelul telefonic este un fenomen care nu rezultă din acţiunea comună a mai multor persoane, ci apare sub o formă individuală De aceea nu este surprinzător faptul că în majoritatea lucrărilor despre sistemele de aşteptare s-a presupus că sosirile sînt individuale Aceasta şi pentru a înlesni calculele matematice; ar fi însă greşit — orice motiv am avea s-o facem să considerăm sosirile ca fiind individuale, în unele cazuri în care este destul de clar că nu este deloc aşa în anumite sisteme de aşteptare, clienţii nu sosesc singuri, ci în grupuri De exemplu, pasagerii pentru cursele aeriene sosesc la aeroport în grupuri familiale, delegaţii de afaceri etc Primul pas ce trebuie făcut întotdeauna în examinarea unui sistem de aşteptare este de a stabili dacă sosirile sînt individuale sau în grup (agregate) Momentele de timp în care clienţii individuali sau grupurile intră în sistemul de aşteptare se numesc momente de sosire, iar intervalele de timp dintre două momente de sosire consecutive se numesc timpii dintre sosiri în studiile privind procesele de aşteptare este necesar a culege, datele privitoare la timpii dintre sosiri într-o formă convenabilă pentru a putea da răspuns la două întrebări: Timpii dintre sosiri sînt oare dependenţi, există oare între ei anumite corelaţii? Timpii dintre sosiri sînt oare corelaţi cu alte elemente ale sistemului de aşteptare? De cele mai multe ori se presupune că răspunsul la ambele întrebări este negativ Cititorul circumspect va crede în mod consecvent că răspunsurile pot fi presupuse negative atunci cînd procesele de aşteptare se referă la apelurile telefonice şi că răspunsurile negative conduc la o tratare matematică mai simplă Cu siguranţă că va avea dreptate Vom da un exemplu de timpi dintre sosiri care nu sînt independenţi Zborurile necesitate de schimburile comerciale sînt planificate de conducerea aeroportului Timpii de sosire sînt planificaţi cu mult timp înainte în perioadele de vîrf ale zilei (la marile aeroporturi) sosirile sînt planificate la intervale egale Dar avioanele sosesc mai devreme sau mai tîrziu Momentul real al sosirii oricărui avion se găseşte în apropierea momentului planificat Dacă întinderea acestei repartiţii a sosirilor este redusă în comparaţie cu intervalele de timp planificate între sosiri, nepunctualitatea poate fi neglijată în cele mai multe dintre cazuri Dacă însă este mai mare, sosirile pot fi considerate — din nou din raţiuni practice — aleatoare însă în cazul intermediar, cînd timpii planificaţi dintre sosiri sînt de ordinul primelor două dispersii ale abaterilor momentelor sosirilor reale de la momentele planificate, sosirile nu pot fi tratate ca regulate sau ca aleatoare Timpii consecutivi dintre sosiri sînt corelaţi între ei, corelaţie ce nu poate fi ignorată fără riscul unei erori Următorul pas ce trebuie făcut constă în stabilirea repartiţiei statistice a timpilor dintre sosiri şi pentru aceasta trebuie mai întîi să ne întrebăm : poate fi considerată oare această repartiţie ca fiind staţionară în timp ? Altfel spus, putem presupune oare că repartiţia timpilor dintre sosiri nu-şi modifică în timp proprietăţile statistice ? Oare probabilitatea unui interval cu durata de x minute, ales în mod aleator, este independentă de momentul zilei în care acest interval a fost măsurat? în studiile teoretice făcute se presupune de obicei că repartiţiile timpilor dintre sosiri sînt staţionare în realitate asemenea cazuri sînt rare Figura 2 1 reprezintă un model zilnic al convorbirilor telefonice dintr-un birou pentru rezervări de bilete de avioane Este limpede pe acest exemplu că media timpilor dintre sosiri este variabilă în timp (modul de variaţie este aproape sistematic cinci zile pe săp- tămînă) Fig 2 1 Variaţia frecvenţelor orare ale apelurilor telefonice intr-o zi tipică (luna mai 1959) Această discrepanţă dintre teorie şi realitate nu este în fapt foarte mare Majoritatea problemelor practice legate de procesele de aşteptare sînt problemele orelor de aglomeraţie Cu alte cuvinte, acestea devin foarte serioase cînd frecvenţa sosirilor este maximă După cum vom vedea, pe parcursul unor asemenea perioade, dacă acestea nu sînt excesiv de lungi, repartiţiile statistice ale timpilor dintre sosiri nu sînt de cele mai multe ori suficient de staţionare pentru a ne putea permite să le considerăm astfel Forma repartiţiei timpilor dintre sosiri, care a fost cel mai mult studiată, este cea corespunzătoare aşa-numitei repartiţii a sosirilor aleatoare, adică funcţiei exponenţiale negative Densitatea de probabilitate a unui timp dintre sosiri de mărime t este în acest caz a(t) = X e Z(, unde X este media sosirilor în unitatea de timp (unitate considerată în această carte a fi minutul, pentru cazul în care nu se specifică altceva) Aceasta este echivalent cu a spune că probabilitatea a n sosiri în unitatea de timp este dată de xn e~Ă (repartiţia Poisson) n! Semnificaţia acestei repartiţii este că probabilitatea XdZ a unei sosiri într-un interval mic de timp dt este constantă Ea nu este afectată de momentul în care a avut loc ultima sosire, nici de numărul sosirilor care au avut deja loc într-o exprimare mai tehnică, intervalele dintre sosirile consecutive sînt variabile aleatoare independente între ele, toate avînd aceeaşi funcţie de repartiţie Cititorul îşi va reaminti că transformata Laplace a unei funcţii exponenţiale negative este (1 - 2/X)"1 Funcţia generatoare corespunzătoare repartiţiei Poisson este Un alt tip de proces de intrare, care are o oarecare importanţă practică, este cel numit de către Cox şi Smith tipul cu intrare reînnoită Acesta poate avea loc atunci cînd intrările intr-un sistem de aşteptare sînt generate de ieşirile unui sistem anterior de aşteptare, sau de emisii de la un număr finit de surse, ce emit fiecare cîte un client la intervale ce urmează o anumită repartiţie De exemplu, să considerăm k asemenea surse A i-a sursă generează clienţi la anumite momente, astfel incit intervalul dintre două generări succesive este repartizat cu o densitate de probabilitate at{t) în acest caz fluxul combinat de sosiri format de ieşirile celor k surse are o repartiţie specială a timpilor dintre sosiri Este cunoscut faptul că în asemenea cazuri, forma statistică a sosirilor cumulate este aproximativ aleatoare pentru perioade mai scurte decît media intervalelor dintre sosirile de la toate sursele în practică, intrările în sistemele de aşteptare pot fi considerate ca fiind aleatoare MECANISMUL DE SERVIRE Intervalul dintre momentul în care clientul începe să fie servit şi momentul în care servirea este terminată este cunoscut sub denumirea de timp de servire Este important să facem distincţie între această mărime şi timpul de blocare al staţiei de servire Vom numi timp de blocare intervalul de timp dintre momentul în care începe timpul de servire al unui client şi momentul în care începe timpul de servire pentru următorul client, în cazul că acest al doilea client a stat la coadă pentru a fi servit Timpul de blocare poate fi conceput ca fiind timpul total în care un client împiedică servirea următorului, prin blocarea staţiei de servire Nu este necesar ca timpul de servire al unui client să coincidă cu timpul său de blocare De exemplu, să considerăm un cititor care împrumută o carte de la bibliotecă El soseşte la ghişeu, cere cartea dorită, prezintă legitimaţia, bibliotecarul îi caută cartea, o înregistrează şi i-o dă împreună cu legitimaţia Din punctul de vedere al cititorului, timpul său de servire este încheiat şi, ca atare, pleacă Dar bibliotecarul nu serveşte imediat cititorul următor El va pune înapoi fişa cărţii şi — dacă este cazul — o va înregsitra într-un dosar special (căci bibliotecarul Alderman Jones vrea să ştie proporţia cărţilor împrumutate din sectorul lucrărilor social-politice) Abia în acest moment el este gata pentru a servi următorul cititor : timpul de blocare s-a încheiat Diferenţa dintre timpul de servire şi timpul deblocare poate apărea, în acest caz, banală S-ar putea defini timpul de servire ca fiind chiar timpul de blocare, şi astfel s-ar putea utiliza formulele standard Este, fără îndoială, o consideraţie practică: munca realizată pe „timpul de servire” a clientului este, în general vorbind, obligatorie pentru ca clientul să elibereze staţia Cititorul nu va pleca fără cartea dorită Restul servirii, efectuat în timpul de blocare, nu este neapărat necesar să fie realizat în acest timp, în scopul de a satisface clientul Poate fi efectuat într-un alt moment al zilei în acest fel, numărul celor ce vor sta la rînd va scădea simţitor Metoda ar fi foarte eficientă dacă s-ar putea aplica tuturor sistemelor de servire, astfel încît, practic, timpii de servire să se confunde cu timpii de blocare Următorul exemplu va reliefa semnificaţia distincţiei ce o facem între aceşti timpi GHIŞEE 1 I 1 1 !xi i i i ■ 1 1 1 1 1 1 1 1 1 1 1 1 Ix1 lxl 1 1 1 1 lxl 1 1 1 1 'XI 1 1 i i i i 1 1 ■ t 1 1 1 i 1 1 1 1 1 1 1 1 X Client in deplasare spre ghişeu/ | liber —o Controlul şirului de osteptore X X X X SOSIRI Fig 2 2 Sistemul de aşteptare al controlului paşapoartelor Procesul de aşteptare ilustrat în fig 2 2 nu este ieşit din comun; el poate fi remarcat în operaţiile privind controlul paşapoartelor, la multe aeroporturi internaţionale Fie c ghişee de control al paşapoartelor, toate realizînd acelaşi tip de servire Cînd un turist ajunge la unul dintre aceste ghişee, ofiţerul de serviciu va începe executarea formalităţilor sau — dacă preferaţi — va începe ,,să-l servească” Timpii de servire pentru fiecare turist şi pentru fiecare ghişeu în parte sînt independenţi, identic repartizaţi Toţi turiştii sosiţi se aşază la coadă în ordinea strictă a sosirii, pe un singur şir Nimănui nu i se permite să înainteze spre ghişeu, pînă cînd turistul anterior nu a fost servit In sfera controlului paşapoartelor această regulă este de obicei controlată de un funcţionar plasat în faţa ghişeelor Dacă clientul ce urmează ar putea să se deplaseze cu viteza luminii la primul ghişeu liber, de îndată ce clientul anterior servit şi-a încheiat timpul său de servire, atunci timpul de servire şi timpul de blocare pentru acest client ar fi acelaşi Aceasta este situaţia în sistemele telefonice, unde „clienţii” sînt convorbirile care, evident, se pot deplasa cu viteză mare la o linie liberă Dar pasagerii ce trec pe la ghişeele de control al paşapoartelor nu pot să se deplaseze cu viteza luminii Ei merg cit mai calm şi mai degajat; graba poate produce suspiciuni ofiţerului de serviciu Mai mult, distanţele din capul şirului pînă la ghişee nu sînt identice Durează, în medie, mai mult a merge la ghişeele laterale decît la cel din centru Putem presupune că există c variabile x( (i = 1,2, , c), reprezentînd timpii de deplasare pînă la fiecare din cele c ghişee Timpul de blocare la ghişeul i este de aceea variabila v(, unde vf = v + x(, v fiind timpul de servire al clientului în majoritatea analizelor proceselor de aşteptare timpii de blocare sînt mai importanţi decît timpii de servire şi trebuie utilizaţi în locul timpilor de servire, orice formulă de aşteptare ar fi fost folosită Aşa stînd lucrurile, procesul de aşteptare de la controlul paşapoartelor nu este unul din cele mai simple; există mai multe staţii de servire şi fiecare are o altă repartiţie a timpului de servire Sperăm că cititorul este convins de semnificaţia deosebirii dintre timpul de servire şi timpul de blocare şi va considera că ei sînt trataţi, în majoritatea modelelor teoretice de aşteptare, ca identici Să mergem mai departe, atribuind deocamdată clienţilor calităţile obişnuite; vom presupune că timpul de servire şi timpul de blocare sînt identici Odată eliminată această dificultate, trecînd la examinarea proceselor de aşteptare, trebuie să ne punem patru întrebări, legate de timpii de servire : Oare toţi clienţii sînt de acelaşi tip, adică timpii lor de servire au aceeaşi funcţie de repartiţie? Oare toate staţiile de servire sînt identice, adică timpii de servire legaţi de staţiile de servire au aceeaşi funcţie de repartiţie? Timpii de servire sînt oare corelaţi cu alte aspecte ale sistemului? Eeparţiţia timpului de servire este oare staţionară? La începerea cercetării se speră întotdeauna că răspunsurile la aceste întrebări vor fi afirmative în lucrările teoretice, de cele mai multe ori aşa se şi presupune ; din păcate, însă, în practică aceste cazuri nu sînt prea întîlnite Cîteva exemplificări o vor dovedi Cargouri ce transportă minereuri intră într-un port pentru a-şi descărca încărcătura de minereu de fier Ele au înzestrări şi capacităţi diferite De aceea, timpii necesari pentru a le descărca — timpii de servire — nu sînt omogeni din punct de vedere statistic Ei au deci repartiţii diferite în cadrul procesului de control al paşapoartelor timpul de blocare, care, în orice analiză a procesului, trebuie să înlocuiască timpul de servire, depinde de repartiţii diferite ce corespund diferitelor ghişee Lingourile de oţel sosesc la laminoare de la cuptoare deschise Ele sînt încă incandescente, dar prezintă variaţii de temperatură Ele trebuie plasate în cuptoare de reîncălzire (puţ infiltrat în pămînt) înainte de a fi laminate, operaţia de laminare necesitînd o temperatură înaltă a lingourilor Dacă cuptoarele de reîncălzire sînt pline în momentul în care soseşte un lingou, acesta trebuie să aştepte în timpul acestei aşteptări, el se răceşte Cu cît aşteaptă mai mult, cu atît va trebui menţinut mai mult în cuptorul de reîncălzire Timpul de servire al acestei operaţii este timpul în care lingoul ocupă cuptorul de reîncălzire De aceea, în acest caz, timpul de servire este direct proporţional cu timpul de aşteptare, iar procesul este autoagravant Agenţii liniilor aeriene civile lucrează opt ore pe schimb, la un aeroport, verificînd pasagerii pentru diferitele zboruri Anumite studii legate de operaţiile de verificare, făcute atît în Anglia cît şi în S U A , au arătat că ele obosesc agenţii şi că timpii deservire sînt mai lungi şi mai variabili spre sfîrşitul schimbului Apare astfel evident că răspunsurile afirmative de la primele pentru întrebări nu pot fi întîlnite întotdeauna în lumea reală Deoarece teoria, în general, presupune contrariul, ne vom ocupa în continuare de anumite procese de aşteptare pentru a reliefa cît de eronat poate fi un raţionament care presupune afirmativ răspunsurile de mai sus şi foloseşte modele teoretice gata făcute, acolo unde este clar că nu sînt adecvate Desigur, cititorul va începe să creadă că, neputînd fi folosite în acest fel, multe modele de aşteptare clasice familiare nu pot fi folosite la nimic altceva, decît în practica telecomunicaţiilor Deci întrebarea „Cît de temeinice sînt aceste modele?” nu este deloc de ignorat Acum, cînd unele aspecte au fost clarificate, vom trece sumar în revistă tipurile de repartiţie a timpului de servire care sînt cel mai des utilizate şi care au fost cel mai mult studiate Vom începe cu repartiţia exponenţială negativă Repartiţia exponenţială negativă a mai fost considerată atunci cînd s-a vorbit despre intrările în sistemul de aşteptare în acest caz, densitatea de repartiţie a timpului de servire de durată t este dată de r elaţia 6(f) = ue-^, unde g este numărul mediu de clienţi serviţi în unitatea de timp sau, altfel spus, b = 1/p este timpul mediu de servire a unui client Din punct de vedere matematic, aceasta este cea mai convenabilă formă a repartiţiei timpului de servire Explicaţia este că timpul rămas pînă la efectuarea deplină a unui serviciu este independent de timpul pentru care serviciul a fost în progres Aceasta înseamnă, de exemplu, că dacă un client găseşte doar o persoană înaintea sa şi această persoană este în curs de servire, atunci timpul de aşteptare este repartizat exact ca şi timpul de servire Acest lucru simplifică foarte mult analiza Repartiţia exponenţială negativă a timpului de servire nu este prea comună, dar nu sînt rare cazurile în care repartiţiile timpilor de servire observaţi nu diferă prea mult de repartiţia exponenţială negativă Cînd apare o deviaţie observabilă, este adesea posibil să o aproximăm cu repartiţia Pearson de tipul III Există suficiente motive matematice pentru alegerea unei repartiţii de acest gen De fapt, în general, este posibil a aproxima astfel de date cu repartiţia logaritmică normală (vezi, de exemplu, cap 15) Dar repartiţia loga- ritmică normală nu este tot atît de convenabil de manevrat, nici nu admite o interpretare mai convenabilă decît repartiţia Pearson de tipul III Repartiţia Pearson de tipul III, care este utilizată în modelele proceselor de aşteptare, are densitatea de repartiţie Atunci cînd k este un întreg, formula poate fi exprimată, în forma sa erlangiană (numită astfel după danezul Karl Erlang, cel care a utilizat această formă în studiile firelor de aşteptare): (^)fc (fc-1)’ în această formă, repartiţia este deci aceea a unei sume de k variabile independente, identic repartizate exponenţial cu media l/(fifc) Aceasta este cea mai utilă legătură între repartiţia Erlang şi cea exponenţială Repartiţia poate fi privită ca o modificare a repartiţiei X2 cu media 1/ji şi cu 2k grade de libertate Cititorul va remarca existenţa a doi parametri în formula repartiţiei, ţi şi k ; astfel, ambii pot fi estimaţi utilizînd un set de observaţii Repartiţia este foarte asimetrică pentru valori mici ale lui k ; pentru k = 1 este tocmai repartiţia exponenţială negativă Pe măsură ce k creşte, dispersia începe să descrească, pînă cînd, pentru un k suficient de mare repartiţia devine aproape constantă Media repartiţiei este 1/ji, iar coeficientul de variaţie este l/]/k Transformata Laplace a repartiţiei este [1 + Al treilea tip important de repartiţie a timpului de servire este repartiţia constantă Aceasta însă nu necesită nici o prezentare 2 2 3 DISCIPLINA DE AŞTEPTARE în teoria aşteptării se presupune ca ipoteză majoră asupra disciplinei de aşteptare faptul că clienţii sînt serviţi în ordinea sosirii Această regulă este numită FIFO, ceea ce înseamnă „primul venit este primul servit" Disciplina de aşteptare este elementul esenţial de corectitudine din punctul de vedere al clienţilor, fiecare fiind preocupat de propriul său timp de aşteptare Ea nu este neapărat esenţială din punctul de vedere al staţiei de servire şi poate fi chiar nesemnificativă în situaţiile în care clienţii sînt indiferenţi faţă de timpul cît aşteaptă Merită să examinăm implicaţiile regulii FIFO în acest scop să considerăm sistemul de aşteptare de la controlul paşapoartelor care a fost descris în secţiunea precedentă Regula FIFO poate fi realizată doar într-un sistem de aşteptare cu mai multe staţii de servire, cînd : Clienţii sînt îndreptaţi, în ordinea sosirii lor, spre un singur şir de aşteptare şi rugaţi să se deplaseze mai departe spre capul şirului Un client poate părăsi şirul numai în clipa cînd a ajuns în faţă şi i-a fost permis printr-un semnal să facă aceasta Ne permitem să comentăm aplicabilitatea acestei proceduri Mai întîi, aşa după cum remarcam şi în capitolul precedent, regula FIFO poate fi foarte bine aplicată la sistemele telefonice, acolo unde clienţii — în acest caz, convorbirile — sînt obligaţi, datorită circuitelor, „să înainteze pe o singură linie” Regula poate fi utilizată, de asemenea, în condiţiile clienţilor- oameni, prin adoptarea unei simple ordonări, întîlnită de obicei în magazinele alimentare nord-americane Fiecare client, la sosire, primeşte un tichet numerotat Atunci cînd un client a fost servit, vînzătorul va striga numărul următor Persoana care are tichetul cu numărul menţionat va înainta pentru a fi servită şi va înmîna tichetul vînzătorului, Regulile FIFO privind disciplina de aşteptare sînt de fapt foarte larg răspîndite Ele tind să apară ori de cîte ori: Clienţii sînt oameni sau fenomene abstracte asociate oamenilor — cum sînt convorbirile telefonice — care sînt foarte afectaţi de durata de aşteptare sau de tratamentul diferenţiat Există restricţii de natură fizică (circuitele telefonice, numărul tichetelor) sau psihică (morală) (de exemplu atitudinea pasagerilor într-o staţie de autobuz din Londra) care obligă să se înainteze în mod ordonat Ar fi însă greşit să considerăm că regulile FIFO ar fi universale, chiar atunci cînd clienţii sînt oameni De exemplu, să considerăm oamenii care aşteaptă autobuzul într-o staţie din Roma Regula de aici ar putea fi comparată cu o altă regulă care a fost studiată destul de amănunţit: servirea în ordine aleatoare, cu alte cuvinte, regula SIRO Regula SIRO a disciplinei de aşteptare poate fi întîlnită în multe situaţii ca, de exemplu, în telefonie, în operaţiile manuale efectuate la un tablou de comutaţie, la vînzarea răcoritoarelor la bufetul unui teatru, în pauzele piesei; la rezervarea de cuşete pentru pasageri (Nu este o glumă referitoare la căile ferate; de fapt, sînt foarte frecvente cazurile în care companiile feroviare nu pot comunica cu pasagerii înregistraţi pentru rezervări, sunîndu-i la numărul de telefon pe care aceştia l-au lăsat, şi vor continua căutările pînă cînd, printr-o şansă, vor găsi pe cineva care mai doreşte o cuşetă ) Există anumite situaţii speciale, aşa după cum vom vedea, în care regula LIFO (ultimul venit, primul servit) a disciplinei cozii ar fi foarte potrivită De exemplu, în cazul lingourilor de oţel ce aşteaptă să intre la reîncălzire, regula LIFO a disciplinei de aşteptare ar putea fi cea mai eficientă Prin utilizarea ei, căldura ar putea fi conservată, s-ar scurta timpul de servire şi în cele din urmă s-ar reduce şi timpul de aşteptare Se pare că aşa este, ţinînd seama de legea lui Ştefan a răcirii Alte tipuri de disciplină a aşteptării sînt asociate unor reguli de prioritate, atunci cînd există mai multe tipuri de clienţi Să presupunem că clienţii fac parte din k clase şi fiecare clasă prezintă o anumită prioritate Clienţii din clasa a i-a sosesc urmînd o repartiţie dată de funcţia A((t) şi sînt serviţi, repartiţia timpului de servire fiind dată de funcţia Bf(t) Clienţii componenţi ai unei aceleiaşi clase sînt serviţi după regula FIFO, dar clienţii cu prioritate i sînt serviţi înaintea celor cu prioritatea i+1 Există două feluri de priorităţi Dacă este în curs de servire un client cu prioritatea i +1, acesta va putea fi înlocuit şi reîntors în şirul de aşteptare în momentul cînd un client cu prioritate mai mare intră în sistem Aceasta se numeşte prioritate absolută Dacă totuşi clientului în curs de servire i se poate termina serviciul, regula de prioritate urmînd să fie operativă după aceasta, avem de-a face cu un serviciu cu prioritate relativă în cadrul unui model, vom nota serviciul cu priorităţi absolute cu codul PSPO (ordinea de servire a priorităţilor absolute), iar priorităţile relative NPPS (ordinea de servire a priorităţilor relativă) Este destul de dificil să tratezi aceste priorităţi din punct de vedere matematic, exceptînd cazurile modelelor simple, şi pînă în prezent nu li s-a acordat prea multă atenţie Acest aspect al priorităţilor ridică probleme interesante legate de cost, care nu a fost incorporat în modelele de aşteptare aşa cum ar fi trebuit De exemplu, costurile pot fi introduse pentru a compara pierderile cauzate la neacordarea unei priorităţi mai mici clientului aflat în serviciu, cu pierderile totale sau parţiale rezultînd din eventuala lui eliminare în concluzie, se pot observa o mulţime de discipline de aşteptare ce operează în multe sisteme reale de aşteptare, dar este dificil a le formula matematic De exemplu, există fenomenul deplasării; deplasări individuale ale clienţilor de la un şir la altul, în scopul găsirii unui avantaj personal (aceasta în cazul cînd există un şir pentru fiecare din staţiile de servire, toate staţiile fiind identice) Pînă nu de mult, nimeni nu a dat prea multă atenţie acestui fenomen, fie şi de o manieră empirică Există, aşadar, multe tipuri de discipline de aşteptare, şi pentru că regulile care guvernează ordinea servirii sînt printre aspectele foarte uşor controlabile şi coordonabile ale proceselor de aşteptare, ele prezintă un mare interes pentru cercetătorul operaţional 2 2 4 NUMĂRUL STAŢIILOR DE SERVIRE Cînd într-un sistem de aşteptare sînt mai multe staţii disponibile pentru a realiza servirea, elementul esenţial devine organizarea acesteia Există două tipuri principale de organizare a servirii: în paralel şi în serie Prin staţii de servire în paralel vom înţelege un număr de staţii avînd facilităţi de servire identice, astfel incit un client să poată fi servit la fel de bine la oricare dintre ele Cînd un client poate selecta sau poate fi îndrumat la oricare dintre staţii, vom spune că staţiile sînt complet abordabile Există, după cum vom vedea, posibilitatea ca să fie mai multe staţii în funcţiune, dar care nu sînt complet abordabile, adică ele sînt complet abordabile doar de către anumiţi clienţi, dar nu şi de către ceilalţi Cînd are loc acest lucru, vom vorbi despre abordarea restrictivă Staţiile de servire în serie se întîlnesc de obicei atunci cind sistemul de servire este defalcat în mai multe elemente Atunci, pentru a realiza servirea, se formează o serie de şiruri de aşteptare, fiecare staţie fiind destinată unui element al servirii Clientul trece pe rînd pe la toate staţiile, începînd cu prima El va avea de aşteptat pentru a fi servit la toate fazele procesului Timpii de aşteptare ai clientului şi timpii totali de aşteptare sînt puternic afectaţi de forma organizării şi este important de a decide corect modul optim de alocare a capacităţii de servire în general, în această carte, ne vom ocupa foarte sumar de sistemele de aşteptare cu staţii de servire în serie Aceasta nu pentru că problema ar fi lipsită de importanţă, ci datorită dificultăţii de a găsi cazuri potrivite pentru a le discuta MODELE CU O SINGURĂ STAŢIE DE SERVIRE 3 1 GENERALITĂŢI Se obişnuieşte să se înceapă expunerile în teoria aşteptării cu o prezentare a diferitelor modele cu o singură staţie de servire şi în particular să se prezinte în detaliu aparatul matematic al celui mai cunoscut şi mai utilizat dintre ele, adică al modelului Jf/Jf/1: (oo/FIFO) Este limpede că motivele pentru aceasta sînt în primul rînd de natură pedagogică însemnătatea practică a modelelor cu o staţie de servire, de orice gen ar fi, este relativ neînsemnată Atunci cînd se descoperă unele neajunsuri în funcţionarea unui sistem cu o staţie de servire, de obicei se încearcă să se evalueze posibilităţile transformării lui într-un sistem cu două staţii de servire Meritele pedagogice ale modelelor cu o staţie de servire nu trebuie însă ignorate Cu toată aplicabilitatea lor practică directă limitată, vom păstra tradiţia şi vom începe analiza cu ele înainte de a trece la discutarea diferitelor modele cu o staţie de servire, trebuie să menţionăm că mulţi autori au studiat modele generale de acest tip şi ,,le-au rezolvat”, din punct de vedere matematic Cu alte cuvinte, ei a considerat modele de tipul GI/G/1 : (oo/FIFO) şi au exprimat repartiţiile timpilor de aşteptare şi a timpilor totali de aşteptare De obicei, formulele obţinute corelează transformatele Laplace ale repartiţiilor timpilor de aşteptare şi a timpilor totali de aşteptare, cu tansformatele Laplace ale repartiţiilor timpilor dintre sosiri şi timpilor de servire Ele nu sînt întotdeauna soluţii în sensul folosit de noi deoarece, cu excepţia formelor celor mai simple ale timpilor de aşteptare şi a timpilor de servire, transformatele Laplace nu pot fi inversate în mod efectiv Este adevărat că sînt posibile evaluări numerice, dar de cele mai multe ori este foarte dificil a le obţine Trebuie să menţionăm totuşi că majoritatea rezultatelor privitoare la modelele cu o staţie de servire ce vor fi discutate în acest capitol pot fi obţinute fără prea mari dificultăţi dintr-un model foarte general Motivul pentru care nu vom proceda astfel constă în dorinţa de a prezenta şi alte tehnici ce sînt foarte utilizate 3 2 MODELE CU TIMPI DE SERVIRE EXPONENŢIALI 3 2 1 MIM/l: (oo/FIFO) Acesta este cel mai simplu model teoretic dintre cele ce au pretenţii de aplicabilitate Vom reaminti caracteristicile lui, precum şi denumirile cîtorva parametri; astfel: M Sosirile sînt aleatoare, dar au loc cu intensitate medie constantă (în această carte, unitatea de timp este minutul, în cazul în care nu este specificat altceva) M Repartiţia timpului de servire este exponenţială negativă cu o valoare medie constantă, 1/p minute 1 Există doar o staţie de servire oo Nu există vreo limită legată de lungimea şirului, clienţii nu sînt respinşi de sistem, iar odată aşezaţi în şirul de aşteptare nu-1 mai pot părăsi FIFO Clienţii sînt serviţi în ordinea strictă a sosirilor Modelul este uşor de mînuit, datorită unei ciudate proprietăţi a repartiţiei exponenţiale negative, ceea ce Saaty numeşte „uitarea” Prin aceasta ne referim la faptul că dacă timpii de servire au această repartiţie, atunci probabilitatea terminării serviciului la un moment dat este independentă de durata serviciului Consecinţa acestei proprietăţi este că, stabilind ecuaţiile pentru a releva modificările survenite în cadrul modelului, nu este necesar să ţinem seamă de cit timp a trecut de la momentul ultimei sosiri, nici cit de mult a trecut pînă cînd a început timpul curent de servire Vom introduce acum simboluri (o listă completă a simbolurilor utilizate se găseşte în Anexa 2) Probabilitatea de a avea o sosire pe parcursul unui interval de timp foarte scurt (Z, t + dZ) este XdZ în mod similar, dacă notăm ţi = 1/b, atunci probabilitatea ca unei persoane să i se termine servirea pe parcursul intervalului de timp ales (presupunînd că au existat încă de la început persoane în curs de servire) este p dZ în fine, fie pn(t) probabilitatea ca la momentul Z să fie n clienţi în sistem Probabilitatea pn{t + dZ) ca sistemul să conţină n persoane la momentul Z + dZ are trei componente : — Probabilitatea pn(t) ca sistemul să conţină exact n persoane la momentul anterior Z, înmulţită cu probabilitatea de a nu se înregistra nici o plecare sau sosire în intervalul de timp dZ dintre aceste momente Cum dZ este foarte mic, această ultimă probabilitate este 1 — XdZ — u dZ iar contribuţia totală a acestei componente este • p„(Z) — Probabilitatea Pn-^t) ca sistemul să conţină exact n—1 persoane la momentul anterior Z, înmulţită cu probabilitatea jjl dZ a unei sosiri în intervalul dZ (Cum dZ este foarte mic, probabilitatea ca să existe mai mult de o sosire este neglijabilă ) Contribuţia acestui termen este XdZ • p„ ! (Z) — Probabilitatea ca sistemul să conţină n + 1 persoane la momentul precedent Z, înmulţită cu probabilitatea jxdt a unei plecări în intervalul dZ (Probabilitatea mai multor plecări este neglijabilă ) Contribuţia acestui termen este în consecinţă, forma generală a relaţiei este pn(t + dt) = pn(t) + Xdt p„ 1(t) + |idtpB+1 (t), şi după o rearanjare a termenilor, împărţind prin dt şi trecînd la limită pentru dt tinzînd către zero, rezultă că J - Pn(t) = - (X + + Xp^t) + fipfl + 1(t) (3 1) dt care este adevărată pentru n > 1 Pentru n = 0, ceea ce înseamnă că staţia este liberă (nu există nici o persoană care să aştepte şi nici una care să fie servită), avem -7-P0W = - x2»o(O + ^i(0- (3 1a) dt Aceste ecuaţii diferenţiale constituie modelul matematic, care poate fi rezolvat, dacă se cunoaşte numărul de persoane aflate la coadă la momentul iniţial Este mai convenabil să se ia momentul iniţial t = 0, deşi acest lucru nu este absolut obligatoriu De exemplu, dacă ar fi m persoane la coadă la momentul 0, condiţiile iniţiale ar fi pm(0) = 1, p„(0) = 0, n m Cu această informaţie, ecuaţiile scrise mai sus permit calcularea probabilităţilor pn(t) pentru orice t şi pentru toate valorile n = 0,1, 2, Un caz important apare atunci cînd m = 0 Aceasta înseamnă că la momentul iniţial nu există clienţi în aşteptare şi deci staţia de servire este liberă Condiţiile iniţiale sînt 2>o(°)= 1, l>n(0)=0, n > 0 Cea mai avantajoasă metodă de obţinere a unor relaţii explicite pentru p„(t) aparţine lui Cox şi Smith* Această metodă învinge dificultăţile legate de faptul că ecuaţia (3 1 a) are o formă diferită faţă de cazul general, reprezentat de ecuaţiile (3 1) înmulţind ecuaţiile (3 1) cu s”, unde s este o variabilă arbitrară, continuă şi nenegativă, însumînd după valorile întregi ale lui n, — oo 2 n — - oo Presupunem că K(s) — £ I rsr Să înlocuim această expresie în re- T~ - QO laţia de mai sus şi să egalăm coeficienţii puterilor lui s Din condiţiile iniţiale po(O) = 1 şi p„(0 = 0 pentru n > 0 şi proprietăţile funcţiilor Bessel [Zo(O) = 1 şi Z«(0) = 0, n >0], rezultă = 1, Â r = 0 pentru r = 1, 2, în consecinţă, pa(t) = exp {- (X + jx) t} [ p“''2 I„(x) + £ kT ?ln+r^In r (x) • Valorile kr sînt determinate astfel incit să asigure satisfacerea condiţiei p0(t) = pp-x(t) înlocuirea directă, urmată de egalarea coeficienţilor funcţiilor Ir(x), arată că fci= p *, k2 — (1 p) p *, k3 = (1 p)p ■ etc Ca urmare, in final vom obţine probabilităţile dependente de timp, de forma pn(t) = exp{— (X + [x)/}x X [>2 I„{x) + p,n-I,/2In+1(^) + P"(l - P) Ş p- ,w+r,/2lB+r(^)j (3 4) Relaţia este adevărată pentru toţi p nenegativi (chiar şi pentru p>l) Pus în faţa acestei relaţii, cititorul, avînd în minte ideea că acest sistem este cel mai simplu dintre toate, va avea poate un sentiment de panică Acest sentiment va creşte cu siguranţă dacă va ieşi la iveală faptul că soluţia dată descrie doar creşterea şirului de aşteptare (exprimată în număr de clienţi ce aşteaptă), nespunîndu-ne nimic despre timpii de aşteptare a clienţilor Pentru a-i obţine şi pe aceştia va trebui să considerăm expresii de formă mult mai complicată Dacă acestea sînt dificultăţile găsirii soluţiei pentru modelele cele mai simple, ce se va întîmpla cu modelele mai complicate ? Răspunsul a această întrebare este dat de prima regulă de lucru din teoria aşteptării: soluţiile dependente de timp pentru modelele de aşteptare sînt sau imposibil de obţinut sau imposibil de manipulat Este însă liniştitor să constatăm că ocaziile în care se apelează la cercetători operaţionali pentru a analiza sisteme de aşteptare de acest fel sînt foarte rare De asemenea, sînt rare cazurile în care se cere să se studieze un proces tocmai în decursul modificării prametrilor de sosire sau de servire De obicei, nu se cere să se facă o analiză a sistemelor de aşteptare pînă nu s-au descoperit dificultăţile procesului, şi de aceea, prin definiţie, a trecut deja o perioadă de timp Lăsînd deoparte acest mod simplist de a privi lucrurile, trebuie să constatăm că de fapt cea mai mare parte a problemelor practice legate de sistemele de aşteptare nu prea au tangenţă cu comportarea tranziţională în cazul în care staţiile de servire şi/sau clienţii sînt oameni, perioadele de început sînt marcate de o asemenea instabilitate a repartiţiilor timpilor de servire, încît este practic imposibil să se apeleze la analiza matematică în mod frecvent, repartiţiile timpilor dintre sosiri sînt şi ele instabile în aceste condiţii Ar fi deci absurd să cerem să aplicăm o soluţie exactă a modelului în astfel de circumstanţe De aceea, în cele ce urmează, exceptînd un singur caz, nu vom mai acorda atenţie soluţiilor dependente de timp şi formelor riguroase ale comportării tranziţionale Iar atunci cînd va fi necesar să studiem astfel de probleme, vom utiliza metoda de aproximare sau tehnici de simulare Problemele majore ce se pun aici sînt de forma : oare cîte staţii de servire sînt necesare, presupunînd că ele lucrează cu o intensitate normală? Intensităţi normale de obicei nu apar deloc în perioadele de început Aceasta, ca şi dificultăţile matematice, explică de ce cea mai mare parte a studiilor făcute asupra sistemelor de aşteptare au fost dedicate sistemelor ce se află în echilibru statistic Un sistem de aşteptare se află in echilibru statistic dacă probabilitatea existenţei a n persoane în sistem (pentru orice n > 0) este independentă de timpul ce a trecut de la începerea procesului Cu alte cuvinte : Pn(t) = Pn- Aceasta nu este, poate, cea mai elegantă definiţie, dar ne va fi suficientă Cititorul va observa că o altă formă de a exprima cele de mai sus este de a spune că în echilibrul statistic (sau de starea stabilă) avem —pn(t) = 0, n > 0 di Astfel, în cazul nostru simplu, modelul se reduce la (3 5) (3 6) (3 6) (3 7) (3 7)(X + p )pn = XpB ! + p pn+1; n > 1 Aceste ecuaţii pot fi uşor rezolvate, dînd : Pn = (1 — p) p”, n = 0,1, , unde p = X/[i este intensitatea traficului măsurată în unităţi internaţionale de măsurare a traficului numite erlangi* Repartiţia lungimii şirului de aşteptare este geometrică Este limpede că aceasta este o funcţie de repartiţie numai în cazul în care p este mai mic decît 1 Cu alte cuvinte, echilibrul statistic există numai în condiţiile în care intensitatea traficului este mai mică decît 1, adică atunci cînd staţia de servire poate, în medie, servi clienţii mai repede decît sosesc ei Aceste aspecte evidente nu necesită o justificare matematică De asemenea, este clar că dacă clienţii sosesc mai repede decît poate să-i servească staţia, şirul se va lungi continuu, iar probabilitatea ca şirul de aşteptare să aibă o lungime dată va depinde de timpul scurs Nu este atît de evident că în cazul în care intensitatea sosirii şi intensitatea de servire sînt egale, adică atunci cînd p = 1, echilibrul statistic este din nou imposibil Dar, în cazul în care p= 1, în expresia de mai sus pn sînt toţi egali cu zero şi, fără îndoială, nu dau o funcţie de repartiţie Următoarea noastră sarcină este de a exprima repartiţiile timpului total de aşteptare şi a timpului de aşteptare în şir Transformata Laplace a repartiţiei timpului de servire este [x ( e lt‘ e~zl di = Jo [X + 2 în consecinţă, transformata Laplace a sumei de r astfel de repartiţii Ca urmare, transformata Laplace a timpului total de aşteptare este - p) (*(i — p) + z’ - p) (*(i — p) + z’ însă aceasta este tocmai transformata Laplace a expresiei (3 8) (3 8)((i - X)e- ->', care este densitatea de repartiţie a timpului total de aşteptare Timpul mediu de aşteptare este obţinut prin diferenţierea transformatei Laplace şi făcînd apoi z = 0 din cauză că transformata Laplace este funcţia generatoare a momentelor; astfel: H(1 — P) în mod similar se obţine repartiţia timpilor de aşteptare în şir Ea are valoarea 1 — p în 0 şi, pentru t > 0, este dată de densitatea de repartiţie p(X — p) e (X 1I)( REPARTIŢIA PERIOADEI OCUPATE Perioada ocupată este, în terminologia matematică a aşteptării, intervalul de timp care începe atunci cînd soseşte un client care nu găseşte pe nimeni în faţa lui şi se sfîrşeşte cu prima ocazie cînd, la plecarea unui client, acesta nu lasă pe nimeni în urma sa Nu este deloc acelaşi lucru cu „ora de aglomeraţie" care semnifică perioada în care intensitatea sosirilor este foarte ridicată Această deosebire trebuie reţinută, deoarece expresia perioadă ocupată va fi utilizată din cînd în cînd şi în acest ultim sens Care este repartiţia perioadelor ocupate? Aproape că am şi determinat-o Este dată de (d/dt)p0(t), căci aceasta măsoară frecvenţa de apariţie a unei situaţii fără aşteptare Ştim că -Ţ-^o( 1, (3 10) dt n +1 -^-îo(0 =XţZi(t) - (p +X) g0(t) (3 11) dt Condiţiile iniţiale sînt g„(0) = 1 pentru toate valorile lui n, deoarece oricît de devreme s-ar aşeza clientul nostru în şir, este sigur că el va sta acolo cel puţin un interval de timp zero Vom presupune de asemenea că intensitatea medie de servire este superioară intensităţii medii a sosirilor, adică p 0, 8(0) = 1 Apare astfel evident că [ă(n)]2 = — 8(n) şi n[l — ă(n)] = n în echilibru statistic, dacă numărul mediu de persoane care aşteaptă este n, atunci n = n' în mod similar D(n) = D(n') Prin urmare, pe parcursul timpului de servire de durată t 8(n) = 1 -r = 1 - kt Şi D(r) = X2/2 + kt Deoarece timpul mediu de servire este 1/p , avem r = X/p = p Probabilitatea a r sosiri este independentă de n şi de 8(n) Aplicînd regulile statistice obişnuite privind valorile medii ale produsului de variabile independente, avem M[r2] = X2 D(i) + p2 + P De asemenea, n’2 = n2 — 2n(l — r) + (r — l)2 + S(2r — 1) care, în stare de echilibru, dă M[n'2] - M[n2] = 2M[n] M[r - 1] + M[(r - l)2] + + M • M , relaţie care poate fi combinată cu relaţiile anterioare şi din care se obţine (3 16) (3 16)■ p2 — X2 D(t) 2(1 - P) pentru numărul mediu de persoane existente în sistem în momentele de după plecări Acesta nu poate diferi cu mai mult de o persoană de numărul mediu de persoane prezentate la un moment dat Este uşor de arătat că media timpului de aşteptare este (3 17) (3 17)w = p2+ X2D(t) P 2X(1 - p) Argumentul care conduce la acest rezultat este acela că numărul mediu de sosiri pe parcursul timpului mediu de aşteptare este X(W9 + 1/p ), pe care îl egalăm cu n Aşadar : de unde rezultă (3 17) x(jya + î/ji) = P x(jya + î/ji) = P p2 + X2 D(t) , 2(1 - p) p2 + X2 D(t) , 2(1 - p) Aceste rezultate sînt cunoscute ca formulele Pollaczek-Hincin Ele se utilizează pentru examinarea efectelor diferitelor repartiţii ale timpilor de servire M/Ekll : (oo/FIFO) Doar modelele cu sosiri aleatoare şi repartiţii exponenţiale negative ale timpului de servire conduc la ecuaţii liniare simple Pentru a' te tipuri de repartiţii, sistemul de ecuaţii este de obicei mult mai complicat şi în general el nu poate fi rezolvat explicit cu ajutorul unor funcţii simple Majoritatea proceselor de aşteptare care au loc nu au repartiţii exponenţiale ale timpului de servire Ce este atunci de făcut, pentru a crea modele uşor de analizat pentru asemenea procese? Karl Erlang, în lucrarea sa de pionierat, sugerează un răspuns la această problemă El observă că cea mai mare parte a repartiţiilor neexponenţiale ale timpului de servire pe care le-a întîlnit pot fi aproximate prin funcţia de repartiţie cunoscută sub numele de repartiţia Euler, repartiţia Pearson de tipul III, repartiţia x2, repartiţia Erlang şi repartiţia gamma (Ar fi interesant de aflat explicaţia creşterii numărului de astfel de nume ) Vom utiliza densitatea de repartiţie, egală cu dL (fc-1)! Aşa după cum am văzut (cap 2) aceasta este repartiţia sumei a fc variabile aleatoare, fiecare repartizată exponenţial cu media l/( p /c) Erlang şi-a dat seama că un mecanism de servire cu o astfel de repartiţie poate fi conceput ca fiind format din k etape în serie, fiecare etapă fiind o staţie de servire exponenţială El nu a afirmat că mecanismul de servire necesită o structură în k faze; tot ceea ce este necesar pentru scopuri analitice este ca acesta ar putea să se comporte ca şi cum ar fi aşa Cititorul va vedea că dacă un model poate fi construit pe baza elementelor exponenţiale, ecuaţiile ce trebuie să fie rezolvate vor fi mult simplificate Regulile care guvernează activitatea unei faze-tip a mecanismului de servire cu k faze (sau etape) sînt: Numai un client poate fi servit la un anumit moment dat de staţia de servire Fiecare client trebuie să străbată fiecare fază a serviciului, începînd cu prima După terminarea unei faze, clientul se deplasează imediat la faza următoare Atunci cind un client părăseşte faza k, următorul client poate intra în faza 1, dai’ nu înainte Cititorul va remarca că aceasta nu este acelaşi lucru cu k servicii în serie, deoarece în cazul de care ne ocupăm, numai un client poate fi servit la un anumit moment dat, în timp ce în celălalt caz puteau fi implicaţi în mecanismul de servire pînăla k persoane Datorită acestei condiţii date de regula 1, putem considera cele k faze echivalente cu o singură operaţie Modelul de aşteptare M/£t/l : (oo/FIFO) poate fi studiat utilizînd metoda Erlang Un mod de a stabili ecuaţiile constă în a nota probabilitatea ca la momentul t să existe n clienţi prezenţi iar clientul care este servit să se afle în faza i, cu simbolul pni(t) Aceasta este totuşi o complicaţie inutilă Este preferabil să notăm cu pn(t) probabilitatea sil existe nfaze de aşteptare şi de servire la momentul t; astfel, fiecare client nou sosit va face să crească numărul de faze existent cu k, şi ori de cîte ori o fază este completată, numărul lor este redus cu o unitate Ecuaţiile stării stabile pentru această situaţie pot fi construite in maniera obişnuită, examinînd toate posibilităţile Aceste ecuaţii sînt : (3 18) (X + = *fipn+1 + Xp„ t, n > 1 (3 19) Funcţia generatoare poate fi dedusă din aceste ecuaţii în mod obişnuit, prin înmulţirea membrului sting al ecuaţiei în pn cu s" şi însumînd după valorile lui n, de la 0 la oo Astfel: F(s) = - 1 F(s) = - 1 (3 20) (3 20)Pa = A>(1 ~ s) ,s 1 ~ > 1 — s - ps (1 — s)fc 1 — s Deoarece P(l) = 1, prin definiţia funcţiei generatoare p0 = 1 — fcp în aceste expresii p este interpretat ca fiind X/(A;p), adică intensitatea traficului fiecărei faze individuale şi nu este intensitatea traficului total X/pt, corespunzător întregului mecanism de servire Nu este dificil, fiind doar calcule algebrice anevoioase, să determinăm expresia lui pn din funcţia ei generatoare De fapt, coeficientul lui s” dă = (3 2!) unde însumarea trebuie făcută după toate partiţiile lui n, conform relaţiei n = j + ik + m După alte calcule se găseşte numărul mediu al clienţilor în sistem la un moment dat: n = — —• (3 22) 2 1 — (X/p ) k Repartiţia Erlang sau Pearson de tipul III, care reprezintă repartiţia timpului total de servire în cele k faze, are evident valoarea medie 1/p şi dispersia (l/fc)(l/[i2) Prin urmare, dacă coeficientul de variaţie este notat cu Vt, atunci VI = Ijk Aşadar, vom avea şi formula pentru numărul mediu de clienţi din sistem poate fi scrisă astfel: n 1 (1 + VI) 2 1- (X/p) care este exact rezultatul dat de formula Pollaczek-Hincin Este o verificare utilă a rezultatului, mai cu seamă că interpretările legate de semnificaţia lui pn sînt diferite în cele două situaţii în deducerea formulei Pollaczek-Hincin, pn a fost ales ca fiind probabilitatea ca la plecarea unui client, acesta să lase pe altcineva în urma sa în situaţia noastră, pn a fost considerat a fi proporţia timpului pentru care vor fi prezenţi n clienţi în sistem (X/ţx) , 1 - (X/fi) ’ (X/ţx) , 1 - (X/fi) ’Timpul mediu de aşteptare este dat de Este limpede că în condiţiile în care Ic tinde la infinit, va tinde către zero, iar formulele vor da rezultate pentru timpii de servire constanţi Repartiţia timpului de aşteptare poate fi obţinută considerînd că n reprezintă numărul fazelor existente în sistem şi că fiecare fază este exponenţială Timpul total de aşteptare pentru ultimul client, considerînd ca disciplină de aşteptare regula FIFO, este suma a n variabile exponenţiale, fiecare avînd media 1/p Deci transformata Laplace a repartiţiei timpului de aşteptare este dată de relaţia Reamintind că P(s) este funcţia generatoare a lui pn, aceasta este echivalentă cu Momentele repartiţiei timpului de aşteptare pot fi calculate din această expresie, prin derivare repetată, în raport cu z şi făcînd z = 1 în rezultate Această inversare a transformatei Laplace ar fi de un risc enorm, chiar folosind tehnicile aproximative De aceea, repartiţia timpului de aşteptare nu este uşor de obţinut, exceptînd cazurile în care numărul fazelor este relativ mic (să spunem, k= 2 sau 3)* 3 4 MODELE CU TIMPUL DE SERVIRE CONSTANT Min/1 : (oo/FIFO) în capitolul următor ne vom ocupa sumar de analizele făcute de Crommelin asupra modelului M/D/c : (oo/FIFO); fiind un caz special al acestui model, anume pentru c = 1, nu vom insista îna demonstra rezultatele ce vor urma Vezi Morse, pp 72 — 82 Funcţia generatoare pentru acest model cu o staţie de servire este (3 25) Prin derivarea de n ori în raport cu «, împărţind apoi prin n şi făcînd s = 0, obţinem (3 26) (3 26) pentru n >2 pentru n >2 p0 = i - p; Pi = (i - p)(tp - i); P, = (1 - P) X (-l)“-fceÂp P, = (1 - P) X (-l)“-fceÂp(k?y-k ţip)-1- (n — Ic)! (n—k—1)! Aceasta este forma generală a soluţiei stării stabile Prin derivarea în raport cu s a funcţiei generatoare, punînd apoi s = 1, obţinem că numărul mediu de clienţi în sistem este n = 1 + - —• 2 1 — p Deci timpul mediu de aşteptare în şir şi timpul mediu total de aştep- tare pot fi deduse uşor : w,- -1- p , w 1 - 1 2(i 1 —p 2(i 1 - p Repartiţia timpului de aşteptare poate fi obţinută şi în forma W5(ţ) = (1 -p) X e^’-*’ ~ (3 27) i“o i' Acest model este deosebit de util, deoarece există numeroase situaţii reale în care coeficientul de variaţie a repartiţiei timpului de servire este destul de mic, astfel încît timpul de servire poate fi considerat ca fiind constant în cadrul acestei cărţi vom avea ocazia să studiem cel puţin un astfel de caz M/D/l : (oo/SIRO) După cum arătam, în cadrul multor fenomene cotidiene, regula FIFO a disciplinei de aşteptare nu este respectată; de exemplu în încercarea de a lua maşina în timpul orelor de vîrf, de a cumpăra alimente dintr-un magazin etc în astfel de circumstanţe, ordinea de servire degenerează în aleatoare şi devine posibilă aplicarea regulii SIRO De aceea, modelul prezent care combină timpii constanţi de servire cu servirea în ordine aleatoare devine uneori foarte folositor Dacă nu se cunoaşte cu exactitate disciplina de aşteptare, este indicat a prelucra rezultatele atît de la modelul Jf/D/1 :(oo/FIFO) cit şi de la M/D/l :(oo/SIRO) şi de a le considera ca fiind limitele extreme ale unui rezultat corect Acestea au fost consideraţiile care l-au condus pe Burke să cerceteze acest model Lucrarea sa nu prezintă numai o analiză teoretică ci, de asemenea, un util set de grafice care indică : Probabilităţile |1 — W#(®/6)| pentru valori mai mici decît 1306 Comparaţia timpilor de aşteptare pentru cele două modele Jf/D/1 : (oo/FIFO) şi Jf/D/1 :(oo/SIRO) Comparaţia timpilor de aşteptare pentru cele două modele (Jf/Jf/1: (oo/SIRO) şi Jf/D/1 :(oo/SIRO) înainte de a prezenta analiza făcută de Burke, vom rezuma concluziile acestor grafice în primul rînd, comparaţia 2 indică, aşa cum ne aşteptăm, că regula SIRO determină o repartiţie a timpilor de aşteptare de alt fel decît regula FIFO Există mai mulţi timpi scurţi şi lungi şi mai puţini în mijlocul domeniului repartiţiei (ceea ce este fără îndoială adevărat pentru aceste reguli) Această deosebire merită luată în considerare Ea este foarte mică pentru durate scurte, adică pentru timpi de aşteptare pînă în punctul în care curbele celor două modele se intersectează Dar deosebirea creşte pronunţat pe măsura creşterii timpilor de aşteptare De exemplu, la o intensitate a traficului de 0,7, probabilitatea unui timp de aşteptare egal cu 10 timpi medii de servire este de zero ori mai mare pentru modelul SIRO, faţă de FIFO La o intensitate a traficului de 0,8, probabilitatea unei întîrzieri similare este mai mică de trei ori pentru SIEO faţă de FIFO Pentru p = 0,9, probabilităţile sînt aproape identice Comparaţia 3 făcută de Burke este extrem de interesantă, deoarece ea conduce la o foarte utilă regulă de lucru Pentru orice nivel dat al probabilităţii, timpul corespunzător de aşteptare pentru un timp constant de servire este aproximativ egal cu jumătatea celui determinat pentru modelul cu timp exponenţial de servire Cîteva exemple sînt date în tabelul 3 1 Tabelul 3 1 Comparaţie intre M/M/l : (ao/SIRO) şl M/D/l : (oo/SIRO) Nivelul probabilităţii p v pentru timpi constanţi de servire v pentru timpi exponenţiali de servire 0,1 0,9 12(i 24(i 0,1 0,6 2,5(i 5[i 0,01 0,9 37(i 78(i 0,01 0,6 14(i 0,001 0,8 32(i 66(i 0,001 0,6 12(i 27(i Deducerea repartiţiei timpilor de aşteptare făcută de Burke începe cu ecuaţia Crommelin pentru c = 1 : ĂB 1 e x = (Po + Pi) • — + • • + Pn e"x (n — 1)! El introduce apoi probabilitatea condiţionată q(t\n), ca timpul de aşteptare să nu fie mai mare decît t în cazul în care clientul care aşteaptă este unul din cei n, la momentul următor unei plecări Fie F(t) repartiţia timpilor de aşteptare Atunci F(t) = Pn-!q(t\n) n = l şi problema este de a evalua q(t | n) Intervalele de timp dintre o sosire şi prima plecare din sistem trebuie să fie repartizate uniform deoarece intrările sînt repartizate Poisson şi trebuie să fie independente de intervalele imediat următoare, deoarece disciplina de aşteptare este SIEO în acest caz, timpul total de aşteptare poate fi exprimat ca o sumă de două elemente V = Vj + v2, unde este multiplu întreg al timpului mediu de servire luat ca unitate de timp şi v2 este partea fracţionară rămasă în mod analog, timpul t poate fi divizat în şi t2 Utilizînd independenţele subliniate mai înainte, avem q(t | n) = P(v 1; m> 0), app2(n — 1, m) + (1 — a)pp2(n, »>-l) = (p + w), (n>Q-, m>l) (3 29) Aceste ecuaţii pot fi rezolvate prin metoda cunoscută a determinării funcţiei generatoare, care în acest caz este în două variabile: în £ alocat pentru n şi s pentru m Subliniem că, în mod necesar, deoarece intensităţile medii de servire sînt identice pentru Cj şi C2, avem Ko = 1 — p Mai mult, n-1 s {pi(n - r, r) p2(r, n - r)} = (1 - p)p" T = 0 Calculele algebrice legate de determinarea funcţiilor generatoare sînt foarte anevoioase şi nu vor fi reproduse aici Din funcţia generatoare se obţin următoarele rezultate : (3 30) (3 30)1 +p —ap 1 — ap m = (1 — ct) m = (1 — ct) 1 — ap + ap2 P(1 — P)(l - «Pi 1 — ap + ap2 P(1 — P)(l - «Pi (3 31) (3 31) mQ mQ nQ nQ (1 ~ «)P2 (1 - p)(l - ap)’ (1 ~ «)P2 (1 - p)(l - ap)’ (3 32) (3 33) (3 32) (3 33) (3 34) (3 34) 92 — 92 — X (g —X)(fi — aX) X (g —X)(fi — aX) (3 35) (3 35) X 9 (±(p — aX) Este interesant de remarcat că raportul timpilor medii de aşteptare pentru cele două clase de clienţi este fx -X care este independent de a Nu este surprinzător să constatăm că atunci cînd intensitatea traficului se apropie de 1, şirul clasei C, va deveni extrem de mare, în timp ce pentru clienţii clasei (\ tinde spre limita stării stabile Cox şi Smith analizează în detaliu situaţia în care există k clase de priorităţi Ei arată că în această situaţie : unde a4 este proporţia clienţilor care aparţin clasei C(, este intensitatea medie de servire a acestor clienţi, iar c( este momentul de ordin doi al repartiţiei timpului de servire pentru membrii clasei Ct De asemenea, ei arată că timpul mediu de aşteptare pentru toţi k clienţii este Wq = £ a jWqj i = l Aceste expresii pot fi folosite în evaluarea sistemelor optime de alocare a priorităţilor în mod special, dacă întîrzierilor în şir le pot fi repartizate penalităţi astfel încît valorile lor să fie proporţionale cu timpii medii de aşteptare ai diferitelor clase, atunci priorităţile ar putea fi repartizate în concordanţă cu valorile rapoartelor : timpul mediu de servire penalitatea pe unitate de timp de aşteptare Cu cît valoarea raportului este mai mică, cu atît va fi mai mare prioritatea Aceasta este tot ce am vrut să prezentăm legat de priorităţile absolute Cititorul care doreşte să studieze problema mai profund este sfătuit să consulte lucrările citate în cadrul bibliografiei selectate la sfîrşitul acestei cărţi Recomandăm cu căldură explicaţiile date în cartea lui Cox şi Smith 3 5 2 PRIORITĂŢI RELATIVE Această secţiune va începe cu un comentariu personal Exceptînd electronica şi telecomunicaţiile, nu am întîlnit situaţii în care să fi fost utilizate priorităţile relative sau să fi existat posibilitatea utilizării lor, chiar atunci cînd acest lucru ar fi fost de dorit Nu cred că priorităţile relative să aibă aplicaţii în procese de aşteptare în care clienţii sînt fiinţe umane Deoarece personal am avut puţine ocazii în care am întîlnit priorităţi relative, s-ar putea să fie foarte dificil pentru cititor să-şi imagineze situaţii în care pot să apară O aplicaţie ma joră este legată de utilizarea sistemelor mari de calculatoare electronice digitale care operează on-line Pentru sistemele on-line nu este deloc neobişnuit să aibă de rezolvat o mare varietate de probleme Unele dintre ele sînt mai importante decît altele şi sînt repartizate priorităţi relative, astfel îneît dacă o problemă prioritară necesită o prelucrare imediată, o rutină specializată a sistemului de operare memorează programul care era rulat în acel moment, întrerupîndu-1, şi-l încarcă pentru execuţie pe cel cu prioritate mai mare Această situaţie este foarte utilizată la sistemele electronice on-line, destinate rezervărilor de bilete pentru cursele aeriene în astfel de cazuri, după ce o problemă cu prioritate înaltă a fost prelucrată, programul întrerupt revine ca şi cum acum prioritatea sa ar fi cea mai înaltă, dar imediat prioritatea sa reintră în nivelul normal în cadrul sistemelor cu prioritate relativă clienţii cu maximă prioritate nu sînt afectaţi deloc de cei cu priorităţi mai scăzute Să presupunem de exemplu că, asemănător secţiunii precedente, există două clase prioritare C± şi C2, în proporţiile a şi 1 — a ale clienţilor Timpii lor de servire sînt repartizaţi exponenţial, cu parametrii [ij şi ;x2 Atunci, probabilitatea, în stare stabilă, ca să existe n clienţi din C1 în sistem este Pin = p”(l pj), n = 0,1, 2, , unde pj = aX/fip Ecuaţiile stării stabile pot fi scrise ca de obicei pentru probabilităţile p(n, m) ca să existe n clienţi ai clasei prioritare C1 şi m din C2 în sistem Funcţia generatoare poate fi dedusă destul de uşor, iar apoi prin derivare poate fi arătat că numărul mediu al clienţilor din C2 este P2 ( 1 - Pl-P2\ P2 ( 1 - Pl-P2\ 1 + Pk 1-P1 p-1) 1 + Pk 1-P1 p-1) De asemenea, este posibil a deduce timpii medii de aşteptare şi timpii medii totali de aşteptare în felul următor : W1 = ! 4 L 1 ; 1Fz = - -1- ( 1 4- A ; 1 - P “1 “zl1 - Pi - Pz) V 1 — Pi P-l / Pi “i(l - Pi) Pi “i(l - Pi) “z(l “z(l 1 Pi ~ Pz) 1 Pi ~ Pz) 3 6 COMENTARIU GENERAL Există un model foarte general al proceselor cu un şir de aşteptare, cu intrări şi servire generale şi disciplina de aşteptare FIFO, elaborat de W L Smith Este posibil — cu unele eforturi însă — a deduce rezultatele FIFO determinate în cadrul acestui capitol din formulele generale ale lui Smith Nu există însă asemenea analize pentru modelele cu mai multe şiruri de aşteptare Cititorul, care este familiarizat cu transformata Laplace şi teoria funcţiilor de variabilă complexă, este călduros sfătuit să citească lucrarea lui Smith, în cazul fericit că ar putea-o procura MODELE CU MAI MULTE STAŢII DE SERVIRE 4 1 MODELE Diferenţa esenţială dintre sistemele cu o staţie de servire şi cele cu mai multe staţii de servire este aceea că ultimele prezintă o mai mare varietate a disciplinelor de aşteptare; trebuie ţinut seama mai mult de posibilitatea unei comportări neobişnuite a unui client Nu este deloc uşor de a formula reguli matematice în cazul modelelor cu tipurile cele mai amorfe de disciplină a aşteptării Pînă în prezent s-au făcut foarte puţine încercări în acest sens Modelele clasice cu mai multe staţii de servire presupun de obicei regula FIFO şi mai rar regulile LIFO sau SIRO în orice caz s-a presupus principial că toţi clienţii sosiţi vor forma un singur şir de aşteptare şi că acest şir va alimenta cu clienţi staţiile de servire, după o anumită regulă bine definită de tipul FIFO, LIFO sau SIRO Posibilitatea dezintegrării — unicul şir să se împartă în mai multe şiruri, poate chiar cite unul pentru fiecare staţie de servire — precum şi fenomenele asociate, ca cel de deplasare, servire preferenţială etc , nu au fost luate în considerare Chiar şi în aceste situaţii, modelele cu mai multe staţii de servire sînt mult mai greu de analizat din punct de vedere matematic decît cele cu o singură staţie Din această cauză au fost cercetate foarte puţin şi dispunem de foarte puţine modele utilizabile 4 1 1 M/M/c : (oo/EIFO) Acesta este cel mai simplu model cu mai multe staţii de servire şi de aceea îl vom prezenta la început Se presupune că sosirea clien- ţilor are loc după o repartiţie Poisson, cu o intensitate medie de X pe minut La sosire, fiecare client are trei opţiuni de comportament: Dacă toate staţiile de servire sînt ocupate, se va duce la unicul şir de la care se alimentează toate staţiile Dacă o staţie este liberă şi toate celelalte ocupate, se va duce imediat la cea liberă Dacă există două sau mai multe staţii libere, se va duce la oricare din ele, fără vreo preferinţă Pe măsură ce staţiile se eliberează, clientul care a stat cel mai mult în aşteptare se duce de îndată la staţia eliberată; dacă nu există clienţi în aşteptare, evident că staţia va rămîne liberă pînă la o nouă sosire, iar atunci, conform cu regulile descrise, clientul nou sosit va putea alege serviciul acesteia Nu există posibilitatea ca o staţie să fie liberă dacă sînt clienţi în aşteptare Toate staţiile realizează acelaşi fel de servicii în cadrul fiecărei staţii, repartiţia timpului de servire este exponenţială cu intensitatea b, astfel încît intensitatea medie a serviciului este de p clienţi pe minut Ca de obicei, vom nota cu pn(t) probabilitatea prezenţei an clienţi în sistem la momentul t Cu această notaţie este uşor să deducem ecuaţiile diferenţiale pentru modelul pe care-1 analizăm, urmînd în linii mari raţionamentul utilizat în capitolul precedent pentru modelul Jf/Jf/1 :(oo/FIFO) Acestea sînt; dZ dZ= - Xp(O, t) + p p(l, t), d#(n —= — (x + nji)p(n, t) + Xp(n — 1, t) + dZ (4 1) (4 1)+ (n + 1) u p(n + 1, t); l C- c 1 în aceste formule p = X/(cg), adică traficul măsurat în erlangi împărţit cu numărul staţiilor de servire Traficul obţinut este uneori notat cu e, iar formulele de mai sus sînt atunci scrise înlocuindu-1 pe p cu e/c (Excelenta carte a lui Fry utilizează această notaţie ) Acest model presupune că nu există limitare a clienţilor care pot exista în sistem la orice moment dat; altfel spus, sistemul are o capacitate infinită Aceasta este o simplificare ce poate fi înlăturată cu uşurinţă, iar modelul rezultat este mai general şi mai folositor Să presupunem că prima din cele trei reguli legate de comportamentul clientului este înlocuită de următoarea : la Dacă la sosirea unui client toate staţiile de servire sînt ocupate, acesta se va aşeza la unicul şir de aşteptare al tuturor staţiilor dacă numărul clienţilor existenţi în sistem este mai mic decît L, iar dacă acest număr este L, el nu se va mai aşeza în şir şi va părăsi sistemul Vom numi aceasta, restricţia asupra capacităţii sălii de aşteptare Capacitatea maximă a sistemului este L, iar capacitatea maximă a sălii de aşteptare este L — c Clienţii care sosesc şi, găsind sistemul supraaglomerat, îl părăsesc, sînt numiţi clienţii respinşi sau, în cazul sistemului telefonic, convorbiri pierdute Pentru acest model, ecuaţiile stării stabile sînt : (Ă + + (n + 1) = 0, (l c Ecuaţiile stării stabile sînt: (l-8omj >P (n-l ; m1? , - 1) 4- + («h + P(n + 1; ml + 1, m2, mk) + A'-l + X (1-80mr) (mr+1+l)k[Lp(n-, mk, , mT—1, mr+1 + l, , w„) - î=i — (X + nky ) p{n ; mk, , mk) = 0 (c > n = S mr), (1 — 80mfc) [Xp(c—1; , mk — 1) + + («h + 1) k[Lp(c + 1; «h +1, m2, ,mk — 1)] + fc-i + S (1 —^omr) (mr+1-l) ktfp — l,mTJrk +1, , mk) - T = 1 —(X + ck\i)p(n ; m1, mk) = 0 (c = n = S mr), ~>P (n - 1; mlf , mk) + + (1 — S0m4) («ii+ 1) ky d(n + 1; m1 + 1, m2, , mk l, mk-l) + fc-1 + £ (1 — 80mr) (mr+1 +1) kp p (n mr—1, mr+1 + l, , mk) — r = l — (X + cky )p (n ; mv , mk) =0 (n > c = S M'r) Morse a rezolvat algebric aceste ecuaţii pentru c = 2, k = 2 Analiza este excesiv de laborioasă şi se pare că este practic imposibil de a trata astfel cazuri mai ample Pentru rezolvarea unui sistem pentru care modelul M/E^c : (oo/FIFO) este cel mai adecvat, este necesar a se aplica unul din următoarele trei procedee : Utilizarea calculatorului electronic pentru a rezolva ecuaţiile stării stabile, pentru cazuri numerice specificate Simularea procesului (absolut necesară dacă se cer informaţii legate de stările de tranziţie — vezi cap 6) Utilizarea unui model aproximativ (vezi cap 8 şi 10) Nu există tabele sau grafice capabile să furnizeze informaţii legate de timpul total de aşteptare sau de timpul de aşteptare pentru astfel de modele, exceptînd cazuri foarte rare, pentru care c = 2 sau c = 3 MID/c : (co/I-Il'O) Acest model a fost studiat pentru prima dată în 1930 de către Crommelin Au fost obţinute rezultate doar pentru starea stabilă Fie pn probabilitatea existenţei a n clienţi în sistem, aşteptînd sau fiind serviţi, şi fie an probabilitatea cumulată : Putem calcula timpul în multipli ai timpului de servire constant, care el însuşi poate fi considerat ca fiind unitatea de timp Astfel putem lega numărul clienţilor prezenţi la începutul unui interval de timp de cel al clienţilor de la sfîrşitul intervalului : Po = ac e"\ (Faptul că la începutul intervalului de timp sînt mai puţini clienţi decît staţii de servire şi faptul că în intervalul de timp considerat nu înregistrează sosiri implică faptul că la sfîrşitul intervalului nu v fi clienţi pentru servit ) Pi = cc Xe z + pc+] e x, X” X”-1 Pn = « , 1—sc eX(1 ” (4 9) (s 1) (g SjJ (s &C-1)> , 1—sc eX(1 ” (4 9) în care 1, s1( s2, • • • > sc-i sînt cele c zerouri ale expresiei situate pe cercul unitate din planul complex Probabilitatea pn poate fi obţinută prin derivarea succesivă a funcţiei P(s) şi punînd apoi s = 0 : Pn Pn Timpul mediu de aşteptare poate fi dedus calculînd pe P'(l) în maniera obişnuită Rezultatul obţinut va fi £ y + X2 ~ c (1 ~ t) ca timpul de aşteptare să depăşească durata t Aceasta are o expresie foarte complexă şi practic inutilizabilă, exceptînd cazurile în care se utilizează calculatorul electronic Din fericire, pot fi găsite grafice foarte utile pentru c = 1, 2,3, 4, 5, 8,10, 20 şi pentru p cr Jackson a dedus soluţia stării stabile a acestui model Sarcina determinării soluţiilor dependente de timp este considerabilă (şi, din punct de vedere practic, neesenţială) Ecuaţiile stării stabile, conform notaţiilor lui Jackson, sînt următoarele : Â + £ â (nr) a(nr) gj p(nn n2, , nk) = k = £ 8(nr + 1) a(nr + 1) y r p(nlf nr + 1, nr+1 - 1, , nk) + r=l + Ap («i-l, «2, «3, • • •, nk)- (5 1) Trebuie subînţeles faptul că dacă un argument de forma nr+1 — 1, de exemplu, devine negativ, atunci probabilitatea corespunzătoare poate fi ignorată, fiind nulă Mai mult, termenul cu —1 şi acela cu 4-1 pot fi neglijaţi în termenul final al ultimei sume Simbolurile speciale folosite au următoarele semnificaţii: a(nr) = \n” 8(nr) = P’ lo, 8 K+1) = 1 Aceste ecuaţii sînt destul de dificil de interpretat şi de aceea, înainte de a trece la prezentarea soluţiei lui Jackson, cititorului îi este recomandat să urmărească aceste ecuaţii pentru k = 2 faze Ele iau forma : Xp(0, 0) = |X2P(O, 1); (X-|-p x) p(l, 0) = Xp(0, 0) + |X2P(1, 1); (X + H2) p(0, 1) = 0) 2p 2 p(0, 2); (X 4- Hi + 1) = W, 1) + HiP(2, 0) + 2[x2p(l, 2); (5-2) (X+c1[i 1+n2ti 2)p(n1, na)=Xp(nx—1, n2)4-cx p xp(nx4-l, «2-l)+ +(n24-l)}x2î>(n1, n2+l) (X 4- nx[ix 4- ^a) P(nu wa) = *P(nj ~ 1, ^2) + 4- (nx 4-1) 4-1, n2 - 1) + 4" CiPaPfol, W24"l) w2^c2) ; (X4-cxnx 4- C2H2) P(ni, = *P(ni — w2) 4- C1H1 P(ni + 1, n2—1) 4- 4-c2[i2p(nx,n24-l) (n>cx; «2>c2) Aceste ecuaţii sînt într-adevăr destul de dificile Totuşi soluţiile pentru n2) sînt relativ simple Să ne întoarcem însă la modelul general al lui Jackson, cu k faze Soluţia este (5 3) (5 3)k p(nY, nz, nk) = p(0, 0, , 0) JJ b(nr), r-1 în care b(nr) = Cititorul care va face o comparaţie cu capitolul anterior va găsi aceste expresii asemănătoare cu soluţiile lui pn pentru modelul Jf/Jf/c: (oo/FIFO) Termenul p(Q, 0, , 0) se determină prin însumare, după n, deoarece suma tuturor termenilor este egală cu 1 Deci ao x oo S I-I m1 = 0 n £ = 0 nk = (5 4) -n[£ 5(%)l=n^ 0 = 1,2, ,*) r=lLi = O J j = l Din nou, o vizualizare a formei acestei soluţii pentru k = 2 ar putea fi utilă De exemplu, pentru şi n2 (n^ • p(2> (n2), (5 6) în care p(i)(n>) pentru j = 1 sau 2 sînt soluţiile corespunzătoare pentru cele două modele de aşteptare A/’(X)/Af(pa)/c1 : (oo/GD) şi respectiv Jf(X)/Jf(p 2)/c2: (oo/GD), după cum se constată imediat (vezi capitolul anterior), adică soluţia găsită cu mulţi ani în urmă de către Karl Erlang pentru modelele cu o singură fază Ca implicaţie a acestei situaţii, ieşirile primei faze trebuie să fie Poisson cu media X Astfel spus, trebuie să fie, din punct de vedere statistic, de aceeaşi formă ca şi intrările Acest lucru va permite cu siguranţă celor două faze să se comporte ca două şiruri independente De fapt, ieşirile modelului M/M/c : (oo/GD) sînt aleatoare, mai precis Poisson, cu aceeaşi medie ca şi intrările Există o demonstraţie foarte riguroasă a acestei informaţii într-un articol al lui Burke, pe care însă nu o vom prezenta aici Este absolut necesar a sublinia utilitatea practică a acestui fapt: dacă nu ar fi fost adevărat, atunci multe sisteme de aşteptare care pot fi studiate (cel puţin la un nivel rezonabil de aproximare) cu ajutorul modelelor independente ale sistemelor mai simple ar fi, din punct de vedere matematic, de nerezolvat Simularea pe calculatorul electronic ar deveni în acest caz absolut necesară Cititorul poate observa că se pot petrece lucruri foarte stranii în cazul în care capacitatea şirului primei faze a unui model serial, să zicem cu două faze, este limitată; de exemplu, dacă prima fază este reprezentată de modelul M/JUI^ : (Z/GD) De fapt, o anumită proporţie a traficului primei faze, cu o intensitate de X(Ci(ii) este înlăturată, datorită sosirii unor clienţi atunci cînd capacitatea sistemului, L, este completată De aceea, traficul care trece prin a doua fază este oarecum diminuat Problema care se pune este următoarea : care este repartiţia intervalelor dintre plecările succesive care constituie ieşirile filtrate de la prima fază? Ca şi în multe alte cazuri, răspunsul la această întrebare poate fi dat prin intermediul următorului exemplu forarte simplu Să considerăm modelul Jf/Jf/1 : (1/GD) — disciplina de aşteptare este de fapt nesemnificativă pentru această situaţie — în care sosirile au loc cu o intensitate medie X iar timpul mediu de servire este 1/p j Este limpede că atunci cînd un client este în curs de servire, toate celelalte sosiri sînt respinse Prin urmare, intervalele dintre plecările consecutive din sistem, aşa-numiţii timpi dintre plecări, sînt fiecare suma a cîte două componente: tx — timpul scurs între o plecare din sistem şi o nouă sosire în sistem; i2 — timpul necesar servirii noului sosit Deoarece repartiţia timpilor dintre sosiri este exponenţială, tv prima din aceste componente, este total independentă de momentul în sare clientul anterior — care poate a fost respins — a sosit, ca şi de momentul celei mai recente plecări (Reamintiţi-vă „uitarea” caracteristică repartiţiei exponenţiale, menţionată în cap 2 ) Prin urmare, din punct de vedere statistic, prima componentă are aceeaşi repatriţie cu cea a timpilor dintre sosiri, adică este exponenţială cu media 1/X A doua, prin definiţie, este de asemenea repartizată exponenţial Ca urmare, timpii dintre plecările din sistem sînt repartizaţi ca suma a două variabile exponenţiale, astfel încît repartiţia lor nu poate fi exponenţială Acest raţionament referitor la faptul că ieşirile din- tr-un proces Jf/Jf/1: (1/GD) nu sînt Poisson a fost făcută de către P G Finch, care a demonstrat că : Intervalele dintre plecări sînt repartizate exponenţial doar atunci cînd capacitatea sistemului este infinită sau, utilizînd notaţiile noastre , L = oo Mai general, o capacitate infinită a sistemului (L = oo) şi o repartiţie exponenţială a timpilor deservire sînt condiţii necesare şi suficiente pentru ca intervalele dintre plecări să fie independente De aceea, dacă aşteptarea din prima fază a unui proces format din două faze este limitată, atunci cele două faze nu sînt independente, chiar dacă între faze este permisă o aşteptare nelimitată (acest rezultat se menţine şi în cazul existenţei mai multor şiruri paralele, pentru fiecare fază) Intrările în cea dea doua fază nu sînt repartizate Poisson în discuţia ce urmează, despre fenomenul de blocare, vom arăta că cele două faze nu sînt independente dacă reciproca este adevărată ; adică, în caz că avem o capacitate infinită a sistemului în prima fază şi o capacitate limitată pentru aşteptarea dintre faze Este important să subliniem că restricţiile făcute capacităţii sistemului sau, altfel spus, a spaţiului de aşteptare la orice fază a unui sistem serial multiplu, determină dependenţa fazelor în practică însă, de obicei se acceptă ipoteza că fazele unui sistem, să spunem, bifazic, sînt independente, în cazul în care este admisă o capacitate suficient de mare înaintea fiecărei faze; o metodă practică de a decide ce înseamnă o capacitate suficient de mare este aceea de a calcula probabilitatea de refuzare a unei sosiri O probabilitate mai mică decît 0,10 înseamnă practic capacitate suficientă Cu cît este mai mică probabilitatea refuzului unei sosiri, cu atît încrederea în independenţa fazelor este mai mare (în cap 15 vom face referiri mai largi legate de aceste aspecte ) Să ne întoarcem la sistemul de aşteptare Jf/Jf/1 : (1/GD) Ar fi util să ni-1 închipuim ca fiind prima fază a unui sistem bifazic, în care faza a doua este formată dintr-o staţie exponenţială cu timpul mediu de servire 1/(a2, cu o capacitate infinită, şi să ne imaginăm condiţiile pentru care şi acest sistem simplu poate fi considerat ca fiind format din două faze independente Există două cazuri speciale de acest fel, ambele extreme Primul corespunde unei intensităţi a traficului p = X/(xx foarte scăzute pentru prima fază şi uneia moderat de mare pentru a doua Deci nu are loc aproape nici o respingere a vreunui client în prima fază şi, cum nici unul nu va fi nevoit să aştepte, repartiţia timpului dintre plecări din prima fază este identică, din punct de vedere statistic, cu repartiţia timpului dintre sosiri Efectul primei faze de servire este acela de a determina timpi exponenţiali de sosire şi nu de a modifica repartiţia timpilor dintre sosiri în acest caz, faza secundă poate fi tratată ca un sistem lf/Jf/1 : (oo/GD) Astfel de situaţii, de obicei, au mai multe staţii paralele pentru fiecare fază, pot fi uneori întîlnite cînd prima fază este o operaţie de servire de genul „agentului de circulaţie” Exemple de acest gen apar în telecomunicaţii şi în restaurante Al doilea caz special corespunde suprasaturării primei faze, adică cazului în care intensitatea traficului depăşeşte unitatea în ciuda proporţiei ridicate a respingerilor, prima fază este aproape permanent ocupată, astfel încît repartiţia intervalelor dintre plecări este, din punct de vedere statistic, aproape identică cu repartiţia timpului de servire Dacă aceasta din urmă este, aşa cum presupunem, exponenţială, înseamnă că faza a doua are intrările Poisson, şi poate fi considerată independentă de prima Asemenea situaţii pot fi remarcate în unele birouri de la companiile aeriene, în care apelurile telefonice au o frecvenţă maximă pentru perioadele în care vremea este foarte nefavorabilă Un alt aspect interesant al proceselor seriale este fenomenul de blocare Acesta are loc doar atunci cînd se admite un şir de aşteptare limitat în faza a doua, iar clienţii nu pot părăsi prima fază fără a trece prin a doua De obicei, lucrările despre teoria aşteptării acceptă restricţia prin care nu este admis nici un şir de aşteptare în faza a doua, dar această simplificare, cu toate avantajele ce le prezintă pentru calcule, nu este necesară De exemplu, să presupunem, că fiecare fază este formată dintr-un şir individual de aşteptare şi a doua fază poate cuprinde doar L2 clienţi; în momentul cînd faza a doua este completată, iar un client a fost servit în prima fază, prima fază este blocată în aceste condiţii, valorile date, lui L2 apar fără însemnătate Dar cu cît L2 este mai mare, cu atît probabilitatea ca blocarea să se producă este mai mică Uneori se presupune capacitatea primei faze L1 ca fiind egală cu 1, astfel încît dacă un client soseşte şi o găseşte ocupată, acesta să fie respins Această presupunere poate fi însă extinsă şi pentru alte valori ale lui Lr Totuşi, cu cît Lx este mai mare, cu atît respingerile vor fi mai puţine Pentru exemplificare, vom analiza cazul convenţional a două faze exponenţiale, de capacitate infinită în prima fază şi 1 în a doua Clientul care doreşte să părăsească prima fază pentru a intra într-a doua, în momentul cînd staţia corespunzătoare fazei secunde este completă, va trebui să rămînă blocat în prima fază Vom defini următoarele probabilităţi ale stării stabile : p(n, 0) — există n clienţi în prima fază, unul care este servit şi n—1 aşteptînd; nimeni nu se găseşte în faza a doua; p(n, 1) — n clienţi în prima fază, unul care este servit şi n — 1 aşt eptînd ; faza a doua completă; pb(n, 1) — n clienţi în prima fază, unul care a fost servit şi blochează staţia de servire, în timp ce n — 1 clienţi aşteaptă Vom presupune că pb(n, 1) se schimbă instantaneu în p(n — 1,1) în momentul eliberării fazei secunde Ecuaţiile stării stabile sînt: Xp (0, 0) = (i2p(0,1); (X + V^)P (0,1) = P-1P(1, 0) + (XjPfiCl, 1); (X + n2) Ps(l, 1) = |X1 p(l, 1); (x + y-i)p(n, 0) = \p(n — 1, 0) + p 2p(n, 1); n>0; (5 10) (X+p 1 + (i2)î’(»> 1) = Xp(n — 1,0) + p !p(n + 1, 0) + + y 2pb(n + 1,1); n> 0; (X+^Lg) pb(n, 1) = Xp6(n—1,1) + p !P(n, 1); n > 1 Aceste ecuaţii pot fi rezolvate în maniera cunoscută, calculul algebric fiind însă destul de anevoios Elementul de interes este acela că blocarea ce apare reduce intensitatea de servire efectivă a primei faze sub valoarea ; de aceea, volumul maxim de trafic ce poate fi păstrat în echilibru statistic este mai mic decît Z/ixj A fost demonstrat că pentru exemplul dat, acest volum este P ClX, + IX,) (5 U) [if + fJ-iP-2 + (X-3 în cazul particular, cînd p 1= p,2, s™ax devine egal cu 2/3 Este fără îndoială o valoare mai mică decît 1 — rj, unde tj este o cantitate pozitivă, valoare ce rezultă dintr-o stare stabilă în care nu apare blocarea Aceasta indică faptul că surplusul de capacitate dintre două faze a sistemului de aşteptare va permite o utilizare superioară a primei faze Se poate arăta că un surplus de capacitate de numai o unitate măreşte pe e“:,x la valoarea 3/4 Modelele seriale nu pot fi uşor analizate, exceptînd cazul în care repartiţiile timpilor de servire sînt exponenţiale Motivul este că dacă aceste repartiţii nu sînt exponenţiale, atunci intervalele dintre plecările consecutive din oricare fază nu sînt independente De obicei însă, ele nu sînt exponenţiale Cele două exemple pe care le-am examinat şi în care şirul de aşteptare a fost limitat fie în prima fază, fie în a doua (cu rezultate diferite), arată că problema se complică foarte mult în astfel de situaţii Se pare că tehnicile de simulare constituie singura cale de determinare a proprietăţilor unui model de sistem de aşteptare, serial de orice grad de complexitate ar fi Este posibil a se realiza unele succese în multe aplicaţii practice făcînd presupunerea că faza a doua a unui sistem, să zicem bifazic, este independentă de prima Dacă prima acţionează ca un filtru din cauza unei capacităţi limitate, va trebui corectată intensitatea medie a intrărilor în faza a doua 5 2 SISTEME CICLICE Cu toate că existau de mulţi ani analize legate de sistemele ciclice de aşteptare — cum ar fi analiza făcută de Jackson şi Taylor — principala contribuţie teoretică a fost adusă de către E Koenigsberg Vom încheia acest capitol printr-o succintă prezentare a modelului lui Koenigsberg* Modelul este format din fc staţii individuale, legate între ele astfel încît formează un circuit închis, ieşirile ultimei faze constituind intrările primei Toate staţiile prezintă o repartiţie exponenţială a timpilor de servire, iar intensitatea medie de servire pentru faza i este [! x2 > x3> > xn , subliniază necesitatea găsirii unor măsuri ale eficacităţii sau a unor funcţii obiectiv şi a optimizării lor în cercetarea operaţională, soluţie începe să însemne optimizarea unui model bine determinat Nemulţumirea analiştilor de cercetare operaţională constă în aceea că în general executanţii lucrărilor nu sînt foarte receptivi la indicaţiile ce li se dau, ba uneori nici nu le dau prea mare atenţie Pentru a explica această stare de lucruri au fost formulate mai multe motive, s-a enunţat o întreagă gamă de opinii, de la autosubapreciere („riu ştim cum să vindem produsele d-voastră”) pînă la abuzivul „conducătorii şi executanţii sînt prea incapabili şi reacţionari” Cu toate că ar fi greşit să-i considerăm pe toţi analiştii de cercetare operaţională de talia liii Albert Laskers sau că toţi conducătorii sînt de talia lui Alfred Sloan, nu acestea sînt cazurile unui fenomen atît de extins, aproape universal Realitatea constă în existenţa unei divergenţe de opinii dintre cercetători şi executanţi, în ceea ce priveşte' înţelesul noţiunii de soluţie în timp ce analiştii cercetării operaţionale se plîngeau că executanţii nu dau atenţie sugestiilor lor, economiştii au avut nemulţumiri asemănătoare, că oamenii de afaceri nu le dau atenţie în ultimii ani, mulţi economişti, mai cu seamă din Statele Unite, au căutat să descopere adevăratele cauze ale neînţelegerilor iscate între ei şi oamenii de afaceri Cu toate că ei continuă să se plîngă că oamenii de afaceri ignoră teoria economică (ceea ce ar putea fi adevărat), economiştii au încetat de a-i mai considera incapabili (ceea ce nu întotdeauna este adevărat) Economiştii au ajuns la două concluzii: prima este că tendinţa academică de a analiza problema, prin intermediul unor modele bine definite, nu este împărtăşită de către oamenii de afaceri; a doua, arată că oamenii de afaceri — chiar şi antreprenorii — nu acţionează şi nu pot acţiona consecvent conform unei singure reguli simpliste, aşa cum este maximizarea beneficiilor Prima concluzie înseamnă următoarele : cei ce iau decizii în administraţie, forţele armate, în afaceri şi industrie nu gîndesc- procesele operaţionale şi tendinţele lor prin intermediul unor modele bine definite : imaginea lor este cea pe care o ai cînd te uiţi prin sticlă mată A doua concluzie arată că executanţii nu acţionează ca pentru a atinge un punct final dorit; cu alte cuvinte, să maximizeze funcţia obiectiv aşa cum este beneficiul A fost lansată ideea că ei ar acţiona de o asemenea manieră, încît să asigure cît mai mult posibil supravieţuirea întreprinderii Baumol a mers mai departe cu această teorie, sesizînd că ei acţionează ca pentru a maximiza beneficiul cu condiţia satisfacerii anumitor restricţii (de exemplu, rata de amortizare a investiţiei să nu fie sub 12%) Aceasta este o idee excelentă, ce şi-a găsit deja locul printre principiile şi metodele cercetării operaţionale, prin intermediul programării liniare în spatele ei stă o impunătoare dezvoltare matematică, ce a început de la Lagrange Se mai poate spune că obiectivele oamenilor de afaceri sînt calei- doscopice Ei au foarte multe scopuri, unele din ele chiar contradictorii, iar balanţa importanţei unuia sau altuia este în permanentă schimbare Faptul că Baumol vorbeşte de o maximizare cu restricţii implică existenţa mai multor factori în cadrul unei situaţii ce au tendinţa de a se opune unul altuia Dacă toţi cresc sau descresc împreună, nu ar mai fi necesară optimizarea prin introducerea unor restricţii Şi totuşi, care este legătura dintre aceste lucruri, cercetarea operaţională şi soluţia problemelor operaţionale? Argumentele aduse de economişti indică motive pentru inaccep- tabilitatea multor rezultate ale cercetării operaţionale, inclusiv aplicaţii ale teoriei aşteptării Şi aceasta ca urmare în primul rînd a preocupărilor exagerate pe care teoreticienii le au pentru teorie şi precizie atunci cînd prima este imposibilă, iar ultima inutilă şi în al doilea rînd insistenţei asupra găsirii unei funcţii obiectiv care niciodată nu este în perfectă corespondenţă cu scopurile de compromis ale executanţilor Cu alte cuvinte, tehnica de obţinere a soluţiilor prin procedeele clasice de optimizare nu este admisibilă în toate cazurile şi aceasta nu pentru că este o tehnică şubredă, ci datorită variabilităţii caracteristicilor funcţiei obiectiv în scopuri practice modelele proceselor de aşteptare sînt aplicate în trei direcţii principale : în planificarea repartizării noilor facilităţi pentru operaţiile în care aşteptarea, sau aglomerarea, este un factor important în evaluarea strategiilor de desfăşurare a forţei de muncă pentru facilităţile existente, astfel încît să se realizeze condiţii de servire standard şi cu preţ redus Pentru compararea diferitelor forme de organizare a facilităţilor de servire Este neîndoielnic că problemele tehnice ale formării şi analizei modelelor de aşteptare, menite să previzioneze necesarul de capacitate a noilor facilităţi, sînt foarte anevoioase în cap 13 vom analiza o asemenea situaţie Dar modelele de aşteptare nu sînt singura sursă de dificultăţi O mare importanţă trebuie acordată previziunii parametrilor modelelor Sigur că prognoza poate fi uneori hazardată, şi este probabil să se utilizeze una ce nu dă rezultate satisfăcătoare De exemplu, chiar şi cu ajutorul unor tehnici cum ar fi ajustarea exponenţială a seriilor de timp (dar nu încă a modelelor econometrice) se consideră a fi satisfăcătoare o prognoză pentru un an, ce dă traficul lunar de curse aeriene de pasageri dintr-un aeroport, cu o abatere de 5% (în plus sau în minus) Prognozele orare au un grad şi mai mare de nesiguranţă Remarcaţi c-ît de dificil este a determina, să spunem cu cinci ani înainte, costul planificării şi realizării unui nou aeroport O eroare de 25 % în prognozele zborurilor lunare ar putea fi, în acest caz, satisfăcătoare în plus, ar fi dificil de prevăzut exact care va fi timpul mediu de servire, mai cu seamă dacă activităţile ce constituie operaţia de servire sînt noi, sau aranjate într-o formă nouă (în cap 10 a fost ales un exemplu de acest fel ) Chiar şi cel mai competent inginer nu este dotat cu puteri supranaturale! în aceste condiţii căutarea unui model de aşteptare elegant şi riguros este greşită Inevitabil, la sfîrşitul lungilor şi dificilelor serii de calcule ale cerinţelor viitoare, în care precizia depinde direct de volumul de date de previzionat, de durata de studiu şi de conştiinţă profesională, inginerii consultanţi sau arhitecţii, bazîndu-se pe experienţa lor, vor adăuga cel puţin 10 % la fiecare calcul Cu alte cuvinte, din punct de vedere practic, este foarte neeconomic să se determine conştiincios şi exact valorile strategice ale modelului, în condiţiile unei planificări pe termen mai lung Este însă esenţial să dăm răspuns la întrebarea : care trebuie să fie gradul de precizie al modelului, sau cît putem noi realiza ? A doua clasă de utilizare a modelelor de aşteptare constă din aplicaţiile menite să determine desfăşurarea facilităţilor existente Se remarcă că aceasta reprezintă o situaţie cu totul diferită de cea prezentată anterior Ne putem imagina în acest caz că staţiile de servire există deja şi problema care se pune este de a determina numărul celor ce lucrează în diferite momente ale zilei, în diferite zile ale săptă- mînii etc , pentru a crea un serviciu standard Acest lucru poate fi realizat cu ajutorul prognozei pe termen scurt a volumului de trafic (utilizînd traficul în sens Erlang) pusă la înde- mîna managerului Trebuie totuşi să menţionăm că de obicei un manager preferă să utilizeze formule foarte simple care să-l ghideze printre multiplele sarcini care le are de îndeplinit în acest caz, simplitatea înseamnă aproximare în cap 15 se va analiza un astfel de caz Un model cu un înalt grad de fidelitate în reprezentarea realităţii va fi foarte complex, cu grad mare de precizie şi deci aproape inutil Un model cu un grad scăzut de fidelitate în reprezentarea realităţii poate fi mai simplu, nu tot atît de precis, dar util Cititorul va decide care ar trebui să fie alegerea făcută de către un analist de cercetare operaţională în sfîrşit, au rămas problemele de comportare a diferitelor forme de organizare a facilităţilor de servire Să presupunem că există trei posibilităţi de combinare (Cititorul este rugat să-şi amintească sugestiile d-lui Murdo Macdonald din cap 1 sau să se uite peste exemplul dat în cap 9 ) Să presupunem în continuare că traficul, exprimat în erlangi, este constant şi independent de natura combinaţiei Pentru simplificare, să considerăm că problema este de a decide care variantă produce cel mai scăzut timp mediu de aşteptare, pentru un număr dat de staţii de servire Cîteodată este posibil să se răspundă la această întrebare fără a face nici un calcul, bazîndu-ne doar pe experienţă De exemplu, presupunînd că repartiţia timpului de servire nu este afectată de natura variantei alese, atunci este limpede că o variantă ce va crea o utilizare deplină a staţiilor va fi preferată altei combinaţii ce nu implică acest lucru Sau o variantă care încurajează mobilitatea clienţilor între şirurile de tipul „bilet-ghişeu” va fi preferată unei variante ce nu realizează acest lucru, deşi nu este atît de favorabilă ca un sistem „maître d’hotel”, cu posibilitatea maximă de servire, în asemenea condiţii, sînt satisfăcătoare şi concluziile generale rezultate chiar din modelele deduse din ecuaţii integro-diferenţiale Problema ce se pune este de a identifica repartiţia timpului de servire, ordinea sosirilor şi disciplina de aşteptare şi nu de a le calcula precis în alte condiţii nu este suficient să cunoşti care combinaţie este mai eficientă; trebuie cunoscut şi cu cît este mai eficientă o variantă decît celelalte, căci costurile asociate variantelor pot diferi Aceasta introduce un nou element al problemei, care poate fi reprezentat în felul următor : Varlantâ Costul forţei de munc&/om Cheltuieli de lntreţlnere/an Timp mediu de aşteptare/ client 1 45 20 14/|i 2 15 40 1,2/fx 3 10 50 0,9/fi Termenul „cheltuieli de întreţinere” semnifică cheltuielile de achiziţie a utilajului (instalaţiei) şi de menţinere în funcţiune Cheltuielile de întreţinere includ deprecierea şi cheltuielile de menţinere care, pentru orice combinaţie, se presupun a fi constante Tabelul arată comparaţiile realizate pentru timpi de aşteptare egali Variaţia cifrelor din ultima coloană este consecinţa faptului că numărul de anga-, jaţi este o variabilă discretă De aceea, este practic imposibil a egaliza timpii medii de aşteptare prin manevrarea forţei de muncă Remarcăm că, din punct de vedere financiar, a doua variantă este cea mai favorabilă Aceasta este o situaţie clasică Totuşi costul forţei de muncă poate varia ca rezultat al unor negocieri asupra salariului, în timp ce cheltuielile de întreţinere rămîn aceleaşi în acest caz problema este : în ce condiţii de modificări ale costurilor, varianta a doua rămîne cea mai bună ? Cititorul va remarca faptul că varianta secundă continuă a fi cea favorabilă pînă cînd costul forţei de muncă se dublează Probabil că aceasta nu se produce curent, şi dacă este aşa, putem spune că „soluţia” — potrivit căreia varianta secundă este cea mai bună — nu este influenţată de costuri Pot apărea însă şi alte efecte De exemplu, rata creşterii forţei de muncă cu volumul traficului, pentru menţinerea unui serviciu standard, poate diferi de la o variantă la alta Acesta este un fenomen des întîlnit, în special cînd timpii medii de servire diferă între variante (vezi cap 9) Astfel, o creştere de 15 % a volumului traficului poate da următoarele rezultate : Variantă Costul forţei de muncă/an Cheltuielile patronul ul/an Timpul mediu de aeteptare/cllent 1 60 25 1,0/ti 2 25 43 0,8/p 3 13 52 0,9/p Varianta a treia este acum cea mai eficientă „Soluţia” pe care am găsit-o este influenţată de modificările volumului traficului Exemplele prezentate sînt artificiale; dar ele se bazează pe o experienţă suficient de convingătoare Cazurile ce vor fi prezentate în capitolele ce urmează vor ilustra fenomenele posibile Cînd se constată că previziunile pot fi eronate, că restricţiile operaţionale se pot modifica şi că orice presupunere făcută asupra valorilor parametrilor poate fi invalidată de evenimente, este lesne de înţeles că „soluţia” ce trebuie căutată este caracterizată prin robusteţe Punctul de vedere ce a stat la baza acestei argumentări este următorul : este puţin recomandabil a construi modelul unui proces şi prin analizarea sa să se tragă concluzii definitive — ca de exemplu, „portul necesită două noi dane pînă în 1970” — fără a verifica senzi- tivitatea modelului Unele modificări făcute valorilor parametrilor au un efect important asupra intensităţilor traficului în aceste condiţii, nu există distincţii clare ce trebuie făcute între unele răspunsuri şi altele Astfel şi-au făcut apariţia mulţimile vagi (fuzzv) ale bussinesmanului De exemplu, dacă curbele caracteristice costului (cheltuielile forţei de muncă plus cheltuielile de întreţinere) ar fi fost proiectate pentru cele trei sisteme diferite din exemplul nostru în condiţiile unui trafic diferit, pentru fiecare s-ar fi găsit un alt punct de minim Dar, deoarece nu se poate cunoaşte exact dacă coeficienţii funcţiei costurilor au fost corect determinaţi, pot apărea anumite oscilaţii verticale Pentru diferite valori ale acestor coeficienţi, curbele se vor intersecta în puncte diferite, iar minimul se va modifica, ca în exemplul considerat Astfel, sînt introduse şi oscilaţii orizontale în aceste condiţii, este puţin probabil ca modelele foarte sofisticate să servească scopului ce ni-1 propunem O soluţie operaţională a unei probleme de aşteptare este orice soluţie suficient de judicios evaluată pentru a conduce la o decizie corectă Ea poate să nu se bazeze pe un model foarte complicat Ea poate fi dedusă şi din anumite constatări evident neadevărate legate de natura repartiţiei timpului de servire Lucrul cel mai important în cercetarea operaţională nu este de a obţine soluţii în sens matematic, ci soluţii operaţionale, ce vor conduce la elaborarea unor decizii valide, executabile TREI TIPURI DE APROXIMĂRI în cadrul capitolelor precedente am prezentat unele problemt teoretice ; cele ce urmează au legătură directă cu aplicaţiile practice, Acest scurt capitol este o punte de legătură între ele şi este un prilej de a relua argumentarea fundamentală a acestei lucrări în principiu, atît din punct de vedere teoretic cît şi estetic, varianta cea mai dorită de rezolvare a problemelor legate de sistemele de aşteptare constă în a construi un model matematic foarte fidel şi a determina soluţii prin analiza acestui model S-a spus suficient de mult despre dificultăţile matematice ce intervin în analiza modelelor, chiar şi pentru cele mai simple sisteme, astfel încît se pune problema : „Tratarea ideală este oare posibilă?” Examinarea analizelor matematice făcute în capitolele 2 şi 5 pare a duce la concluzia că numai modelele elementare pot fi supuse unei analize matematice exacte şi că aceste modele nu au legătură cu procesele ce se întîlnesc frecvent în lumea reală Cazurile ce vor fi prezentate în capitolele următoare vor demonstra că sistemele cotidiene de aşteptare nu au o structură atît de complicată ca cea a modelelor de aşteptare ce au fost amplu analizate pînă în prezent Putem conchide că tratarea ideală, tradiţionala procedură ultraştiinţifică a cercetării operaţionale, în general, nu poate fi realizată Se pare că reprezentarea sistemelor de aşteptare cu ajutorul modelelor matematice foarte fidele nu poate fi realizată adesea, iar dacă da, atunci modelele nu pot fi analizate profund Aceasta a dus la crearea tehnicilor de simulare Am vorbit în cap 6 despre aceste tehnici pentru a sugera că, folosindu-le, nu se pot obţine chiar atît de uşor pe cît se pare rezultate exacte privind sistemele de aşteptare Apar mai multe dificultăţi: comportarea tranzitorie, mărimea selecţiilor (lucrarea lui Gebhart, prezentată în anexa 1, este recomandată celor care nu cred că aceasta ar reprezenta o dificultate) şi proiectarea experimentală Atît procedura tradiţională cît şi tehnica simulării tind spre a atinge precizia maximă Procedura tradiţională nu este adesea utilizabilă ; simularea însă da, şi se pot obţine răspunsuri suficient de precise, dar cu un cost de calcul relativ ridicat Din păcate, tocmai în astfel de situaţii, care au o mare importanţă practică — ca de pildă, cele legate de intensităţile mari ale traficului — ciclurile de simulare trebuie să fie foarte lungi, pentru a atinge performanţele dorite în capitolul precedent ne întrebăm „Oare este atît de necesară o precizie ridicată F’ Şi răspunsul la această întrebare era că în multe situaţii nu este necesară, iar pentru a indica de ce trebuie să fie aşa a fost introdusă ideea de soluţie operaţională Conceptul de soluţie operaţională nu este analizat in extenso în această lucrare : cititorul va fi capabil să o raporteze la cercetările întreprinse în domeniul proceselor de decizie executivă şi al teoriei deciziei Odată acceptată ideea că a căuta modele exacte şi riguroase, ce vor determina rezultate numerice sau analitice neperturbate, este o greşeală, şi că o aproximare cu un grad acceptabil de precizie este calea de urmat, atunci putem privi în cu totul alt mod atît procedura de construire a modelului matematic cît şi tehnicile de simulare Argumentul acestei lucrări este că teoria aşteptării este utilă, în pofida complexităţii matematice, pe de o parte, şi unor inadvertenţe practice, pe de altă parte, teoria în forma ei actuală poate determina soluţii operaţionale pentru multe situaţii Aceasta cere însă ca să se accepte necesitatea aproximaţiilor de primul tip, adică aproximarea prin utilizarea modelelor foarte simplificate Alt mod de redare a acestei idei este că primul tip de aproximare constă în a obţine soluţii operaţionale prin apelarea modelelor care sînt, apriori, eronate, dar sînt suficient de simple pentru a le putea manipula Capitolele ce urmează, şi în special capitolele 10, 11 şi 12, conţin cazuri practice ce vor arăta aplicabilitatea acestei proceduri Nu o să mai insistăm acum asupra acestor lucruri, exceptînd faptul că pentru un analist de cercetare operaţională procedura ar fi deosebit de utilă, poate cea mai utilă dintre toate în cazul în care aproximarea de primul tip este considerată a fi inacceptabilă, este necesar a utiliza simularea Aceasta ar fi necesară în următoarele condiţii : în cazul că riscul luării unei decizii eronate datorită modelului foarte simplificat nu este cunoscut, sau în cazul că penalizările pentru o decizie eronată sînt foarte mari Prima dintre aceste condiţii ţine de experienţa cercetătorului : în cercetarea operaţională continuă să aibă o pondere importantă arta şi îndemînarea cu care lucrezi A doua este un criteriu mai obiectiv Modelele de simulare sînt o formă a aproximării, dar aproximarea este de un tip cu totul diferit aceleia rezultate din procedurile clasice de construire a modelelor Utilizînd un mare volum de muncă, este posibil a simula procesul de aşteptare de o manieră foarte apropiată de realitate Fidelitatea modelului de simulare poate fi mai ridicat decît a oricărui alt model construit prin metode matematice pure Calitatea aproximativă a simulării este o consecinţă, aşa cum am văzut în capitolul precedent, a generării de numere aleatoare, adică a introducerii erorilor statistice de selecţie Dar, prin multiple încercări, acestea pot fi reduse şi în principiu simulările pot avea un grad de precizie oricît de înalt Din păcate, nu este întotdeauna posibil a aloca sume imense corespunzătoare costului unui număr foarte mare de cicluri de simulare şi de aceea este recomandabil să se utilizeze al doilea tip de aproximaţie Şi deoarece pentru a găsi o soluţie operaţională nu este întotdeauna absolut necesar un rafinament extrem, aproximarea de tipul al doilea este validă şi folositoare Principiul aproximării de tipul al doilea este următorul : un sistem este simulat în anumite condiţii, astfel încît se obţin puncte pe suprafeţele de răspuns ce interesează Acestea devin apoi obiectul analizelor statistice cu scopul determinării unor relaţii empirice dintre diferiţii factori (staţii de servire, traficul, priorităţile, timpii de servire etc ), ce pot fi apoi trataţi ca un model de ordinul doi al sistemului Este o înlocuire a experimentării sistemelor din lumea reală Un exemplu de utilizare al acestei tehnici este prezentat în cap 15 Al treilea tip de aproximare este cel clasic, utilizat în matematică, fizică şi în mai toate domeniile Acesta este exemplificat prin două metode : Dezvoltarea unor funcţii imposibil de manipulat în serii de puteri, în care termenii de rang mare pot fi ignoraţi, pentru anumite domenii de variaţie ale variabilelor Substituirea valorilor iniţiale în formulă, urmată de procesul iterativ Aceste proceduri clasice au fost utilizate în studiul modelelor de aşteptare, dar nu atît de mult cum s-ar putea crede Avantajul acestui tip de aproximare (comaparativ cu celelalte două) este relativa uşurinţă în estimarea erorilor rezultatelor Probabil că cea mai specta- culoasă utilizare a acestui tip de aproximare a fost realizată de Eiordan, care a ajuns la formulele, prezentate în capitolele 3 şi 4, pentru repartiţiile timpului de aşteptare în cazuî modelelor : (oo/SIEO) şi MjMIc: (oo/SIEO), formule care au constituit baza calculelor făcute pentru realizarea graficelor necesare determinării rapide a rezultatelor acestor modele Merită să mai menţionăm alte două exemple de utilizare al primului tip de aproximare matematică şi unul pentru al doilea tip; toate trei se referă la modelul de aşteptare M/D/c : (oo/FIFO) Molina a determinat următoarele formule aproximative pentru timpul mediu de aşteptare şi pentru probabilitatea ca toate staţiile de servire să fie ocupate: wq = P(t> o)— -1 [ic(l-p)c+l 1 - pc Xce z c P« > 0) = „ C~\Ce T c (8-1) 1 - e-“ V (X”/n!) H n = c C ! C X 1 Pce : cf2^c (1-p)2 1 Pce : cf2^c (1-p)2Pollaczek a determinat formule mai simple pentru aceleaşi mărimi ce sînt folositoare în special cînd există un număr mare de staţii de servire Eezultatele obţinute sînt : W9 = 1 o'e*1-^ (8 2) (8 2)P(t > 0) = 1 — W9(0) = 4= i — p p y 2kc Cititorul va remarca că utilizînd aproximarea factorialelor dată de Stirling este posibil a ajunge la rezultatul lui Pollaczek pentru Wa plecînd de la rezultatul lui Molina în sfîrşit, în 1953 Everett publică o procedură iterativă de calcul a probabilităţilor stării stabile, p„, pentru acest model de aşteptare Procedura sa, care este prea vastă pentru a o prezenta complet, depinde de observaţia că pentru valori mari ale lui n Pn + c Pn- Această relaţie a fost combinată cu ecuaţiile lui Crommelin pentru p,,, deducîndu-se o relaţie de determinare a lui ce pare a juca un iol în modelul M/D/c: (oo/GD) similar cu cel jucat de p în modelul MIM/c : (oo/GD) Ca urmare, Everett a introdus valoarea sa pentru 4» în formulele lui Erlang legate de modelul anterior, pentru a obţine noi estimaţii pentru pn Aceste valori sînt apoi utilizate în ecuaţiile modelului cu timpi constanţi de servire, procesul fiind iterat pînă la atingerea unor estimaţii stabile Apare astfel necesitatea ca cercetătorii să dezvolte aproximaţiile matematice de acest tip şi este surprinzător că pînă acum nu s-a acordat mai multă atenţie acestor probleme Cu aceasta am încheiat partea teoretică şi „filozofică-’ a lucrării Următoarele capitole, exceptînd concluziile, tratează exclusiv probleme ale aplicaţiilor practice SISTEME DE REZERVARE DE LOCURI Primul dintre sistemele reale de aşteptare pe care le vom analiza este cel observat pe parcursul studiului sistemului de transport de mărfuri al unei companii aeriene Fără îndoială că există sisteme similare şi în alte situaţii, iar problemele legate de acest sistem sînt, probabil, printre cele mai comune pe care le întîlneşte un analist de cercetare operaţională Sistemul în sine este destul de adecvat,pentru apliearea teoriei aşteptării, astfel încît este util să începem analiza cu o astfel de problemă 9 1 DESCRIEREA PROBLEMEI într-unul din oraşele principale pe care le serveşte, Compania aviatică Europa Airways Corporation (EAC) întreţine un vast departament pentru organizarea transportului mărfurilor Atunci cînd un expeditor doreşte să obţină spaţiu în cadrul unui zbor EAC, spre un alt oraş, îl poate obţine chemînd la telefon Centrul de rezervări EAC, existent în cadrul aeroportului (detaliile privind procedura de rezervare nu ne interesează) Toate apelurile provenite de la clienţi sînt direcţionate, prin intermediul unui tablou de distribuţie, spre un grup de zece linii telefonice, iar prin intermediul acestor linii apelul poate ajunge la cinci angajaţi Adică oricare dintre cei cinci angajaţi pot rezolva un apel sosit pe oricare dintre cele zece linii Dacă apelul soseşte atunci cînd tabloul de distribuţie este blocat, apelul va fi respins sau pierdut; dar acest caz are loc suficient de rar şi-1 vom putea ignora Dacă toate cele zece linii sînt ocupate, apelul este reţinut la tabloul de distribuţie; dacă există linii libere, dar toţi agenţii sînt ocupaţi, apelul este reţinut pe una dintre liniile libere în ambele cazuri apelul aşteaptă Departamentul de relaţii cu publicul din cadrul EAC a primit o serie de plîngeri din partea clienţilor că, într-o anumită perioadă a zilei, ei trebuie să aştepte excesiv de mult înainte ca apelurile lor să fie preluate de către agenţi în acelaşi timp în cadrul departamentului au apărut anumite probleme de utilizare a forţei de muncă, care implică necesitatea unei reorganizări radicale Ca urmare, directorul general a elaborat un proiect de reorganizare care, printre altele, afectează şi serviciul de rezervări telefonice în proiectul său se arată că activitatea serviciului de rezervare poate fi clasificată în funcţie de destinaţia geografică a coletelor şi, după o scurtă analiză, el a stabilit cinci categorii de trimiteri, anume spre : 1) Insulele britanice; 2) nordul şi estul Europei; 3) sudul Europei; 4) Orientul Apropiat; 5) Africa în proiectul său va exista cîte o secţie pentru fiecare dintre aceste categorii de trimiteri, iar în cadrul fiecărei secţii angajaţii vor efectua rezervări, înregistrări, întocmiri de facturi etc Centrul de rezervări actual va fi desfiinţat şi fiecare secţie va avea propriile linii telefonice în timpul discuţiilor purtate directorul general al Companiei EAC a relevat faptul că traficul telefonic legat de fiecare dintre cele cinci zone este diferit (de exemplu există de două ori mai multe apeluri pentru trimiteri de colete către nordul şi estul Europei, decît cele spre Orientul Apropiat) şi că existenţa unui singur telefon în fiecare secţie nu va asigura un serviciu egal tuturor clienţilor El susţine că probabil anumite seeţii vor necesita două sau chiar trei linii telefonice şi angajaţi Astfel, numărul total de angajaţi în cele cinci secţii poate fi mai mare decît cinci Se pune întrebarea : „Poate să-şi permită compania o asemenea creştere de personal” ? în acest moment putem considera că este cerută părerea unui analist de cercetare operaţională (pe care-1 vom numi dr Wiley) 9 2 PRIMA FAZĂ Primul pas făcut de dr Wiley a fost de a vizita Centrul de rezervări în scopul de a observa operaţiile curente El a constatat astfel că în faţa fiecărui angajat există un sistem de vizualizare în care fiecare dintre cele zece linii este reprezentată de un bec alb şi unul roşu Atunci cînd un apel aşteaptă pe una dintre linii, se aprinde becul roşu corespunzător; iar cînd apelul a fost preluat, se aprinde becul alb în consecinţă, este destul de uşor să culegem date privind timpii totali de aşteptare la telefon; timpul total de aşteptare al unui apel se poate cronometra între momentul aprinderii becului roşu şi momentul aprinderii becului alb corespunzător La apeluri nu se răspunde în ordine, căci angajaţilor le este greu să reţină ordinea de aşteptare Apelurile sînt rezolvate, aparent, în ordine aleatoare Dr Wiley a decis că trebuie să aibă anumite date înainte de a continua analiza; din păcate însă liniile telefonice nu au dispozitive de măsurare; cu toate acestea, el găseşte posibilitatea de a aranja cu şeful centralei telefonice a aeroportului ţinerea evidenţei tuturor apelurilor către centrul de rezervări în timp de o săptămînă La sfîr- Fig 9 1 Apeluri Înregistrate la centrul de rezervare (8 martie 1957) şitul săptămînii primeşte datele şi constată că numărul apelurilor în zilele de luni pînă vineri este repartizat în acelaşi mod (v fig 9 1) Sîmbătă şi duminică apelurile sînt mai puţin numeroase iar repartiţia lor este alta Dr Wiley şi-a dat seama că nu sfîrşitul de săptămînă trebuie avut în vedere, deoarece atunci sînt probleme mai puţin numeroase Trebuie acordată atenţie doar celor cinci zile lucrătoare; în special trebuie concentrată atenţia asupra comportamentului sistemului în perioadele de aglomeraţie, pentru că în aceste perioade este cel mai probabil ca sistemul să fie supraîncărcat Diagrama arată că perioadele ocupate din zi sînt între 10,30—12 şi 14,30—15,30 Următorul pas este de a observa sistemul în timpul perioadelor de aglomeraţie Dr Wiley şi asistentul său reuşesc să înregistreze următoarele informaţii legate de apelurile din perioadele de aglomeraţie, folosind sistemul de vizualizare al angajaţilor : intervalele dintre apelurile consecutive; timpii de aşteptare pentru mai multe apeluri; timpii de servire pentru mai multe apeluri Intervalele dintre apelurile consecutive, pe parcursul a patru intervale de jumătate de oră din perioadele ocupate, sînt testate prin testul indicelui de dispersie Poisson bazat pe gruparea în 15 intervale de cîte 2 minute Toate cele patru eşantioane ale apelurilor dau dreptul să se facă ipoteza că sosirile sînt repartizate Poisson Există două tipuri de apeluri telefonice : 10 % sînt informative şi sînt rezolvate foarte rapid (de exempiu : „Mîine există vreo cursă spre Praga?”), celelalte sînt apeluri pentru rezervări de spaţii de transport Repartiţiile timpilor de servire pentru ambele tipuri de apeluri sînt unimodale în apropierea lui zero şi, după cum a putut remarca cu uşurinţă Wiley, pot fi suficient de bine aproximate cu repartiţii exponenţiale, cu mediile de 3,2 şi respectiv 3,7 minute Dr Wiley a remarcat de asemenea un fapt semnificativ în legătură cu activitatea depusă de un angajat atunci cînd răspunde la un apel de rezervare în aproape toate cazurile, agentul trebuie să consulte unul din fişierele, dealtfel excelent proiectate, în care să verifice dacă există spaţiu disponibil la zborurile pentru care este interesat clientul său Un desen al poziţiei acestor fişiere şi al angajaţilor este prezentat în fig 9 2 Cititorul va remarca, aşa cum a făcut-o şi dr Wiley, eă în mod frecvent agenţii trebuie să-şi părăsească locurile pentru a consulta fişierele şi aceasta este o componentă importantă a timpului de servire în acest moment analistul nostru se retrage în biroul său pentru a analiza cele observate Este clar că sistemul existent poate fi clasi- ficat ca fiind de tipul Jf/Jf/5 : (oo/SIRO); deci un model care nu creează probleme deosebite şi care a fost foarte mult studiat Care este sistemul propus ? în loc de un grup de zece linii cu cinci angajaţi, disponibili oricînd, vor fi cinci secţii, fiecare de tipul Jf/Jf/c : (Z/?); din punctul de vedere al analizei, un rol important îl joacă c, L şi 12 3 4 5 □ □ □ □ □ Angajaţi 12 3 4 5 □ □ □ □ □ Angajaţi Fişiere Fişiere Fig 9 2 Schiţa oficiului de rezervare regula încă nedeterminată a disciplinei de aşteptare Nu se poate încă stabili nimic precis Dr Wiley mai face cîteva verificări asupra sistemului existent, utilizînd datele culese Tabelul 9 1 prezintă cîteva date legate de timpii de aşteptare, înregistrate pe parcursul celor patru intervale de jumătate de oră în timpul perioadelor de aglomeraţie, atunci cînd frecvenţa apelurilor este de aproximativ unul pe minut (Nu toate datele culese pot fi utilizate, deoarece frecvenţele medii de sosire sînt prea diferite între ele pentru a putea fi comparate ) Prima coloană a tabelului prezintă rezultatele unui model teoretic de tipul Jf/Jf/5 : (oo/SIRO) Corespondenţa este suficient de corectă şi ne putem imagina satisfacţia dr-lui Wiley în faţa obţinerii unei aproximări satisfăcătoare a realităţii Tabelul 9 1 Timpi observaţi şi teoretici de aşteptare Af/Af/5; ţoc/SIRO) 1 Intervalul 1 Intervalul 2 Intervalul 3 Intervalul 4 P(v> 5 min) j 1,30 0,17 1,10 ! 1,80 rezultate comasate 0,13 1,43 1,35 Analistul nostru consideră că acum poate să se gîndească la noua propunere Dr Wiley începe prin a culege mai multe date (Atît conducătorii cît şi personalul de execuţie au întotdeauna o mare nedumerire asupra volumului mare de informaţii numerice, pe care îl necesită cercetarea operaţională ) Făcînd observaţii asupra liniilor telefonice, pe baza unui plan de culegere eşantionată a acestora, pe o perioadă de două săptămîni, el obţine următoarele rezultate privind apelurile telefonice relative la zona geografică a destinaţiei transportului solicitat : Insulele britanice Nordul şi estul Europei Sudul Europei Orientul Apropiat 17,0% 21,6% 28,6% 12,2% 5 Africa 26,6% Zona 3 primeşte mai mult decît dublul comenzilor pentru zona 4 Dr Wiley verifică apoi dacă acest aspect al repartizării procentuale a apelurilor pe zone de destinaţie este în concordanţă cu transporturile efectuate şi găseşte că este Mai mult, el constată că repartizarea cererilor de transport spre diferite zone de destinaţie este aproximativ constantă în timpul perioadei de vîrf a verii şi de asemenea în timpul perioadei de vîrf premergătoare anului nou Aceasta îi întăreşte convingerea că procentele determinate pentru numărul de apeluri aferente diferitelor zone de destinaţie sînt aplicabile şi perioadelor de vîrf ale anului (Nu e nici o îndoială că această procedură are la bază presupuneri foarte „tari”, dar ce altceva poată să facă?) Dr Wiley mai face încă o verificare, căci este posibil ca procentele determinate să fie corecte în sensul mediei, dar numărul de apeluri pentru diferitele zone de aşteptare pot să atingă maximul la ore diferite în timpul zilei, astfel încît, la un moment dat al zilei, pot fi mai multe apeluri pentru zona 4, decît pentru zona 3 Aceasta nu ar fi avut importanţă în vechiul sistem, dar va fi foarte important dacă va fi adoptat noul sistem Analiza înregistrărilor făcute de el însuşi şi de asistentul său îi arată că deşi sînt oarecari neconcordanţe a momentelor de vîrf, rapartizarea procentuală obţinută pentru apelurile către diferite zone, relativă la perioadele de vîrf ale primelor cinci zile ale săptămînii, este în esenţă corectă Dr Wiley va putea face o serie de calcule numai dacă rezolvă problema disciplinei de aşteptare ce va exista în noul sistem Directorul companiei ce răspunde de telecomunicaţii îi spune să se prevadă înlocuirea tabloului central de distribuţie cu unul automat Un echipament automat va funcţiona pentru fiecare linie, şi nici un apel nu va mai putea fi menţinut în tabloul de distribuţie Deci regula FIFO va determina disciplina de aşteptare Dacă va fi menţinut centrul de rezervări, adică sistemul actual, el va dispune de o capacitate de 20, avînd cel mult zece linii telefonice Dacă va fi adoptat noul sistem, atunci fiecare secţie va dispune de cinci Dr Wiley ştie că modelul cel mai apropiat pentru noul sistem ce a fost propus este, pentru fiecare din cele cinci zone, M]M/c : (5/FIFO) La un nivel al traficului de 3,7 erlangi, intrările în modelul Jf/Jf/5 : (20/FIFO) reprezentînd sistemul existent sînt de 3,7 erlangi, iar intrările pentru fiecare zonă propusă sînt respectiv 0,64; 0,80; 1,06; 0,45 şi 0,75 erlangi, corespunzînd procentelor de apeluri ce au fost stabilite Se pot face comparaţii la acest nivel al traficului, deoarece acesta este o reprezentare rezonabilă a încărcării ce trebuie rezolvată; putem compara timpii medii de aşteptare, probabilităţile de pierdere şi probabilităţile ca un timp de aşteptare să depăşească un minut Scopul constă în a calcua aceste cantităţi pentru sistemul actual cu cinci staţii de servire şi apoi de a determina necesarul de staţii de servire pentru diferitele zone ale sistemului propus, capabil să asigure aceeaşi calitate a serviciului Pentru modelul : (10/SIRO) nu există o formulă convena bilă, astfel că Wiley va trebui să obţină rezultate prin simulare Pentru modelul M/M/S : (20/FIFO) ca şi pentru JH/M/c : (5/FIFO) există formula (4 6) Cu aceste mijloace el a obţinut rezultatele prezentate în tabelul 9 2 în raportul prezentat, dr Wiley scrie : „Reorganizarea propusă pentru sistemul de rezervare va conduce la un serviciu mult mai nesatisfăcător pentru clienţi Pentru a egala actualul standard de servire, care este suficient de bun, vor fi necesari cel puţin doi salariaţi în fiecare zonă, în timpul orelor de aglomeraţie Aceasta va însemna o creştere de 100% a personalului Odată cu creşterea traficului, sistemul propus devine şi mai nesatisfăcător decît cel actual, căci creşterea numărului de angajaţi tinde să se facă cu multipli de 5 (numărul zonelor) Sistemul actual, cu toate telefoanele concentrate şi capabile să preia orice apel, admite o dezvoltare mai graduată odată cu creşterea traficului Trabuie renunţat la ideea separării pe zone geografice şi trebuie reţinut din sistemul actual modul de Tabelul 9 2 Comparaţie între sistemul existent şl cel nou Sistemul existent Sistemul propus Numărul de angajaţi 5 11 Procentul de apeluri ce găsesc toate liniile ocupate 1,3 1,6 Timpul mediu de aşteptare (min) 1,1 1,2 Procentul de apeluri ce aşteaptă mai mult de 1 min 31 39 grupare a telefoanelor Dar pentru a face un sistem de rezervare mai eficient, este de dorit să se revadă modul de obţinere a informaţiilor de către angajaţi şi de asemenea modul de păstrare a informaţiilor” Raportul a fost prezentat Argumentele sînt bine primite : cu toate că trebuie făcut ceva în organizarea sistemului de rezervare, centrul de rezervare este menţinut la fel cum a fost; ei bine, să zicem că aproape cum a fost Dr Wiley, cu sentimentul că a făcut o treabă bună, pleacă de la şedinţă, de la Companie şi din această carte INTERLUDIU Timpul trece şi au loc schimbări în sistemul de transport Anterior departamentul era împărţit din punct de vedere funcţional în sectoarele : contabilitate, preţuri, documentaţie şi vamă-expediţie Anumiţi angajaţi aveau sarcina de a ţine contabilitatea încasărilor, precum şi de a înregistra telefoanele prin care se cereau informaţii sau se anula vreo comandă Această împărţire nu mai este acum valabilă Departamentul de transport este împărţit acum în secţia de rezervări elefonice şi alte sectoare care se ocupă de preţuri, documentaţie, vamă^expediţie şi urmărirea transporturilor cu o anumită destinaţie Pe lingă acestea, telefoanele prin care se solicită informaţii sînt preluate de către sectorul corespunzător zonei geografice respective Sistemul telefonic a fost de asemenea îmbunătăţit Anterior, orice angajat putea să răspundă la orice apel telefonic ce solicita informaţii de orice fel (adică era un sistem de informare complet) Acum cele cinci perechi de circuite telefonice interne de la centrala aeroportului sînt repartizate la cele cinci sectoare Clienţii permanenţi au fost sfătuiţi să ceară o anumită linie atunci cînd îşi rezervă locuri pentru o anumită destinaţie De exemplu, dacă se solicită o rezervare pentru Geneva şi dacă una sau ambele linii care deservesc sectorul 3 nu sînt ocupate, se conectează direct pe legăturile 31 sau 32 Dacă ambele sînt ocupate, apelul telefonic este conectat la orice altă linie liberă Acest sistem s-a dovedit nesatisfăcător pentru departament, deoareee clienţii s-au plîns că secţia de rezervări telefonice funcţionează defectuos Atunci s-a încercat să se amelioreze situaţia, rein- troducîndu-se angajaţi în fiecare sector care să se ocupe de rezervările telefonice Directorul îl însărcinează pe dl Staple să întreprindă o nouă anchetă FAZA A DOUA Din fericire Staple este un om cu vederi largi care verifică totul în practică Cele zece linii telefonice din secţia de rezervări vor fi şi pe viitor înpărţite în cinci grupuri a cîte 2 linii, dar fiecare grup va recepţiona numai acele apeluri telefonice care sînt strict legate de activitatea sectorului respectiv Telefoniştii de la centrală vor reţine acele apeluri care nu pot fi conectate la nici una din cele două linii corespunzătoare în timpul orelor de vîrf, o colaboratoare va controla la întîmplare diferite linii pentru a vedea cum decurge activitatea Telefoniştii de la centrală (care, deşi comentează, sînt mulţumiţi că participă la un proiect special) au fost de acord să numere apelurile telefonice dintr-o oră pentru fiecare sector Dl Staple a cronometrat şi el la întîmplare aceste apeluri Experienţa a durat două săptămîni, iar imensul volum de date strînse a fost analizat Iată cîteva dintre rezultatele obţinute Din tabelul 9 3 rezultă că timpul mediu de servire al fiecărui apel telefonic este mai scăzut faţă de cel obţinut cu 12 luni înainte, dar diferenţa nu este semnificativă, la nivelul de 5% Tabelul 9 3 Comparaţie intre timpul de funcţionare şi cel de aşteptare Studiul I Studiul II Timpul mediu de servire al unui apel telefonic 3,7 3,7 Timpul mediu de aşteptare al unui apel telefonic 1,35 1,58 însă nu s-a mai răspuns la apelurile scurte care solicitau informaţii Dl Staple a tras concluzia în mod just că angajaţii, fiind specializaţi pentru o zonă geografică relativ restrînsă, pot să prelucreze mai repede apelurile telefonice care solicită rezervări Timpul mediu de aşteptare în orele de vîrf este mai mare decît era anterior, fapt datorat : unui număr sporit de apeluri în timpul orelor de vîrf; unei discrepanţe mai accentuate între procentele apelurilor corespunzătoare diverselor sectoare; introducerii grupurilor cu acces limitat Este adevărat că numărul mediu al apelurilor telefonice care sosesc la centrală pe oră în timpul perioadelor de vîrf a crescut, atingînd circa 75 apeluri pe oră Nu există nici un motiv să presupunem că repartiţia lor nu este Poisson Procentele apelurilor pentru diverse sectoare înregistrate în timpul experienţei au fost următoarele : 1) Insulele britanice 17% (17,0) 2) Europa de nord şi est 22% (21,6) 3) Europa de sud 29% (28,6) 4) Orientul Apropiat 12% (12,2) 5) Africa 20% (20,6) Numerele din paranteze reprezintă rezultatele analizei dr Wiley Aparent nu s-a produs nici o schimbare şi putem să tragem concluzia, la fel ea dl Staple, că principala cauză a deteriorării serviciului o reprezintă faptul că nu există un acces complet Adepţii teoriei aşteptării sînt satisfăcuţi, dar dl Staple este impresionat de reducerea timpului de servire datorită specializării Mai mult, el a mers mai departe decît dr Wiley (care observase că angajaţii trebuiau să-şi irosească o parte din timp cu căutarea informaţiilor adecvate şi sugerase să se facă ceva în această direcţie) şi a efectuat măsurări care au arătat că pentru acele convorbiri telefonice care nu-i obligă pe angajaţi să îşi părăsească locul, timpul mediu de servire pentru un apel este de 2,3 minute faţă de 3,5 minute în general Dl Staple consideră că dacă rezervările pentru diferite zboruri s-ar face telefonic de către angajaţii împărţiţi pe sectoare, atunci informaţiile privind aceste sectoare pot fi păstrate în stative special amplasate între cele două posturi telefonice care pot recepţiona apelul telefonic Ambii angajaţi pot ajunge cu uşurinţă la ele fără să-şi părăsească locul El face o nouă paralelă, avînd în vedere trei sisteme posibile : aceiaşi angajaţi se ocupă de rezervări şi de informaţii privind toate zonele geografice, avînd acces la toate liniile telefonice; un grup de angajaţi se ocupă de rezervări pentru toate zonele geografice, informaţiile însă sînt date de angajaţi speciali în fiecare sector în parte; rezervările şi informaţiile pentru fiecare sector sînt date separat de către angajaţii sectorului corespunzător Timpul mediu de servire care este luat în considerare în calcule este dat în tabelul 9 4 Tabelul f> 4 Timpii de servire Sistemul 1 (a) (b) (c) Apeluri telefonice pentru rezervări 3,5 3,5 2,3 Apeluri telefonice pentru informaţii 1 3,2 1 2,8 2,8 Aceste date se bazează pe observaţii directe în cazul sistemelor (a) şi (b) şi pe o ipoteză suplimentară în cazul sistemului (c) Această ipoteză constă în faptul că buna cunoaştere : informaţiilor de către angajaţii care se ocupă de rezervări într-un anumit sector reduce durata medie a convorbirilor pentru rezervări la cea a apelurilor telefonice pentru rezervări de care se ocupă un angajat în sistemul (a) atunci cînd acesta nu trebuie să-şi părăsească biroul (adică nu pierde timp să caute informaţiile cerute în altă parte) Sistemul (a) are meritul că fiecare angajat poate răspunde la apelurile pentru orice zonă, dar presupune o risipă mai mare de timp din partea angajaţilor care trebuie să-şi părăsească locurile pentru a consulta informaţiile Sistemul (c) impune funcţionarilor să-şi părăsească locurile mai puţin dar impune o diferenţiere mai mare a tipurilor de apeluri şi întîrzieri datorite transferului apelurilor dintr-un sector în altul Acest fapt este important deoarece 21% dintre apelurile pentru rezervări includ tranzacţii multiple (doar 10 % în 1957) Sistemul (b) este un compromis Pentru fiecare sistem a fost calculat (tabelul 9 5) numărul personalului necesar în timpul unei ore de vîrf pentru a asigura servirea clienţilor (măsurat cu timpul mediu de aşteptare de circa 1 min/apel) Sistemul (b) este mai puţin recomandabil şi alegerea se va face între sistemul (a) şi sistemul (c) Concluziile doctorului Wiley au fost răsturnate deoarece este evident că cea mai bună soluţie este ca angajaţii care se ocupă de rezervări telefonice să fie amplasaţi în sectoarele respective; De fapt dl Staple a folosit acelaşi model ca şi dr Wiley, modelul clasic Tabelul 9 5 Necesar de personal în orele de vîrf Sistemul Sistemul (a) | (b) (a) | (b) (c) (c) Angajaţi care lucrează permanent (care nu se ocupă clecît de convorbirile telefonice) Angajaţi care lucrează un timp limitat Angajaţi care lucrează permanent (care nu se ocupă clecît de convorbirile telefonice) Angajaţi care lucrează un timp limitat Folosirea angajaţilor per- manenţi Folosirea angajaţilor pe timp limitat Folosirea angajaţilor per- manenţi Folosirea angajaţilor pe timp limitat Sectoarele 2, 3 : 2 Sectoarele 1, 4, 5 : 1 Total = 7 70% 40% Sectoarele 2, 3 : 2 Sectoarele 1, 4, 5 : 1 Total = 7 70% 40% i Sectoarele 2 3 : 2 Sectoarele } 1 4, 5: 2 Total = 12 i Sectoarele 2 3 : 2 Sectoarele } 1 4, 5: 2 Total = 12 38% 38% M/MIc: (oo/FIFO) Deosebirea a apărut deoarece dl Staple şi-a dat seama că sistemul poate fi ameliorat inventînd un mijloc de reducere a timpului de servire, lucru evident dar, de cele mai multe ori, neglijat! COMENTARIU Atît Wiley cît şi Staple au dreptate şi se înşeală Wiley a preferat (în mod just) metoda servirii pe mai multe canale specializării înguste, iar Staple a insistat corect asupra scurtării timpului de servire Cititorul va lua în consideraţie ambele concluzii Nu este foarte greu să se inventeze un sistem de vizualizare, sau bazat pe computere, care să permită funcţionarilor de la rezervări, care, se află într-o sală, să răspundă la toate întrebările şi să facă toate operaţiile de rezervare repede şi precis fără să-şi părăsească locurile Această soluţie a fost adoptată pe scară largă acum, după mulţi ani, de către multe companii aeriene în limbajul capitolului 7, configuraţia optimă a timpului de aşteptare este influenţată de anumiţi factori — în acest caz, de timpul mediu de servire al apelurilor telefonice Pentru o anumită valoare, soluţia lui Wiley a fost cea mai bună; pentru alta cea a lui Staple Staple a făcut aproximări cam grosolane folosind tabelele lui Molina pentru Jf/3f/c : (oo/FIFO) pentru a deduce valorile aproximative din tabelul 9 5 Ar fi fost mai bine dacă ar fi fost folosite modelele M/M/c : (20/FIFO) şi Jf/Jf/e : (5/FIFO) ca mai bune ilustrări ale comportării viitoare a sistemului Deşi rezultatele estimative ale oamenilor diferă, totuşi clasificarea celor trei sisteme posibile este aceeaşi Această concluzie a mai fost trasă în cap 8 10 METODE DE CONTROL AL CĂLĂTORILOR Al doilea caz pe care-1 vom examina este sistemul de verificare a pasagerilor la aeroport Majoritatea cititorilor cunosc sistemul de verificare, iar pentru cei care nu sînt familiarizaţi vom da o descriere a acestei operaţii Atunci cînd un pasager soseşte la aeroport să ia avionul, mai întîi trebuie să se prezinte la un ghişeu de verificare unde angajatul liniei aeriene, numit şi agent de pasageri, cîntăreşte bagajele călătorului, înregistrează rezultatul cîntăririi pe talonul de zbor al călătorului (acea parte a biletului care indică despre care zbor este vorba), etichetează bagajele călătorului, le trimite la camera de bagaje, îi dă călătorului contramarca bagajelor şi, în unele cazuri, îi ia cuponul de zbor şi îi dă în schimb un permis de trecere După ce au fost îndeplinite toate aceste operaţii, călătorul este admis pentru zbor Cu excepţia aeroporturilor mici, în anumite perioade ale zilei verificarea călătorilor pentru mai multe zboruri se face în acelaşi timp şi compania aeriană trebuie să asigure mai multe ghişee care să servească pasagerii Există numeroase moduri de a organiza aceste ghişee pentru a asigura o servire eficientă, însă majoritatea celor folosite în mod curent sînt amestecuri judicioase (mai mult sau mai puţin) a două tipuri extreme Primul tip poate fi denumit sistem cu restricţie completă sau sistem de zbor în această formă de organizare fiecărui zbor îi este distribuit un grup de ghişee (adesea unul singur) Toţi pasagerii pentru un anumit zbor trebuie să fie verificaţi la ghişeele special desemnate, care nu se ocupă decît de pasagerii acestui zbor Cel de-al doilea tip este cunoscut fie ca sistem liber, sau comun, sau ca sistem cu acces liber în acest sistem fiecare pasager care soseşte poate să fie verificat la orice ghişeu al companiei aeriene Primul sistem are dezavantajul că pasagerii aşteaptă de regulă la un ghişeu în timp ce altul, repartizat pentru un alt zbor, este liber; se poate admite ca acest fenomen nu apare în cazul unui sistem cu acces liber Lipsa principală a sistemului cu acces liber este că funcţionarul de la ghişeu nu ştie pentru ce zboruri stau la coadă pasagerii în faţa ghişeului său Călătorul care vine tîrziu pentru un anume zbor poate sta la coadă după alţii care sînt verificaţi pentru zboruri ulterioare şi astfel poate întîrzia pînă la a pierde avionul; funcţionarul nu poate să anticipeze asemenea situaţii în anii 1950, în Anglia, liniile aeriene foloseau exclusiv metoda tradiţională de verificare şi nu aveau experienţă în ceea ce priveşte caracteristicile operaţionale ale sistemelor obişnuite în 1955 s-a prevăzut că creşterea rapidă a traficului aerian şi introducerea iminentă a avioanelor mai mari face ca sistemul de verificare să fie neeconomic şi incomod, cel puţin pentru aeroporturile mari şi aglomerate 10 1 PRIMA FAZĂ în 1956 secţia de cercetări operaţionale a companiei BEA a întreprins un studiu asupra sistemului de verificare a călătorilor la aeroportul din Londra Acest studiu avea un scop dublu : să stabilească în mod obiectiv cum funcţionează sistemul de verificare existent şi să-i stabilească perspectivele viitoare Sistemul de verificare existent era un sistem de tip mixt Toţi pasagerii care soseau să fie verificaţi cu mai mult de o oră înainte de plecare erau serviţi de către unul dintre cele două ghişee speciale de înregistrare Pasagerii care soseau mai tîrziu erau verificaţi la celelalte ghişee înregistrările amănunţite ale momentelor sosirilor pasagerilor au arătat că aproximativ 30% din totalul călătorilor au fost verificaţi la ghişeele speciale de înregistrare Figura 10 1 prezintă graficul sosirii călătorilor exprimat prin proporţia F(z, t) a călătorilor sosiţi în diferite momente pînă la ora închiderii (adică pînă la cel mai tîrziu moment posibil înainte de plecare) pentru verificare După închidere nu se mai verifică nici un călător (zborul a fost „închis”) Din fig 10 1 reiese că 70% din călători au fost verificaţi la alte ghişee De fapt călătorul care soseşte mai tîrziu decît cu o oră înainte de ora zborului său, dar cu 25 min înainte de zbor, poate să fie verificat la un ghişeu special purtînd numărul zborului său (Cititorul va observa că compania aeriană, în dorinţa de a-şi satisface clienţii, I-'ig 10 1 Modul de sosire a pasagerilor verifică şi călătorii care sosesc mai tîrziu de ora închiderii ghişeului, adică cu 25 min înainte de zbor Cu toate acestea, principiul rămîne în vigoare ) Figura 10 2 indică faptul că clienţii sosesc de obicei în grup organizat Ca urmare, ca unitate de măsură a intrărilor s-a luat grupul de călători S-a examinat modul de intrare în sistemul de aşteptare la ghişeele speciale de înregistrare Numerele grupurilor de călători care au sosit la interval de cinci minute pe parcursul celor două ore de vîrf ale celor cinci zile luate sub observaţie sînt prezentate în tabelul 10 1 Din experienţă, este de aşteptat ca intrările să nu fie repartizate Poisson; în fond, au fost observaţi călători provenind de la im număr fix de surse (zboruri) care nu au toate aceleaşi efecte Cititorul trebuie să-şi reprezinte în mod clar repartiţia totală obţinută prin suprapunerea mai multor modele de tipul indicat în fig 10 1, cu momentele de „închidere” a zborurilor decalate cîteodată cu 5 şi alteori chiar Fig 10 2 Repartiţia numărului pasagerilor pe grup cu 10 minute unul de altul Tabelul 10 1 indică faptul că totuşi ipoteza că intrările sînt repartizate Poisson nu este prea absurdă Evident, n-am avut un motiv pentru a face un test statistic de semnificaţie a acestor date, căci ştiam că sosirile nu pot fi în Tabelul 10 1 Repartiţia numărului grupurilor de pasageri eare au sosit la interval de 5 mln la ghişeele speciale de înregistrare (ore de aglomeraţie); îîr grupurilor pe interval Frecventa observată (%> Repartiţia Poisson teoretică Nr grupurilor pe interval! Frecventa observată (%) Repartiţia Poisson teoretică 0 7 3,3 6 7 7,2 1 12 11,3 7 4 3,5 2 20 19,3 8 3 1,5 3 24 21,9 9 2 0,6 4 12 18,6 10 1 0,2 5 7 12,6 11 1 (0,05) Media 3,4 3,4 mod real repartizate Poisson Altfel spus, repartiţia intrărilor, deşi nu este apriori Poisson, poate fi totuşi aproximată practic printr-o repartiţie Poisson A fost observată repartiţia timpului de servire la ghişeele speciale de înregistrare în timpul orelor de aglomeraţie (fig 10 3) Nu poate prelucrării grupurilor de pasageri fi repartiţie exponenţială dar, pe de altă parte, repartiţia timpului de servire chiar şi în teoriile clasice privind traficul telefonic nu se asemăna prea mult cu repartiţia exponenţială Grupul de cercetători a admis că repartiţia timpului de servire este exponenţială Trebuie să se sublinieze faptul că repartiţia timpului de servire pentru ghişeele obişnuite de verificare a fost identică cu cea pentru ghişeele speciale de înregistrare Aceasta nu este surprinzător, deoarece activitatea cerută pentru grupuri de pasageri este identică pentru ambele tipuri de ghişee Figura 10 3 este compusă din suma datelor obţinute de la toate ghişeele în acelaşi timp cu înregistrarea timpului de servire şi a orelor de sosire a călătorilor, cercetătorii au înregistrat timpul de aşteptare al călătorilor la cele două ghişee speciale de înregistrare Aceste observaţii au fost verificate şi în alte perioade, nu neapărat de aglomeraţie Apoi, timpii de aşteptare înregistraţi, în perioade de 1/2 oră, au fost reprezentaţi grafic în raport cu utilizarea funcţionarilor în aceleaşi perioade Figura 10 4 prezintă aceste rezultate (cititorul îşi va aminti că u = 100 p, v cap 4) Sînt sigur că şi cititorii consideră aceste rezultate ca fiind încurajatoare Observaţiile empirice sînt în concordanţă cu previziunile modelului teoretic Jf/JI/2 :(co/FIFO) şi cercetătorii au considerat, în mod just, că în analiza ulterioară toate ghişeele trebuie tratate ca şi cum ar fi un sistem cu mai multe staţii de servire în paralel, sistem ce poate fi studiat cu ajutorul modelelor de tipul M/JH/c : (co/FIFO) Această extrapolare îndrăzneaţă a experienţei s-a dovedit a fi greşită Sistemul de aşteptare la ghişeele obişnuite la care au fost serviţi călători în proporţie de 71% este foarte neobişnuit Cercetătorii nu au considerat rentabil să-l analizeze folosind modele matematice Observaţiile empirice şi experienţa dobîndită de ghişeele obişnuite au stabilit cîteva lucruri importante Mai întîi, datorită ritmului variabil al sosirilor pasagerilor în timpul perioadei de verificare (de la 60 la 25 min înainte de momentul plecării), utilizarea unui singur funcţionar la un asemenea ghişeu a fost foarte variabilă (tabelul 10 2) Tabelul 10 2 Utilizarea funcţionarilor la ghişeele de verificare Perioada %) min înainte de plecare 60 — 50 53 50-40 1 53 40-31 i '6 30-25 ' 15 în timpul perioadei de 40—31 minute înainte de plecarea cursei, de obicei este necesar ca un al doilea funcţionar să vină să ajute; cei mai experimentaţi funcţionari chemau automat în ajutor un asemenea auxiliar cu 45 min înainte de plecare pînă cu 25 min înainte de plecare (în timpul acestor ultime 5 min, 30—25 înainte de plecare, deşi erau puţine sosiri ale călătorilor, cele care se înregistrau erau foarte critice; în special în timpul ultimelor sosiri de după ora de închidere, al doilea funcţionar nu era deloc în plus) Cu ajutorul acestor stratageme, funcţionarii de la ghişeele obişnuite de verificare au reuşit să menţină timpul de aşteptare al călătorilor la un nivel tolerabil, comparabil cu cel obţinut la ghişeele speciale i-'ig 10 4 Creşterea timpului de aşteptare a pasagerilor la ghişeele speciale de Înregistrare odată cu creşterea ocupării angajaţilor * Observaţii efective se bazează pe medii calculate pe perioade de 1/2 ori - Curba datA de M/M/2; (co/FIFO) de înregistrare Timpul mediu de aşteptare în orele de vîrf la ambele tipuri de ghişee varia între 1,2 şi 2,5 minute S-a pus întrebarea : este această valoare satisfăcătoare ? Conducerea companiei a considerat că nu şi s-a recomandat ca standard de calitate ca timpul mediu de aşteptare în orele de aglomeraţie să fie de cel mult un minut O altă problemă care se pune : este posibil ca acest sistem să fie extins pentru a face faţă unui volum sporit de călători? Se părea că ghişeele speciale de înregistrare nu pun nici o problemă Extinderea la trei sau patru va asigura acoperirea solicitărilor pentru următorii cinci ani Se estima dublarea numărului călătorilor între anii 1956 şi 1961 Unul din principalii factori care determină această creştere îl constituie introducerea avioanelor de capacitate sporită Desigur, se estima că va creşte şi numărul zborurilor, dar nu cu mult, însă creşterea numărului de pasageri trebuia să aibă loc în special prin înlocuirea vechilor tipuri de avioane cu tipuri noi, de capacitate mult mai mare decît cele vechi De exemplu, avioanele de tipul Vickers Van- guard au o capacitate dublă faţă de tipul mai vechi Viscount 701S Acest lucru va influenţa considerabil situaţia la ghişeele obişnuite Astfel, conform calculelor, cu 40—25 min înainte de plecarea avioanelor mari, vor fi necesari pînă la trei funcţionari pe zbor, iar cu 60—40 min înainte de zbor cîte doi, pentru a se evita aglomeraţia, lucru care a fost considerat inadmisibil Una dintre soluţiile propuse a fost de a scurta timpul de transfer de la ghişeele speciale, de la 60 min înainte de plecare la 40 min înainte de plecare, astfel incit ghişeele obişnuite să fie afectate doar cîte 15 min unui zbor Examinîndu-se această soluţie, s-a constatat că se poate economisi timp dirijîndu-se cît mai mulţi călători spre ghişeele speciale de înregistrare în acest scop s-au folosit ca modele Jf/Jf/5 : (oo/FIFO) pentru ghişeele speciale şi Jf/Jf/2 : (oo/FIFO) pentru perechile de ghişee obişnuite Ultimul model este evident foarte grosolan înainte ca acest sistem să poată fi pus în practică trebuiau să se rezolve cîteva probleme de funcţionare Deoarece sistemul de verificare existent nu presupunea mai mult de doi angajaţi — sau ghişee — pentru fiecare zbor, nu a existat o problemă a controlului Angajatul (angajaţii) putea(u) să verifice pasagerii urmărind lista, obţinîndu-se un control strict al numărului şi identităţii acestora Este clar că în cazul noului sistem apare necesitatea confruntării biletului intr-un punct central pentru toţi pasagerii unui zbor, deoarece un număr mai mare de angajaţi nu poate utiliza aceeaşi listă, iar identitatea moment cu moment a mai multor liste este greu de realizat Transferul biletelor de zbor intr-un punct central de control pe o bandă transportoare sau prin poştă pneumatică generează un alt şir de aşteptare la camera de control, unde biletele ar mai trebui să fie sortate în funcţie de destinaţia zborului Centrul de control ar trebui să îndeplinească deci următoarele funcţii: să sorteze biletele după indicativul zborului pe măsură se acestea sosesc pe bandă de la ghişeele obişnuite; să ia fiecare set ie bilete la închiderea zborului respectiv, să adune greutatea bagajelor şi să stabilească numărul pasagerilor, să comunice aceste date secţiei de control al încărcăturii, să verifice biletele după lista pasagerilor şi să observe neconcordanţele (pasageri care nu s-au prezentat, alţi pasageri suplimentari etc ); cînd este necesar, să expedieze datele care prezintă interes la aeroportul de destinaţie Unele dintre aceste operaţii erau efectuate de secţia de control, iar în sistemele de verificare existente nu erau necesară sortarea (biletele erau deja clasificate), adunarea greutăţii bagajelor era făcută de către funcţionarii delaghişee, iar verificarea biletelor pe baza listei pasagerilor la fel Trebuia să se evalueze mîna de lucru necesară pentru îndeplinirea tuturor acestor operaţii în noul centru de control Astfel, a fost stabilit timpul mediu de servire care va rezulta pentru ghişeele obişnuite dacă o parte din aceste operaţii vor fi transferate spre centrul de control S-a prevăzut că va fi în medie de 1,6 minute pentru un grup de călători şi s-a convenit ca repartiţia sa să fie considerată exponenţială De asemenea s-a calculat numărul de angajaţi necesari în centrul de control, sintetizîndu-se şi calculîndu-se volumul total de muncă, aplicîndu-se un factor de utilizare S-a apreciat că în centrul de control este nevoie de n persoane, unde n a fost calculat după cum urmează: n = 2 + 0,2 x (numărul funcţionarilor de la ghişeele obişnuite) (10 1) Fină acum studiul s-a bazat pe fapte Următoarea problemă pusă ţinea mai mult de imaginaţie : ce fel de disciplină a aşteptării se va obţine în cazul unui sistem de verificare obişnuit ? Răspunsul era greu de dat După cum vă amintiţi, cele două ghişee speciale de verificare au funcţionat ca şi cum s-ar fi aplicat regula FIFO Părea valabil ca pasagerii să se răspîndească la ghişee, deplasîndu-se liber de la unul la altul, pentru a-şi reduce timpul de aşteptare Deşi nu se respectă strict disciplina clasică de aşteptare FIFO: clienţii sînt serviţi în ordinea sosirii lor, totuşi s-a apreciat că se întîmplă aproximativ acest lucru Standardul de servire luat drept bază de calcul era un timp mediu de aşteptare care nu este influenţat de disciplina de aşteptare în modelele particulare cu mai multe staţii de servire care păreau că aproximează sistemul real de verificare ; cercetătorii au fost de acord să admită regula FIFO ca regulă a disciplinei de aşteptare Astfel, ca model pentru sistemul de verificare la ghişeele obişnuite a fost stabilit M/M/c : (oo/FIFO) cu p = 0,625 şi X (în orele de aglomeraţie) iniţial 3,1 pe minut, momentul de transfer fiind cu 40 min înainte de plecare Graficele lui Molina arată că pentru a ajunge la un timp mediu de aşteptare de 1,0 min în aceste condiţii, numărul de angajaţi trebuie să fie c = 6 în mod similar, c poate fi calculat pentru alte valori ale lui X care corespund altor ore ale zilei Numărul ghişeelor de zbor, numite „ghişee finale de verificare”, a fost stabilit pentru fiecare oră, determinîndu-se numărul de zboruri maxim care pleacă în timpul fiecărui interval de 15 min (fiecare ghişeu final trebuie să fie deschis timp de 15 min) şi repartizîndu-se un angajat dacă avionul este de tipul DC—3, 1,5 pentru un avion de tipul Elisabethan sau Viscount, 2 pentru avion de tipul DC—6B Astfel a fost calculat numărul personalului de control Necesarul de angajaţi şi ghişee astfel calculat a fost ceva mai mic decît cel necesar pentru sistemul vechi de verificare a zborurilor ; şi astfel personalul a început să fie sigur că se va adopta noul sistem de verificare cît mai curînd posibil în acest punct al studiului s-a întîmplat ca cinci membri ai personalului de verificare să se îmbolnăvească cam în acelaşi timp într-o perioadă foarte aglomerată a anului Supraveghetorii au fost de acord că este imposibil să se prelucreze fluxul de călători cu vechiul sistem de verificare, fără personal suficient, şi au hotărît să încerce noul sistem de verificare : un sistem fără ghişee finale Principalele probleme ce apăreau au fost următoarele : timpul era prea scurt pentru instruirea personalului rămas şi lipsea un sistem de transport al taloanelor într-un punct central de control şi sortare; ultima problemă a fost rezolvată ad-hoc, folosindu-se drept curieri angajaţii tineri care duceau taloanele la centrul de control format în grabă Noul sistem de verificare şi prototip au intrat în funcţiune cu două zile înainte de un week-end aglomerat specific englezesc La început totul a decurs conform planului, dar pe măsură ce a crescut traficul, au apărut deficienţe în serie Timpul de aşteptare al clienţilor era cu mult mai mare decît cel anticipat Principala problemă privea influenţa sistemului asupra punctualităţii plecării avioanelor Mulţi pasageri care soseau cu 5 — 6 min înainte de închiderea zborului nu aveau timp să ajungă la ghişee în timp util Pe bună dreptate aceştia reclamau, blocînd ghişeele şi întîrziindu-i şi pe pasagerii din spatele lor, astfel ajungîndu-se la întîrzierea plecării mult mai multor zboruri decît în mod obişnuit Apariţia acestor blocări a fost observată de aproximativ 20 de ori, din care una a durat chiar 18 min în a patra zi a experienţei, supraveghetorii au iniţiat anunţarea pasagerilor cu cinci min înainte de momentul de închidere a fiecărui zbor, astfel încît cei care nu au fost verificaţi încă să treacă imediat în faţă, cerînd prioritate de servire Această măsură nu a avut rezultate semnificative pentru că englezii, care sînt foarte disciplinaţi, nu au vrut să treacă în faţa altora, şi în cazul în care au făcut-o s-a produs dezordine : călcîiul lui Ahile s-a dovedit a fi problema bagajelor Importanţa acestui experiment rezidă în faptul că a scos la iveală problema călătorilor care sosesc tîrziu şi a celor cu probleme speciale De asemenea s-a stabilit că momentul la care trebuie să se fi încheiat verificarea ultimului pasager (momentul de completare a zborului) nu trebuie să depăşească cu cinci min momentul de închidere, altfel în 25% din cazuri plecările avioanelor vor fi întîrziate Cheltuielile şi neajunsurile provocate de aceste plecări întîrziate companiei aeriene sînt aşa de mari încît conducerea acestor companii nu acceptă sisteme de verificare care pot să le favorizeze în acelaşi timp, conducerea (care a fost informată despre gradul de aglomeraţie constatat în orele de vîrf) a fost de acord să îmbunătăţească standardul de servire la ghişeele obişnuite; de acum înainte, timpul mediu de aşteptare al călătorilor va trebui să fie cel mult 30 s în orele de aglomeraţie în acest moment era necesar să se redefinească problema sistemului de verificare optim Acest sistem va fi cel care va necesita numărul minim de salariaţi astfel încît: 1) să se asigure un timp mediu de aşteptare de 0,5 min în orele de vîrf; 2) să se asigure ca cel mult pentru un zbor din 20 momentul de completare să fie cu 5 min mai tîrziu decît momentul de închidere în aceste condiţii stabilirea unui sistem de verificare optim este acelaşi lucru cu stabilirea timpului de transfer optim pentru un sistem compus Aceasta a reprezentat baza pentru o nouă serie de studii 10 2 FAZA A DOUA Cercetătorii care au întreprins această a doua serie de studii aveau experienţă mai multă şi erau mai multilaterali Ei au întreprins o examinare metodică a celor patru componente ale unui sistem de verificare compus : Verificarea pasagerilor la ghişeele obişnuite Aceasta cuprinde : cîntărirea şi etichetarea bagajelor, înscrierea greutăţii bagajelor pe taloanele de zbor, întocmirea unor chitanţe pentru eventuale plăţi suplimentare, emiterea de permise de trecere şi transmiterea taloanelor de zbor spre camera de control Verificarea călătorilor la ghişeele finale cuprinde toate operaţiile îndeplinite la un ghişeu obişnuit, plus totalizarea numărului de călători, plus a greutăţii bagajelor, împreună cu alte cîteva operaţii minore Sortarea taloanelor care sosesc pe bandă de la ghişeele obişnuite, in funcţie de indicativul zborului Controlul, adunarea şi tipărirea în mod automat a numărului călătorilor şi greutăţii bagajelor acestora, confruntarea taloanelor cu listele de călători etc Cercetătorii au hotărît să studieze acest proces complex în anumite condiţii specifice de trafic şi anume: în timpul unei ore de aglomeraţie se va programa ca plecarea zborurilor să aibă loc din cinci în cinci minute Pentru fiecare zbor în parte, 20 de grupuri de pasageri vor sosi înainte de momentul de închidere care este cu 25 minute înainte de momentul plecării Trebuie remarcat faptul că numărul de pasageri dintr-un grup variază, după cum se indică şi în fig 10 2, şi că verificarea pentru zbor nu se face numai la aeroport, deoarece unii pasageri preferă să fie verificaţi la sediul din oraş 10 2 1 VERIFICAREA LA GHIŞEELE COMUNE Modelul MIM Ic : (oo/FIFO) al sistemului de aşteptare la ghişeele obişnuite nu a fost infirmat în practică Totuşi, în timpul care s-a scurs de la studiile iniţiale, la un alt aeroport englez a intrat în funcţiune un sistem compus Din nefericire, proiectul conceptual şi fizic al acestui sistem a fost atît de particular, încît experienţa dobîndită aici nu a putut fi extinsă, cu o singură excepţie, anume repartiţia timpului de servire Activităţile executate la un ghişeu obişnuit la aeroport erau identice cu cele executate la aeroportul din Londra Ca urmare repartiţia prezentată în fig 10 5 a înlocuit repartiţia teoretică folosită anterior Dificultatea pe care o pune această repartiţie constă în faptul că se ştia prea puţin despre modelele de aşteptare cu mai multe staţii de servire avînd repartiţii ale timpului de servire de tip Pearson III Cercetătorii, cel puţin în Anglia, sînt tentaţi să facă următorul raţionament în cazul unui sistem cu o singură staţie de servire ştim că timpul mediu de aşteptare produs de o repartiţie a timpului de servire de tip exponenţial negativă este mai mare decît cel produs prin repartiţii de dispersii mai mici, toate celelalte condiţii fiind aceleaşi Formula lui Pollaczek asigură aceasta Este logic să se presupună că l-’ig 10 5 Repartiţia timpilor de verificare la ghişeele obişnuite Dreptunghiurile reprezintă repartiţia observată Curba de tip Peareon III - — - - Curba exponenţială negativă l-’ig 10 5 Repartiţia timpilor de verificare la ghişeele obişnuite Dreptunghiurile reprezintă repartiţia observată Curba de tip Peareon III - — - - Curba exponenţială negativă lucrurile se întîmplă la fel şi în cazul mai multor staţii de servire Dacă calculăm totul ca şi cum timpul de servire ar fi exponenţial, acţionăm cu prudenţă şi conservatorism Un membru al echipei de cercetători englezi a construit o simulare modestă şi a făcut-o să funcţioneze manual (nu avea la dispoziţie un calculator) pentru a stabili timpul mediu de aşteptare care corespunde timpului de servire observat El a comparat aceste rezultate cu cele obţinute din modelul M/M/c :(oo/FIFO) cu ajutorul tabelelor lui Molina Comparaţia s-a dovedit a fi atît de interesantă, îneît rezultatele ulterioare au fost calculate, folosind diferitele repartiţii Pearson III ale timpului de servire care au fost observate Concluzia a fost că rezultatele foarte bune se pot obţine repede calculîndu-se timpul mediu de aşteptare pentru c staţii de servire, fără a recurge la simulare, după formula (10 2) (10 2)(1 + D2), unde J/Ic)este timpul mediu de aşteptare în cazul c staţii de servire dacă repartiţia timpului de servire ar fi exponenţială negativă cu aceeaşi medie, şi D este dispersia repartiţiei Pearson de tipul III, Această formulă a fost iniţial obţinută de Pollaczek pentru cazul unei singure staţii de servire (c = 1) Concluzia la care s-a ajuns are o mare valoare deoarece este uşor de calculat Mic} din formula lui Erlang c • c !(1 — p)2 (i c • c !(1 — p)2 (i (10 3) (10 3) unde 1/p este timpul mediu de servire iar (10 4) (10 4) i = 'j,1 (cp)r + (cp)c 1 P(0) r! c! 1 - p ’ după cum cititorul îşi aminteşte din cap 3 David Owen a sugerat o aproximare a timpului mediu de aşteptare pentru un model din familia :(oo/FIFO) în tabelul 10 3 sînt comparate rezultatele obţinute, folosindu-se această formulă de aproximaţie cu cele obţinute prin simularea Monte Carlo Cele două serii de rezultate sînt în concordanţă Este clar că pentru repartiţiile Pearson III ale timpului de servire care se întîlnesc în practică în sistemele de verificare (pentru care de obicei 0,70 2) este suficient de 2 corectă Această metodă de aproximaţie a rezolvat principala dificultate a calculelor şi s-a putut calcula numărul ghişeelor, adică al staţiilor de servire necesare unui sistem compus de verificare pentru a prelucra fluxuri diferite de pasageri, asigurînd o servire standard echivalentă cu un timp mediu de aşteptare de pasager de 0,5 min în orele de vîrf Cercetătorii au formalizat apoi estimaţia ritmului mediu de sosire a grupurilor de călători la ghişeele obişnuite Presupunînd că sînt programate să plece în fiecare oră k zboruri la intervale de 60/fc min, ritmul mediu de sosire trebuie să fie (10 5) (10 5)F (g0, t) unde «0 este momentul de transfer şi n, numărul grupurilor de pasageri cu al i-lea zbor F(z, t) este proporţia cumulată a sosirilor cu z minute înainte de închiderea zborului, care are loc cu t minute înainte Tabelul 10 3 Comparaţie între diferite evaluări ale timpului mediu de aşteptare Nr de rtatll de servire Cazul A Cazul B Cazul C 6 = 2 50 Formula de aproximare -0 = 0,75 Simulare 6 = 1 54 Formula de aproximare P = 0 85 Simulare 6=1 59 Formula de aproximare 0=0 96 Simulare 7 2,59 2,25 1,70 1,64 8 — — 0,48 (0,27) 0,40 (0,51) 9 — — 0,17 0,15 0,17 0,12 10 21,80 (8,4) 0,07 0,07 0,06 0,05 11 1,09 0,92 0,03 0,04 0,03 0,03 12 0,58 0,45 — — — — 13 0,16 0,16 — — — — 14 0,07 0,07 — — — — de închiderea zborului, care are loc cu t minute înainte de plecarea avionului Pentru valorile particulare date parametrilor modelului examinat, numărul ghişeelor de verificare necesare, corespunzătoare diferitelor momente de transfer, a fost găsit ca în tabelul 10 4 Tabelul 10 I Numărul angajaţilor necesari pentru ghişeele de verificare primite Timp de transfer Înainte de Închidere, In min i 5 10 15 20 25 30 c ! 8 8 8 7 6 5 La ghişeele de verificare obişnuite condiţiile legate de momentul completării unui zbor nu sînt importante decît dacă z0 — t = 0 (adică dacă nu există ghişee finale) sau dacă probabilitatea ca un pasager să aştepte mai mult de z0 — t este considerabilă O altă restricţie de care trebuie să se ţină seama : riscul ca un pasager care a sosit înainte de închiderea zborului să ajungă prea tîrziu la ghişeu Reformulînd aceste două observaţii prin luarea în consideraţie a unor limitări ale riscului, ele sună astfel: Riscul ca un pasager care a sosit înainte de momentul de transfer să nu poată fi verificat înainte de închidere trebuie să fie atît de mic, încît să apară în cel mult 1 din 20 de zboruri Situaţia unei verificări a unui pasager începută înainte de închidere şi care nu s-a terminat după 5 min de la închidere trebuie să apară în cel mult 1 din 20 de zboruri Examinînd prima restricţie, s-a constatat că calculul riscului situaţiilor de acest tip implică cunoaşterea repartiţiei timpului de aşteptare Pe cînd utilizarea modelului Jf/Jf/c :(oo/FIFO) apare justificată pentru calculul timpului mediu de aşteptare mai ales după corecţiile realizate folosindu-se ipoteza lui Owen, cercetătorii nu au avut tot atîta încredere că se vor obţine estimaţii corecte ale probabilităţilor asociate părţii superioare a repartiţiei timpului de aşteptare Această dificultate a fost greu de rezolvat, deoarece se părea că nu era nici un mod de a se deduce repartiţia timpului de aşteptare pentru un model mai realist Simulările Monte Carlo făcute pentru un sistem obişnuit complet cu z0 — t = 0 au generat totuşi repartiţii ale timpului de aşteptare, care sînt prezentate în fig 10 6 pentru diverse numere de ghişee obişnuite Este evident că într-un sistem compus în care momentul de transfer este zQ, probabilitatea q{z0, c; v) a unui timp de aşteptare ce este mai mare decît v în condiţiile în care există c angajaţi, trebuie să satisfacă condiţia c; v) 0,5} 0,05 Aceasta este acelaşi lucru cu aflarea celei mai mari valori a lui p care satisface relaţia în care apare o funcţie beta incompletă : IpţlO, 21) (C^1)N - • (13 4) Transformata Laplace a acestei funcţii este ,f*(z)=(—r~e+i (13 5) \ Cjx + z / Modelul utilizat pentru a doua fază presupune de asemenea încărcări aleatoare Cititorul va identifica cu siguranţă ecuaţiile modelului, care sînt: —p(c - l,t) - (c — l)up(C - 1,0, di — p(r, i) = — rnp(r, t) + (r + 1) np (r + 1, t) (1 r di -^-7>(0, i) = p p(l, i) di Soluţia se găseşte uşor : 70 %, numărul de angajaţi necesari pentru o utilizare mai puţin obositoare O variabilă pe care nu am discutat-o încă este timpul mediu de servire al angajaţilor Acesta poate să difere de la un oficiu la altul, dar a trebuit să fie conţinut în tabele Efectul modificărilor acestei variabile a fost punctul central al cercetărilor efectuate în legătură cu introducerea sistemului electronic de rezervare S-a sugerat că introducerea acestui sistem va afecta timpul necesar unui angajat pentru a trata cu clientul Aceasta ar fi putut fi adevărat, astfel că problema a fost cercetată iar rezultatele acestor cercetări, realizate la un oficiu cu tablou de distribuţie automat, sînt prezentate în tabelul 15 3 Tabelul 15 3 Efectul variaţiei timpilor medii de servire asupra numărului de linii şi angajaţi Timpul mediu de servire (min) 1 5 2 3 4 NumArul necesar de linii In orele de vîrf 19 24 33 40 Apeluri pe oră Numărul angajaţilor 200 10 12 15 17 400 20 25 33 41 600 24 33 45 56 Numărul de linii a fost astfel ales, încît procentul de apeluri respinse să nu depăşească standardul în orele de vîrf; acelaşi număr de linii a fost utilizat pentru două frecvenţe mai scăzute ale intrărilor Numărul angajaţilor a fost dat de rezultatele simulărilor, ajustate pentru a înlătura o utilizare mai mare de 70% COMPARAŢII CU TEORIA AŞTEPTĂRII Whitton a tras concluzia, corectă dealtfel, că o continuare ce constă dintr-un număr mare de cicluri de simulare, producînd volume imense de date, nu era cea mai bună soluţie El a considerat că dacă se putea găsi o cale de utilizare a modelelor teoretice existente pentru a produce aceleaşi rezultate ar fi mult mai bine chiar dacă de fapt aceste modele diferă foarte mult în structura lor de ceea ce se întîmplă în oficiile de rezervare de locuri De aceea, el a realizat o comparare a rezultatelor obţinute în cîteva cicluri de simulare cu cele rezultate din unele modele matematice Vom continua prezentarea acestui stadiu al primei faze cu propriile sale cuvinte „Cel mai simplu sistem de aşteptare cuprins în sistem este unul în care apelurile sosesc aleator şi sînt tratate în ordinea sosirii de către oricare dintre angajaţi, dimensiunea şirului de aşteptare fiind limitată de numărul fix de linii telefonice Dacă timpii de servire sînt presupuşi a fi repartizaţi exponenţial, atunci pot fi obţinute valorile teoretice pentru diferite statistici utile, prin metode standard Valorile deduse dintr-un astfel de sistem vor fi valabile numai pentru tipul cel mai simplu de oficiu, în care apelurile ajung direct la angajaţi, fără a trece printr-o fază intermediară de aşteptare la tabloul de distribuţie Cele mai multe oficii însă au tablouri de distribuţie, astfel că cea mai evidentă generalizare a unui sistem simplu constă în introducerea unui al doilea stadiu de servire, adică următorul sistem ce trebuie luat în consideraţie este un sistem ce conţine două faze în serie, fiecare cu mai multe staţii de servire, cu o limitare superioară a numărului de apeluri acceptate în sistem Formulele aplicabile unui astfel de sistem nu pot fi găsite în literatură, dar dificultăţile legate de deducerea lor sînt numai algebrice Cea mai facilă ilustrare a modului în care expresiile teoretice corespund rezultatelor simulării constă într-o comparaţie a estimaţiilor procentului de apeluri respinse O astfel de comparaţie poate fi făcută cu ajutorul următorului tabel, care dă estimaţiile acestor procente pentru o frecvenţă a intrărilor de 600 apeluri pe oră, cu 40 de linii telefonice şi un timp mediu de servire de 3 min pentru sistemele cu o singură fază şi respectiv timpii medii de servire de 6 s şi respectiv 3 min pentru sistemele formate din două faze (tabelul 15 4) Diferen- Tabelul 15 4 Comparaţie a estimaţiilor teoretice ale procentului de apeluri respinse, cu estimaţiile produse prin simulare NumArul fazelor Numărul angalatllor NumArul fazelor Numărul angalatllor 1 1 20 i 40 • Prima fază A doua fază 2 20 2 2 40 1 1 20 i 40 • Prima fază A doua fază 2 20 2 2 40 Estimarea teoretlcA Estimarea prin simulare (media a trei cicluri) 33,3 36,9 1,4 1,8 33,3 38,6 2,2 2,3 Estimarea teoretlcA Estimarea prin simulare (media a trei cicluri) 33,3 36,9 1,4 1,8 33,3 38,6 2,2 2,3 ţele ce apar în cele două cazuri cu 40 de linii şi 20 de angajaţi sînt prea mari pentru a fi atribuite în întregime variaţiilor aleatoare : acestea ar putea fi cauzate fie de caracterul nexeponenţial al repartiţiilor timpului de servire, fie de caracterul prea grosolan al teoriei Repartiţiile teoretice ale timpului de aşteptare pentru cazul cu o singură fază au aceeaşi formă cu cele produse prin simulare (adică acelaşi grad de asimetrie) dar valorile medii ale celor două tipuri de repartiţie diferă Şi aici probabil că o îmbunătăţire a teoriei ar putea furniza aproximaţii mai bune” SFÎRŞITUL PRIMEI FAZE în acest moment, Compania Air Canada posedă o simulare care va arăta, pentru diferite frecvenţe de intrare repartiţii ale timpului de servire şi regula FIFO a disciplinei de aşteptare, numărul de linii telefonice şi de angajaţi, procentul de apeluri respinse şi procentul de apeluri acceptate, care trebuie însă să aştepte mai mult decît un timp specificat Simularea a fost programată pentru un calculator IBM 7070 în 1962 au fost simulate un număr mare de cicluri pentru a reprezenta unul din cele mai mari oficii de rezervări de locuri ale companiei Acest sistem utilizează un tablou de distribuţie automată, cu acces direct, şi este completat cu circuite automate care asigură o perfectă regulă FIFO a disciplinei de aşteptare în urma ciclurilor de simulare s-au completat cîteva tabele Imediat după ce aceste tabele au fost completate, s-a efectuat o reorganizare în cadrul companiei Scopul acestor modificări a fost realizarea unei mai bune funcţionări în faţa cererii crescînde Ca rezultat, Whitton, a cărui experienţă pentru sistemele de mari dimensiuni era considerabilă, a primit noi responsabilităţi, iar legătura sa cu proiectul pe care-1 discutăm a fost întreruptă 15 2 FAZA A DOUA Sarcina continuării proiectului şi a colaborării cu departamentul de analiză a încasărilor i-a revenit luiR W Binder şi autorului Este, cred, un fapt incontestabil — deşi regretabil şi neconvenabil — că realizarea unui proiect de cercetare operaţională depinde mai curînd de persoanele angajate direct decît de tehnicile care se pot utiliza, în cazul prezentat, era inevitabil ca modificările ce au survenit în echipa de studiu să influenţeze modul de abordare a problemei şi într-adevăr, aşa a şi fost împreună cu Linder ne-am axat pe două probleme : Standardele de servire pe care compania doreşte să le impună ca minime în întregul sistem Principalul rezultat al proiectului, anume calculul necesarului de forţă de muncă, cerut de modificările survenite, pentru oficiile de rezervare de locuri Prin aceasta, noi am privit lucrurile în cu totul alt mod Standardele minime ale companiei pentru rezolvarea apelurilor telefonice erau : să nu fie respinse mai mult de 15% din apeluri, în orice oră; să nu fie mai mult de 18% apeluri acceptate, în orice oră, care să stea în aşteptare mai mult de 20 s, înainte de a fi preluate de către un angajat Tabelele obţinute prin simulare puteau fi imediat aplicate, pentru a obţine numărul de linii şi angajaţi necesar a satisface aceste condiţii, la diferite nivele de trafic Exemple de astfel de tabele sînt reprezentate în tabelele 15 5 şi 15 6 Să analizăm aceste tabele, în scopul identificării unor reguli de lucru, uşor de aplicat de către conducătorii locali şi în acelaşi timp să găsim o metodă aproximativă pe care să o putem utiliza în alte studii, în locul unor cicluri de simulare costisitoare După o scurtă căutare, am stabilit că numărul minim de linii şi de angajaţi necesari asigurării standardului minim al companiei, poate fi calculat astfel: Numărul necesar de angajaţi în fiecare oră este E + 3, unde E este numărul de erlangi ai traficului acceptat în ora respectivă, sau 85% din traficul oferit Numărul necesar de linii este E + 5 Tabelul 15 7 are ca scop compararea numărului de linii şi de angajaţi necesari, rezultaţi în urma simulării şi în urma aplicării regulii „trei plus cinci”, pentru a face faţă volumelor de trafic oferit, la nivelul standardelor minime ale companiei Comparaţia realizată arată abateri mici pentru valori scăzute ale nivelelor de trafic şi abateri nesemnificative pentru nivele ridicate Ni s-a părut rezonabil să recomandăm utilizarea acestei reguli de aproximare de către conducătorii locali, în limita condiţiilor precizate şi pentru valori ale traficului dintre 5 şi 25 de erlangi (în afara acestor limite, se pare că regula nu este bună) Aceasta era scopul simulării: obţinerea unei reguli de aproximare pentru calculul necesarului de linii şi de angajaţi în realitate, anumite oficii prezintă standarde de servire mai ridicate decît cele minime impuse de companie (un fenomen des întîlnit în transporturile aeriene); de aceea pentru aceste oficii este necesară determinarea unei alte reguli (aceasta a fost regula „patru plus şapte”) Fără îndoială, că prin rezolvarea acestor probleme s-a făcut un pas important în realizarea proiectului Procentul de apeluri respinse Frecvenţa: 250 apeluri/oră Timp de servire: 4 min Nr de angajaţi 22 21 20 19 18 17 16 15 14 13 12 11 0 0 0 2 6 10 14 18 23 27 32 39 0 0 0 2 6 10 14 18 23 27 32 39 0 0 0 2 6 10 13 17 22 27 32 39 0 0 0 2 5 9 13 16 22 26 31 38 0 0 0 2 5 9 12 16 21 26 31 38 0 0 0 2 5 9 12 16 21 26 31 38 0 0 0 2 5 9 12 15 21 26 31 38 0 0 0 1 4 8 11 15 20 25 31 38 0 0 0 1 4 8 11 15 20 25 31 38 Nr de Huli 18 19 20 21 22 23 24 25 26 Frecvenţa : 400 apeluri/oră Timp de servire : 2,5 min Nr de angajaţi 22 21 20 19 18 17 16 15 14 13 12 11 12 13 14 15 16 17 18 20 23 25 32 36 9 10 12 13 14 15 17 19 23 25 32 36 7 8 10 11 12 14 15 18 22 25 31 36 4 6 8 9 10 12 13 17 21 24 31 36 2 3 6 7 10 12 16 20 24 31 36 2 3 5 6 7 9 12 15 20 24 31 36 2 2 4 5 6 9 11 15 20 24 30 35 1 2 3 4 6 8 11 14 19 23 30 35 1 1 2 3 5 8 11 14 19 23 30 35 Nr de linii 18 19 20 21 22 23 24 25 26 Frecvenţa : 500 apeluri/oră Timpul de servire : 2 min Nr de angajaţi 22 6 5 5 4 4 4 3 2 1 21 7 7 6 5 5 5 4 3 2 20 9 8 7 6 6 5 5 4 3 19 10 9 8 7 7 6 5 4 4 Nr de angajaţi 18 11 10 9 8 8 7 7 6 •G 17 13 12 11 10 10 9 9 8 8 16 15 15 14 13 13 12 12 11 11 15 19 18 17 16 15 15 15 15 15 14 22 21 21 20 20 20 20 20 20 13 26 25 24 23 23 23 23 23 ■ 22 12 32 32 31 31 31 31 31 31 31 11 38 38 38 37 37 37 37 37 : 37 Nr de linii | 18 | 19 | 20 | 21 j 22 | 23 | 21 25 { ' 26 Principala problemă era următoarea : ca să fie capabil să folosească aceste reguli, un conducător trebuia să poată să procure sau să est imeze volumul traficului pentru viitor Cu alte cuvinte, el trebuia să prog- nozeze valori ale intensităţii intrărilor şi ale timpilor medii, de servire Acestea din urmă sînt constante în timp Frecvenţele intrărilor însă sînt fluctuante în timpul zilei, săptămînal şi anual Trebuia deci să rezolvăm într-un fel această problemă După o examinare a înregistrărilor anterioare, am dezvoltat o metodă de previzionare a frecvenţelor orare ale apelurilor, deduse din previzionarea numărului de pasageri Aceasta utiliza o relaţie între frecvenţele apelurilor la un oficiu de rezervări dintr-un anumit oraş şi numărul pasagerilor ce utilizau cursele aeriene S-au obţinut astfel relaţii foarte utile O altă problemă pe care am întîlnit-o în această fază a studiului a fost aceea de estimare a numărului de clienţi ce cheamă telefonic oficiul de rezervări, în condiţiile în care singurele date înregistrate erau numărul de apeluri acceptate şi numărul celor respinse O scurtă reflecţie a cititorului îl va convinge că prin adunarea acestora nu vom obţine numărul total de apeluri Un client, al cărui apel este respins prima oară, în general nu renunţă Aşa cum făcea d-na Overlander, el mai încearcă de cîteva ori Fiecare încercare nereuşită se adună la numărul apelurilor respinse Considerînd acceptarea sau respingerea apelurilor la prima, la a doua şi la următoarele încercări, a fost posibil să obţinem o estimaţie brută a cererii adevărate, în fiecare oră S-a presupus că un apel se Procentul de apeluri acceptate, intirzlate mal mult de 20 s Frecvenţa : 250 apeluri/oră Timpul de servire : 4 min Nr :de angajaţi 22 0 1 1 2 2 2 2 2 2 21 1 2 2 2 3 3 3 4 4 20 2 2 2 3 4 5 7 8 19 3 4 5 6 8 11 14 17 20 ? 18 6 8 10 13 15 20 25 30 36 17 14 17 20 23 26 37 48 63 79 16 26 32 38 44 50 60 70 80 91 15 i 36 61 65 70 74 79 84 89 94 14 75 79 83 87 91 92 94 96 97 13 89 91 93 95 97 97 98 98 99 12 95 96 97 98 99 99 99 99 99 11 1 98 98 98 99 99 99 99 99 99 Nr de linii j ÎS 1 19 j 20 1 21 1 22 | 23 | 24 | 25 26 Frecvenţa : 400 apeluri/oră Timpul de servire 2,5 min Nr de angajat! ' 22 0 0 1 1 2 4 5 7 8 21 1 2 2 3 3 5 8 11 13 20 2 4 5 6 7 11 14 18 21 19 3 6 8 11 14 19 24 29 34 18 6 11 16 21 27 33 39 44 49 17 12 20 29 37 46 50 55 60 64 16 ' 25 34 43 52 62 66 70 74 78 15 50 56 63 70 76 80 83 87 91 14 ! 67 73 78 84 89 91 93 95 96 13 84 86 89 92 95 96 97 97 98 12 92 93 95 97 98 98 99 99 99 11 98 98 98 99 99 99 99 99 99 Nr de linii 1 18 19 20 21 22 23 24 25 26 Frecvenţa : 500 apeluri/oră Timpul de servire : 2 min Nr de angajaţi I 22 0 0 1 1 2 3 5 6 8 21 0 1 2 4 5 7 9 11 13 20 0 2 4 6 8 11 14 18 21 ' 19 1 4 7 11 14 19 24 29 34 N’ 0) • Q( > R), unde Wî( > 0) şi Q( > R) sînt funcţiile prezentate în graficile A3-6 şi respectiv A3-7 FORMULE UTILE Această parte este destinată exclusiv prezentării unor formule foarte utilizate in teoria aşteptării Pentru anumite modele exiBtă formule generale, dar ele nu sînt prezentate aici, deoarece sînt dificil de înţeles şi de calculat 1 O staţie de servire, sosiri repartizate Poisson M[G/l: (oo/GD) Formulele lui Pollaczek — Hincin pot fi utilizate pentru estimarea numărului mediu de clienţi existenţi in sistem şi a timpului mediu de aşteptare Alte elemente nu se obţin cu uşurinţă , exceptînd cazurile particulare în care repartiţia timpului de servire este exponenţială, constantă sau erlangiană : n=p+pt(1+rg), Wt = P(1 + 2(1 - p) ’ 2(1 - p)|i Jf/Jf/1 : (oo/GD) Repartiţiile numărului de clienţi în aşteptare (în şirul de aşteptare sau în sistem) nu sînt dependente de disciplina de aşteptare Soluţia în starea de tranziţie este p(n, 0 = exp { - (X + n)0 x X p" 2ZB(o;)4-p(n 1)/2-fn+i(a;)+p"(l—p) Ş p-'"+'>'2Z„+,(*)] , unde x = 2t /xjî Soluţia în starea stabilă este P» = (1 p) p", îl = 0,1, 2, Timpul total mediu de aşteptare şi timpul mediu de aşteptare nu sînt afectaţi de natura disciplinei de aşteptare Ei sînt Repartiţiile timpilor de aşteptare şi a timpilor totali de aşteptare sînt afectaţi de disciplina de aşteptare Avem : ® = 0, v >0, ® = 0, v >0,-FIFO w4(«) = w4(«) =1 ~P> p(l —p);ie- 0 - SIRO Probabilitatea ca un client care trebuie să aştepte mai mult decît t este dată de aproximarea lui Riordan (v graficul A3-7): exp {— a(l - p) (X t} + -| exp {- P(1 - p) [Lt}, 2 2 unde a = 1 + Kp/2, iar p = 1 — /p/2 - LIFO Repartiţia timpilor totali de aşteptare este identică cu repartiţia perioadelor ocupate pentru modelul FIFO Densitatea de repartiţie este w(u) = exp {— (1 — p) [x v} -1^— M/Ek/1 : (oo/FIFO) Repartiţia numărului de clienţi existenţi în sistem este foarte uşor exprimată cu ajutorul funcţiei generatoare; din aceasta diferitele valori ale lui pn pot fi obţinute prin derivare repetată : J- [9 1 — S — ps(l — s)* Transformata Laplace a repartiţiei timpilor de aşteptare este dată de (1 + O tratare analitică a modelului este destul de anevoioasă Exceptînd cazul k = 1, calculul inversiei transformatei Laplace necesită anumite funcţii speciale, studiate de Luchak, care se aseamănă într-o oarecare măsură cu funcţiile Bessel incomplete Repartiţia perioadei ocupate ridică aceleaşi probleme şi de aceea nu există o formulă utilizabilă pentru modelul MIEk/l : (oo/LIFO) M/D/l :(oo/FIFO) Probabilităţile în starea stabilă sînt date de Repartiţia timpilor de aşteptare este dată de We(®) = (1 — p) y exp rp(ix v — i)l —P^ — Pentru a obţine mediile, se pot utiliza formulele lui Pollaczek — Hincin Valorile medii cele mai utilizate sînt 1 ,u 1 ,u M/D/l : (oo/SIRO) Formulele acestui model nu s-au obţinut încă; vezi graficele A3-1 - A3-4 SERVICIU SERVICIU PRIOR ITAR PRIOR ITAR Singurele formule utilizabile sînt acelea ale mediilor Pentru două clase de priorităţi relative avem : "1 "1 qp(l +p — ap) qp(l +p — ap) 1 — qp 1 — qp (1 — p)(l — ap) ’ (1 — p)(l — ap) ’ W91 = W91 = ap2 1 — qp ap2 1 — qp «î2= «î2= (1 ~ «) P2 (1 - p)(l - ap) ’ (1 ~ «) P2 (1 - p)(l - ap) ’ X n(H — °cX) ’ X n(H — °cX) ’ Wg2 = (H —X) (h — aX) Wg2 = (H —X) (h — aX) Pentru două clase de prioritate, din care una are prioritate absolută faţă de cealaltă, avem : ’ Hi 7 \ Hi 7 ’ Hi 7 \ Hi 7 ni — 0,1, 2, , ni — 0,1, 2, , n n M— P1 ni — > 1 — Pi unde pj = oqX/p !, V,-(l +-i—IA, V 1 — Pl ) Hi V 1—Pl/ Hi W — Pi w rrtfl— „ , 5 rrff2 Pa(l — pi) MAI MULTE STAŢII DE SERVERE, SOSIRI REPARTIZATE POISSON , M/M/c: (oo/FIFO) Pentru acest model există grafice foarte utile Graficele A3-5 şi A3-6 sînt ambele utilizabile în plus, există graficele lui Molina, mult mai detaliate, ce pot fi găsite în lucrarea : Molina, E C Application of the Theory of Probability to Telephone Trunking Problems în : Beli System Technical Journal, July, 1927 De aceea formulele nu sînt de mare folos practicienilor Principalele formule sînt : Pi = Pi = Cc — pnp0, n>c c! Cc — pnp0, n>c c!1 y («P)r , (cp)c 1 ; ,“0 r! c! 1 - P Probabilitatea ca un client să sosească şi să găsească toate staţiile de servire ocupate este 00 XP,= r-c 00 XP,= r-c (cp)c c!(l-p) (cp)c c!(l-p) Po" Po" Această cantitate este frecvent notată în literatura prin simbolul P( > 0) în carte ea este notată $(0) şi este egală cu complementul lui Ws(0), care este probabilitatea unui timp de aşteptare pozitiv Similar Q(v) este probabilitatea unui timp de aşteptare ce depăşeşte v, care este complementul lui Ws(®) Această cantitate este notată de obicei prin P( > t) Atunci Q(v) = exp [— cp, (1 - p) ®] • Q(0) P(cp)c c! (1 -p) P(cp)c c! (1 -p) Poi Poi M/M/c: (Z/FIFO) Nu există grafice utile pentru acest model Principalele formule sînt (cp)" cc „ • — p Po, c Probabili la I cn aşteptării (a ocupării simultane a tuturor slaliilor de servire) pentru 1/(o?/CD) l-'ig A2-7 M/M/c : (oo/SIBO): repartiţia timpilor de aşteptare pentru clienţii nevoili să aştepte Timpul de aşteptare (v/b) în multipli ai timpului (constant)mediu de servire l'ig A2-4 Coinparaţic intre repartiţiile timpilor de aşteptare pentru M/D/l : (co/SIRO) şi M/A//1 : (co/SIRO) 3 4 5 6 8 10 15 20 30 W 60 100 Numărul staţiilor de servire, C 3 4 5 6 8 10 15 20 30 W 60 100 Numărul staţiilor de servire, C Fig A2-5 Timpul mediu de aşteptare pentru M/M/c : (co/SIRO) Fig A2-5 Timpul mediu de aşteptare pentru M/M/c : (co/SIRO) oL=gradul de ocupare o staţiei de servire =Aţpc oL=gradul de ocupare o staţiei de servire =Aţpc l ig A2-2 A//D/: oo/ţFIFO): repartiţia timpului de aşteptare Fig A2-3 Comparaţie intre repartiţiile timpilor de aşteptare corespunzătoare modelelor M/D/l : (oo/SIRO) şi M/D/l : (oo/FIFO) BIBLIOGRAFIE Lucrări citate de către autor: Cox, D R ( Smith, W L Queues London, Methuen, 1960 Saaty T L Elements of Queueing Theory with Applications New York, McGraw- Hill, 1961 Morse, P Queues, Inventories and Maintenance New York, Wiley, 1958 Feller, W An Introduciion to Probability Theory and its Applications New York, Wiley, 1957 Kendall, M G , Stuart, A The Adoanced Theory of Statistics I New York, Hafner, 1958 Widder, D V The Laplace Transform Princeton, Princeton Univ Press, 1949 Tocher, K D The Art of Simulation London, English Univ Press, 1963 Oxenfeldt, A E (Edj Models of Markets New York-London, Columbia Univ Press, 1963 Edie, L C Trafic Delays at Toii Booths în : Op Research, voi 2, br 2, 1954 Friend, J E Two Studies in Airport Congestion în : Oper Res Quarterly, voi 9, nr 3, 1958 Fry, J , Jackson, R R P , Welch, J D Appointment Systems in Hospitals and General Practice în : Oper Res Quarterly, voi 15, nr 3, 1964 Gebhart, R F A Limiting Distribution of an Estimate of Mean Queue Lenght în : Op Research, voi 11, nr 6, 1963 Lee, A M Some Aspects of a Control and Communication Systems în : Oper Res Quarterly, voi 10, nr 4, 1959 Lee, A M , Linder, R W Two Problems in Scheduling Manpower în: J Canad Oper, Res Soc , voi 3, nr 1, 1965 Lee, A M , Longton, P A Queueing Processes Associated with Airline Passenger Check-in în : Oper Res Quarterly, voi 10, nr 1, 1959 Molina, E D Appltcation of the Theory of Probability to Telephone Trunking Pro blems în: Beli System Teleph J , voi 6, July, 1925 Shelton, J R Solution Methods for Waiting Line Problems în: J Industrial Eng , voi 2, nr 4, 1960 Stillson, P Implementation Problems in O R în : Op Research, voi 11, nr 1, 1963 Tocher, K D The Role of Models in Operaţional Research în : J Royal Stat Soc , Part II, 1961 Wilkinson, R I Working Curves for Delayed Exponenţial Calls Served in Random Order în : Beli System Teleph J , voi 34, 1960 Cititorul poate consulta, pentru explicarea noţiunilor teoretice utilizate, următoarele lucrări în limba română: Mihoc, Gh , Ciucu, G Introducere in teoria aşteptării Editura tehnică, Bucureşti, 1967 Mihoc, Gh , Ciucu, G , Muja, A Modele matematice ale aşteptării Bucureşti, Edi tura Academiei R S R , 1973 Ciucu, G Elemente de teoria probabilităţilor şi statistică matematică Bucureşti, Edi tura didactică şi pedagogică, 1963 Iosifescu, M , Mihoc, Gh ş a Teoria probabilităţilor şi statistica matematică Bucureşti, Editura tehnică, 1966 Mihoc, Gh , Ciucu, G , Craiu, V Teoria probabilităţilor şi statistică matematică Bucureşti, Editura didactică şi pedagogică, 1970 Kaufmann, A Metode şi modele ale cercetării operaţionale Voi I, Bucureşti, Editura ştiinţifică, 1967 CUPRINS Prefaţă Introducere Problema oţelăriilor Problema casei de bilete Problema soţiei furioase Natura sistemelor de aşteptare Modele ale sistemelor de aşteptare Generalităţi Caracteristicile modelelor de aşteptare Modele cu o singură staţie de servire Generalităţi Modele cu timpi de servire exponenţiali Modele cu timpi de servire neesponenţiali Modele cu timpul de servire constant Reguli de priorităţi pentru disciplina de aşteptare Comentariu general Modele cu mai multe staţii de servire Modele Fenomene de deplasare Modele ale sistemelor seriale şi ciclice Sisteme seriale Sisteme ciclice Tehnici de simulare Metode Monte Carlo Numere aleatoare şi pseudoaleatoare Generarea variaţiilor aleatoare Utilizarea simulărilor Experimente de simulare Concluzie Natura soluţiilor operaţionale Trei tipuri de aproximări Sisteme de rezervare de locuri Descrierea problemei Prima fază Interludiu 118 Faza a doua 119 Comentariu 123 Metode de control al călătorilor 124 Prima fază 125 Faza a doua 134 Studiul de validare 147 Comentariu 152 Problema sortării cupoanelor 153 Patru situaţii selectate 158 Proiectarea reţelei de piste în aeroport 158 Problema autobuzelor în aeroport 160 Întîrzieri ale circulaţiei în punctele de vamă 161 Sisteme de intîlnire 164 Proiectarea sediului din londra al companiei bea 166 Prezentarea problemei 166 Capacitatea locului de parcare 168 împrăştierea pasagerilor ce părăsesc sediul 171 încărcarea autobuzelor spre aeroport 180 Probleme de aşteptare într-un sistem de control 186 Prezentare a problemei 187 Studiul preliminar 188 Studiul principal: controlul plecărilor 190 lo O problemă a standardelor de servire 195 Prima fază 196 Faza a doua 211 Organizarea schimburilor 217 Repartizarea personalului 219 Comentariu 221 16 Concluzii 223 Anexa 1 Lista principalelor simboluri 228 Anexa 2 Formule şi grafice utile 230 Bibliografie 245 Redactor: Valentina Creţu Tehnoredactor: Elena Geru Coperta şi supracoperta seriei: Valentin Vişan Bun de tipar: 14 06 1976 Coli de tipar: 15,50 Tiraj: 5850 + 50 exemplare legate C Z : 517 8:510 9 O 182 I P INFORMAŢIA Str Brezolanu 23 — 25, Bucureşti